Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913986,5600AG38.4% 86.4 / 147.4 139G141G (GGT→GGC) ompFouter membrane porin 1a (Ia;b;F)
Reads supporting (aligned to +/- strand):  ref base A (39/46);  new base G (24/29);  total (63/76)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.72e-02
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 152 total aligned reads displayed.

GTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCC  >  NC_000913/986461‑986659
                                                                                                   |                                                                                                   
gTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATcacca                                                                                                     >  1:1790419/1‑100 (MQ=255)
   cAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCAccacca                                                                                                  >  2:1584849/1‑100 (MQ=255)
   cAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCAccacca                                                                                                  <  2:1743569/100‑1 (MQ=255)
     ggCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCaaa                                                                                                >  1:2249482/1‑100 (MQ=39)
      gCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCaaac                                                                                               >  1:389008/1‑99 (MQ=39)
      gCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCaaac                                                                                               <  2:1688076/100‑2 (MQ=38)
       ccATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAAtt                                                                                              >  2:920872/1‑100 (MQ=255)
       ccATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAAtt                                                                                              >  1:2263174/1‑100 (MQ=255)
       ccATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAAtt                                                                                              >  2:2562050/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)
          tCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTg                                                                                           <  2:666862/100‑1 (MQ=255)
          tCAACCAGACCAAAGAAGTTGGAGTTACGAGAGGTAGCAACGCCGCCAACGCGACCAACGAAGAAGTCAGCGCTGTATGCAGTATCACCACCAAATTCTg                                                                                           <  1:23240/100‑1 (MQ=255)
            aaCCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGgc                                                                                         <  2:78577/100‑1 (MQ=255)
             aCCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAAAGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGgca                                                                                        >  1:476450/1‑100 (MQ=255)
               cAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGgcagc                                                                                      <  1:947140/100‑1 (MQ=255)
                aGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGgcagca                                                                                     <  2:184754/100‑1 (MQ=255)
                aGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGgcagca                                                                                     <  2:2310366/100‑1 (MQ=255)
                   ccAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCatat                                                                                  >  1:1025916/1‑100 (MQ=255)
                     aaaGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCg                                                                                >  1:60616/1‑100 (MQ=255)
                     aaaGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCg                                                                                >  1:2376687/1‑100 (MQ=255)
                      aagaagTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCgg                                                                               <  1:2376512/100‑1 (MQ=255)
                      aagaagTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCgg                                                                               <  2:1167459/100‑1 (MQ=255)
                      aagaagTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCgg                                                                               <  1:871177/100‑1 (MQ=255)
                        gaagTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGgtg                                                                             <  2:375956/100‑1 (MQ=255)
                                 gTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGt                                                                    >  1:512209/1‑100 (MQ=38)
                                  ttACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTg                                                                   <  2:1663794/100‑1 (MQ=255)
                                   tACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGc                                                                  >  1:415469/1‑100 (MQ=255)
                                    aCGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGca                                                                 <  2:1465256/100‑1 (MQ=38)
                                     cGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGcat                                                                <  1:1881806/100‑1 (MQ=255)
                                      gATAGGTAGCAACGCCGCCAACACGACCAACGAAGTAGTCGTCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGcatc                                                               >  2:2532933/1‑100 (MQ=255)
                                      gATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGcatc                                                               <  2:1659208/100‑1 (MQ=255)
                                       aTAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGcatca                                                              >  1:1167522/1‑100 (MQ=38)
                                          ggTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATaa                                                           >  1:845596/1‑100 (MQ=255)
                                          ggTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATaa                                                           <  2:544030/100‑1 (MQ=255)
                                           gTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATaaa                                                          >  1:75202/1‑100 (MQ=38)
                                            tAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATaaan                                                         <  1:2113864/100‑2 (MQ=255)
                                             ggCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAAcc                                                        >  1:370123/2‑100 (MQ=38)
                                             aGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGGTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAAcc                                                        <  2:2357353/100‑1 (MQ=37)
                                             aGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAAcc                                                        <  2:603457/100‑1 (MQ=38)
                                             aGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAAcc                                                        <  2:232579/100‑1 (MQ=38)
                                             aGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAAcc                                                        <  2:682629/100‑1 (MQ=38)
                                               cAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACcac                                                      <  2:918968/100‑1 (MQ=255)
                                               cAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACcac                                                      >  1:1795486/1‑100 (MQ=255)
                                                naCGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTAGCACCACCAAATTCTGGCAGCATATCGGGGTAACCCAGTGCATCAGAAACcaca                                                     >  1:2113983/2‑100 (MQ=255)
                                                aaCGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACcaca                                                     >  2:1786649/1‑100 (MQ=255)
                                                 aCGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCacac                                                    <  2:2580217/100‑1 (MQ=255)
                                                  cgccgcCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACAcc                                                   >  1:2386289/1‑100 (MQ=38)
                                                  cgccgcCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACAcc                                                   >  1:2149892/1‑100 (MQ=38)
                                                  cgccgcCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACAcc                                                   >  1:1484546/1‑100 (MQ=38)
                                                   gccgccAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCg                                                  >  1:1702312/1‑100 (MQ=38)
                                                   gccgccAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCg                                                  >  1:449368/1‑100 (MQ=38)
                                                   gccgccAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCg                                                  >  1:1955816/1‑100 (MQ=38)
                                                    ccgccAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGt                                                 <  2:1791476/100‑1 (MQ=38)
                                                      gccAACGCGACCACCGAAGAAGTCTTCGCTGTATGCCCTGTAACCGCCAAACTCAGGCAGCATAGCGGTGTAACCCAGTGCATCATAAACCACACCGTAg                                               <  2:1604335/100‑1 (MQ=255)
                                                      gccAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAg                                               <  2:2208070/100‑1 (MQ=38)
                                                      gccAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAg                                               <  2:1487797/100‑1 (MQ=38)
                                                      gccAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAg                                               <  2:1529888/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)
                                                           aactgccAACGAAAAGGTCATCGCTGTATGCAGTATCACCACCAAATACTGGCAGCATATCGGGGTAACCCAGTGTATCATAAACCACACCGTAGTTACg                                          <  2:2075618/95‑1 (MQ=255)
                                                             cGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGc                                        >  2:1861456/1‑100 (MQ=255)
                                                              gACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGcc                                       <  1:1686432/100‑1 (MQ=255)
                                                                ccAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGt                                     >  2:2424157/1‑100 (MQ=255)
                                                                ccAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGt                                     <  2:1320542/100‑1 (MQ=255)
                                                                ccAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGt                                     >  2:322852/1‑100 (MQ=255)
                                                                    cGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATc                                 >  2:1404601/1‑100 (MQ=255)
                                                                    cGAAGAAGGCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATc                                 <  2:1321000/100‑1 (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)
                                                                        gaagTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGaaa                             <  1:2241415/100‑1 (MQ=255)
                                                                             cATCGCTGTATGCCGTGTCACCACCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAAcc                        >  1:2207538/1‑100 (MQ=255)
                                                                              aTCGCTGTATTCAGTATCACCACCAAAGGCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACGCCGTAGTTACGGCCGTAATCGAAAGAACCa                       <  2:1417284/100‑1 (MQ=255)
                                                                                 gCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTCCGGCCGTAATCGAAAGAACCAACg                    >  1:1169936/1‑100 (MQ=255)
                                                                                 gCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACg                    >  2:877330/1‑100 (MQ=255)
                                                                                  cTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGt                   <  1:1921459/100‑1 (MQ=255)
                                                                                     tATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAg                <  1:2495391/100‑1 (MQ=255)
                                                                                      aTGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGc               >  1:998172/1‑100 (MQ=38)
                                                                                      aTGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGc               >  1:1988973/1‑100 (MQ=38)
                                                                                      aTGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGc               >  1:683724/1‑100 (MQ=38)
                                                                                      aTGCCGTGTCACCGCCAAACGCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCAAACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAgt               <  2:1058746/100‑2 (MQ=255)
                                                                                      aTGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGc               >  1:421804/1‑100 (MQ=255)
                                                                                       tGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCg              >  1:1960501/1‑100 (MQ=38)
                                                                                       tGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCg              <  1:1478391/100‑1 (MQ=255)
                                                                                           gTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTAtt          >  1:1009060/1‑100 (MQ=255)
                                                                                           gTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTAtt          <  1:930731/100‑1 (MQ=255)
                                                                                               caccgcCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAg      >  1:691436/1‑100 (MQ=255)
                                                                                               caccacCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAg      >  2:2235593/1‑100 (MQ=255)
                                                                                                accaccAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGa     >  2:585856/1‑100 (MQ=255)
                                                                                                 ccaccaAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGAc    >  1:200388/1‑100 (MQ=255)
                                                                                                  caccaAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGAcc   >  1:728836/1‑100 (MQ=255)
                                                                                                   accaAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGAccc  <  1:2161553/100‑1 (MQ=255)
                                                                                                   accaAATTCTGGCAGCATATCCGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGCCAGCGTATTTAAGaaaa  >  2:248614/1‑97 (MQ=255)
                                                                                                   |                                                                                                   
GTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCC  >  NC_000913/986461‑986659

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: ‑