Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_pgi 3,440,243 T→G T263P (ACT→CCT)  rpoA ← RNA polymerase, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_pgi3,440,2430TG98.8% 299.2 / NA 86T263P (ACT→CCT) rpoARNA polymerase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (41/44);  total (42/44)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACGAATCGCCTC  >  NC_000913_3_pae_pgi/3440097‑3440387
                                                                                                                                                  |                                                                                                                                                
tCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCaa                                                                                                                                               >  1:342138/1‑150 (MQ=255)
   aCGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTc                                                                                                                                            >  2:301361/1‑150 (MQ=255)
   aCGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTc                                                                                                                                            >  1:389884/1‑150 (MQ=255)
       gAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGa                                                                                                                                        <  2:341777/150‑1 (MQ=255)
        aaGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGAt                                                                                                                                       <  2:499010/150‑1 (MQ=255)
         aGCCAGCGCGTCTTTAATCTGAGTAATATATTTTTTACCAAGGTTAGGCGGTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCTGTTAGCAGAGCTGTCAGGCAATTCCAGATc                                                                                                                                      <  2:581199/150‑1 (MQ=255)
          gCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCg                                                                                                                                     <  3:140305/150‑1 (MQ=255)
           ccAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACTAAGGTTATGCGTTTTATTGAGCTCAACCCCGGTACGCTGTACCAGATCACCTATCTAGTGTATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCACATccc                                                                                                                                    >  2:597361/1‑148 (MQ=255)
            cAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTc                                                                                                                                   <  2:144970/150‑1 (MQ=255)
              gCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCaa                                                                                                                                 <  4:636920/150‑1 (MQ=255)
                   tCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATCTAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAggg                                                                                                                            >  1:791660/1‑150 (MQ=255)
                       tAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCgca                                                                                                                        <  1:195633/150‑1 (MQ=255)
                       tAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCgca                                                                                                                        >  2:195633/1‑150 (MQ=255)
                              aGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAggat                                                                                                                 >  1:390084/1‑150 (MQ=255)
                                  agagATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAggatcgga                                                                                                             <  2:111189/150‑1 (MQ=255)
                                   gagaTTATTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAggatcggat                                                                                                            <  1:605911/150‑1 (MQ=255)
                                    agaTTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAggatcggatc                                                                                                           >  2:1090167/1‑150 (MQ=255)
                                    agaTTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAggatcggatc                                                                                                           <  2:1064218/150‑1 (MQ=255)
                                        tttttACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAAc                                                                                                       >  1:962063/1‑150 (MQ=255)
                                              ccAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGt                                                                                                 >  4:398634/1‑150 (MQ=255)
                                                      aGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCtcttc                                                                                         <  3:489489/150‑1 (MQ=255)
                                                      aGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCtcttc                                                                                         <  1:616564/150‑1 (MQ=255)
                                                       ggCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTc                                                                                                                                   <  2:41849/107‑1 (MQ=255)
                                                       ggCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTc                                                                                                                                   >  1:41849/1‑107 (MQ=255)
                                                       ggCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCtcttct                                                                                        >  3:330933/1‑150 (MQ=255)
                                                        gCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTcttctt                                                                                       >  3:448704/1‑150 (MQ=255)
                                                         cctggcaaGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACACGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCttt                                                                                      <  2:714017/144‑1 (MQ=255)
                                                           ttttAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCa                                                                                    <  2:337913/150‑1 (MQ=255)
                                                           ttttAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCa                                                                                    <  3:398634/150‑1 (MQ=255)
                                                                    gCTCAACCCCGGTACGCTGTACGAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTGCAGATCGTCAACAGGGCGTAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCt                                                                           <  2:33526/150‑1 (MQ=255)
                                                                      tCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGa                                                                         >  1:1123577/1‑150 (MQ=255)
                                                                       cAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGAc                                                                        <  4:539209/150‑1 (MQ=255)
                                                                           cccGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTAc                                                                    >  3:32316/3‑150 (MQ=255)
                                                                           cccGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTAc                                                                    <  2:342138/148‑1 (MQ=255)
                                                                                tACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACAt                                                                  >  4:402615/1‑147 (MQ=255)
                                                                                tACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACAt                                                                  <  3:402615/147‑1 (MQ=255)
                                                                                  cGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGt                                                             <  2:1099113/150‑1 (MQ=255)
                                                                                      gTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGt                                                         <  2:1117724/150‑1 (MQ=255)
                                                                                       tACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTc                                                        >  2:818210/1‑150 (MQ=255)
                                                                                             aTCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGaa                                                  >  2:451675/1‑150 (MQ=255)
                                                                                              tCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGaaa                                                 >  2:407868/1‑150 (MQ=255)
                                                                                                  cGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCtt                                             >  3:454674/1‑150 (MQ=255)
                                                                                                     tataGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCa                                          <  4:330933/150‑1 (MQ=255)
                                                                                                         gTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCttt                                                                                      >  1:264809/1‑102 (MQ=255)
                                                                                                         gTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCttt                                                                                      <  2:264809/102‑1 (MQ=255)
                                                                                                          tGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAgttgt                                     >  4:52505/1‑150 (MQ=255)
                                                                                                            gATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTc                                   >  2:900350/1‑150 (MQ=255)
                                                                                                             aTAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTAc                                                                    >  4:602756/1‑116 (MQ=255)
                                                                                                             aTAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTAc                                                                    <  3:602756/116‑1 (MQ=255)
                                                                                                              tAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAggg                                                                                                                            <  1:618099/59‑1 (MQ=255)
                                                                                                              tAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAggg                                                                                                                            >  2:618099/1‑59 (MQ=255)
                                                                                                               aGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGc                                <  1:1090167/150‑1 (MQ=255)
                                                                                                                gCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGcc                               <  1:900350/150‑1 (MQ=255)
                                                                                                                 cTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCa                              <  4:32316/150‑1 (MQ=255)
                                                                                                                  ttCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAg                             <  4:479889/150‑1 (MQ=255)
                                                                                                                  ttCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAACAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACg                                                              <  3:295377/117‑1 (MQ=255)
                                                                                                                  ttCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAACAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACg                                                              >  4:295377/1‑117 (MQ=255)
                                                                                                                   tCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCaa                                                                                                                                               >  2:197600/1‑35 (MQ=255)
                                                                                                                   tCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCaa                                                                                                                                               <  1:197600/35‑1 (MQ=255)
                                                                                                                   tCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCTAACTCTGGTTTCTCTTCTTTCACTTCa                                                                               <  1:817677/99‑1 (MQ=255)
                                                                                                                   tCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCa                                                                               >  2:817677/1‑99 (MQ=255)
                                                                                                                   tCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGa                            >  1:378325/1‑150 (MQ=255)
                                                                                                                     tGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAAt                          >  2:297754/1‑150 (MQ=255)
                                                                                                                       cTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTc                                                                                     <  4:434480/89‑1 (MQ=255)
                                                                                                                       cTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTc                                                                                     >  3:434480/1‑89 (MQ=255)
                                                                                                                       cTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATgg                        >  4:683314/1‑150 (MQ=255)
                                                                                                                        tttAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAggatcggat                                                                                                            <  4:157194/65‑1 (MQ=255)
                                                                                                                        tttAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAggatcggat                                                                                                            >  3:157194/1‑65 (MQ=255)
                                                                                                                         ttAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGtt                      >  2:508054/1‑150 (MQ=255)
                                                                                                                          tAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGGTCGTCAACAGGGCGCTTGACGATGGGATGGAGCTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTg                     <  1:676117/150‑1 (MQ=255)
                                                                                                                          tAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGaaa                                                 >  1:917853/1‑122 (MQ=255)
                                                                                                                          tAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGaaa                                                 <  2:917853/122‑1 (MQ=255)
                                                                                                                            aGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGaa                           >  3:473688/1‑142 (MQ=255)
                                                                                                                            aGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGaa                           <  4:473688/142‑1 (MQ=255)
                                                                                                                            aGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGcc                   >  1:851639/1‑150 (MQ=255)
                                                                                                                            aGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGcc                   >  2:642795/1‑150 (MQ=255)
                                                                                                                                 gTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCa                                                                               >  1:73449/1‑85 (MQ=255)
                                                                                                                                 gTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCa                                                                               <  2:73449/85‑1 (MQ=255)
                                                                                                                                    aGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCacgccg           >  2:308730/1‑150 (MQ=255)
                                                                                                                                     gCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCAcgacga          >  1:640495/1‑150 (MQ=255)
                                                                                                                                          gCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACGAATCGc     <  3:119432/150‑1 (MQ=255)
                                                                                                                                           cGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACGAATCGcc    <  4:454674/150‑1 (MQ=255)
                                                                                                                                            ggACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACGAATCGCct   <  2:390084/150‑1 (MQ=255)
                                                                                                                                             gACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACGAATCGCctc  >  1:404314/1‑150 (MQ=255)
                                                                                                                                             gACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACGAATCGCctc  <  4:448704/150‑1 (MQ=255)
                                                                                                                                             gACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACGAATCGCctc  <  2:348345/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
TCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCCAGAATGGTTGCCGCACGACGAATCGCCTC  >  NC_000913_3_pae_pgi/3440097‑3440387

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: