Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913986,5510AC37.6% 115.0 / 173.5 149T144T (ACT→ACG) ompFouter membrane porin 1a (Ia;b;F)
Reads supporting (aligned to +/- strand):  ref base A (41/52);  new base C (23/33);  total (64/85)
Fisher's exact test for biased strand distribution p-value = 7.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.21e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 149 total aligned reads displayed.

ACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGT  >  NC_000913/986453‑986648
                                                                                                  |                                                                                                 
aCAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAg                                                                                                  <  2:414055/100‑1 (MQ=255)
   gCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTAt                                                                                               >  1:1663561/1‑100 (MQ=255)
    cGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATc                                                                                              <  2:29822/100‑1 (MQ=255)
     gAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATTCATTATca                                                                                             <  2:25759/100‑1 (MQ=255)
        gTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTcaccg                                                                                          <  2:1663663/100‑2 (MQ=255)
        gTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATcacca                                                                                          >  1:1790419/1‑100 (MQ=255)
         ttCAGGCCATCAACCAGACCAAAGAAGTGGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTcaccgc                                                                                         <  2:1862364/100‑3 (MQ=38)
           cAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCAccacca                                                                                       <  2:1743569/100‑1 (MQ=255)
           cAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCAccacca                                                                                       >  2:1584849/1‑100 (MQ=255)
             ggCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCaaa                                                                                     >  1:2249482/1‑100 (MQ=39)
              gCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCaaac                                                                                    >  1:389008/1‑99 (MQ=39)
              gCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCaaac                                                                                    <  2:1688076/100‑2 (MQ=38)
              gCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCaaac                                                                                    <  2:800424/100‑2 (MQ=38)
               ccATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAAtt                                                                                   >  2:920872/1‑100 (MQ=255)
               ccATCAACCAGACCAAAGAAGTCGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAAtt                                                                                   >  1:2211088/1‑100 (MQ=255)
                cATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTc                                                                                  <  2:2542813/100‑1 (MQ=255)
                 aTCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCt                                                                                 <  2:1546912/100‑1 (MQ=255)
                  tCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTg                                                                                <  2:1640860/100‑1 (MQ=255)
                  tCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTg                                                                                <  1:628115/100‑1 (MQ=255)
                  tCAACCAGACCAAAGAAGTTGGAGTTACGAGAGGTAGCAACGCCGCCAACGCGACCAACGAAGAAGTCAGCGCTGTATGCAGTATCACCACCAAATTCTg                                                                                <  1:23240/100‑1 (MQ=255)
                   cAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCTGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTgg                                                                               <  1:1984325/100‑1 (MQ=255)
                     aCCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGgca                                                                             <  2:58797/100‑1 (MQ=255)
                     aCCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAAAGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGgca                                                                             >  1:476450/1‑100 (MQ=255)
                      ccAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGgcag                                                                            >  1:2182259/1‑100 (MQ=255)
                       cAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGgcagc                                                                           <  1:2219079/100‑1 (MQ=255)
                       cAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGgcagc                                                                           <  1:947140/100‑1 (MQ=255)
                        aGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCCCCACCAAATTCTGgccgca                                                                          >  2:2547403/1‑100 (MQ=255)
                        aGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGgcagca                                                                          >  2:668462/1‑100 (MQ=255)
                        aGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGgcagca                                                                          <  2:1371203/100‑1 (MQ=255)
                        aGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGgcagca                                                                          <  2:2310366/100‑1 (MQ=255)
                           ccAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCatat                                                                       >  1:1025916/1‑100 (MQ=255)
                           ccAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCatat                                                                       >  1:2254619/1‑100 (MQ=255)
                             aaaGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCg                                                                     <  2:384969/100‑1 (MQ=255)
                             aaaGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCg                                                                     >  1:2376687/1‑100 (MQ=255)
                             aaaGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCg                                                                     >  1:60616/1‑100 (MQ=255)
                             aaaGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCg                                                                     >  2:1024821/1‑100 (MQ=255)
                              aagaagTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCgg                                                                    <  1:2376512/100‑1 (MQ=255)
                                  agTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTa                                                                <  2:919445/100‑1 (MQ=38)
                                     cggAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAAcc                                                             >  1:1710895/2‑100 (MQ=38)
                                         gTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGt                                                         >  1:512209/1‑100 (MQ=38)
                                           tACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGc                                                       >  1:415469/1‑100 (MQ=255)
                                            aCGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGca                                                      <  2:1465256/100‑1 (MQ=38)
                                             cGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGcat                                                     <  1:1881806/100‑1 (MQ=255)
                                              gATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGcatc                                                    <  2:1659208/100‑1 (MQ=255)
                                               aTAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGcatca                                                   >  1:1167522/1‑100 (MQ=38)
                                                  ggTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATaa                                                >  1:845596/1‑100 (MQ=255)
                                                  ggTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATaa                                                >  1:1865114/1‑100 (MQ=255)
                                                   gTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATaaa                                               <  2:2509698/100‑1 (MQ=38)
                                                    tAGCAACGCCGCCAACACGACCAATGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATaaan                                              <  1:711453/100‑2 (MQ=255)
                                                    tAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATaaan                                              <  1:2113864/100‑2 (MQ=255)
                                                     ggCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAAcc                                             >  1:370123/2‑100 (MQ=38)
                                                     aGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAAcc                                             <  2:207671/100‑1 (MQ=38)
                                                     aGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAAcc                                             <  2:682629/100‑1 (MQ=38)
                                                     aGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAAcc                                             <  2:603457/100‑1 (MQ=38)
                                                      gCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACca                                            >  1:1944292/1‑100 (MQ=255)
                                                        naCGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTAGCACCACCAAATTCTGGCAGCATATCGGGGTAACCCAGTGCATCAGAAACcaca                                          >  1:2113983/2‑100 (MQ=255)
                                                         aCGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCacac                                         <  2:2580217/100‑1 (MQ=255)
                                                          cgccgcCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACAcc                                        >  1:1484546/1‑100 (MQ=38)
                                                          cgccgcCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACAcc                                        >  1:2149892/1‑100 (MQ=38)
                                                          cgccgcCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACAcc                                        >  1:2386289/1‑100 (MQ=38)
                                                           gccgccAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCg                                       >  1:1955816/1‑100 (MQ=38)
                                                           gccgccAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCg                                       >  1:449368/1‑100 (MQ=38)
                                                           gccgccAACACGACCAACGAAGAAGTCATCGCTGTATACCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCg                                       <  2:2355914/100‑1 (MQ=37)
                                                            ccgccAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGt                                      >  1:795293/1‑100 (MQ=38)
                                                            ccgccAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGt                                      <  2:1791476/100‑1 (MQ=38)
                                                             cgccAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTa                                     <  2:1687743/100‑1 (MQ=38)
                                                             cgccAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTa                                     <  2:1261235/100‑1 (MQ=255)
                                                              gccAACGCGACCACCGAAGAAGTCTTCGCTGTATGCCCTGTAACCGCCAAACTCAGGCAGCATAGCGGTGTAACCCAGTGCATCATAAACCACACCGTAg                                    <  2:1604335/100‑1 (MQ=255)
                                                              gccAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAg                                    <  2:1529888/100‑1 (MQ=38)
                                                              gccAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAg                                    <  2:2311402/100‑1 (MQ=38)
                                                              gccAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAg                                    <  2:2208070/100‑1 (MQ=38)
                                                              gccAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAg                                    <  2:1487797/100‑1 (MQ=38)
                                                              gccAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAg                                    <  1:2311597/100‑1 (MQ=255)
                                                              gccAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAg                                    <  1:647588/100‑1 (MQ=255)
                                                                 aaCACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTa                                 <  2:2502603/100‑1 (MQ=255)
                                                                  acacGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTAc                                <  2:1747999/100‑1 (MQ=38)
                                                                   aactgccAACGAAAAGGTCATCGCTGTATGCAGTATCACCACCAAATACTGGCAGCATATCGGGGTAACCCAGTGTATCATAAACCACACCGTAGTTACg                               <  2:2075618/95‑1 (MQ=255)
                                                                     cGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGc                             >  2:1861456/1‑100 (MQ=255)
                                                                      gACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGcc                            <  1:1686432/100‑1 (MQ=255)
                                                                        ccAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGt                          >  2:322852/1‑100 (MQ=255)
                                                                        ccAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGt                          >  2:2424157/1‑100 (MQ=255)
                                                                        ccAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGt                          <  2:1320542/100‑1 (MQ=255)
                                                                            cGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATc                      >  2:1404601/1‑100 (MQ=255)
                                                                             gaagaaGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCg                     <  2:2341143/100‑1 (MQ=38)
                                                                               agaagTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTACCCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGag                   <  2:1374350/100‑2 (MQ=38)
                                                                               agaagTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGaa                   >  2:2034569/1‑100 (MQ=255)
                                                                                taaatCATCGCTGTAGGCAGTCTCCCCTCCAAATTCTAGCAGCAGGTCGGCGTAACCCAGTGCATCATAATCCCCACCGTAGTTACGGCCGTAATCGaaa                  <  2:2210630/96‑1 (MQ=255)
                                                                                gaagTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGaaa                  <  2:1979841/100‑1 (MQ=255)
                                                                                 aagTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTCACCCAGTGCATCATAAACCATACCGTAGTTACTGCCGTAATCGAAAg                 <  2:2428632/100‑1 (MQ=37)
                                                                                      aTCGCTGTATTCAGTATCACCACCAAAGGCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACGCCGTAGTTACGGCCGTAATCGAAAGAACCa            <  2:1417284/100‑1 (MQ=255)
                                                                                         gCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTCCGGCCGTAATCGAAAGAACCAACg         >  1:1169936/1‑100 (MQ=255)
                                                                                         gCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACg         >  2:877330/1‑100 (MQ=255)
                                                                                            gTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCa      >  2:1847328/1‑100 (MQ=255)
                                                                                             tATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAg     <  1:2495391/100‑1 (MQ=255)
                                                                                              aTGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGc    >  1:1988973/1‑100 (MQ=38)
                                                                                              aTGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGc    >  1:683724/1‑100 (MQ=38)
                                                                                              aTGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGc    >  1:998172/1‑100 (MQ=38)
                                                                                              aTGCCGTGTCACCGCCAAACGCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCAAACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAgt    <  2:1058746/100‑2 (MQ=255)
                                                                                              aTGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGc    >  1:421804/1‑100 (MQ=255)
                                                                                               tGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCg   <  1:1478391/100‑1 (MQ=255)
                                                                                                  |                                                                                                 
ACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGT  >  NC_000913/986453‑986648

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