Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913986,5540AG37.7% 124.8 / 144.0 147D143D (GAT→GAC) ompFouter membrane porin 1a (Ia;b;F)
Reads supporting (aligned to +/- strand):  ref base A (41/50);  new base G (24/31);  total (65/82)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 148 total aligned reads displayed.

AGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTT  >  NC_000913/986455‑986652
                                                                                                   |                                                                                                  
aGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTCACAATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTa                                                                                                    <  2:2476043/100‑1 (MQ=255)
  cGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATc                                                                                                  <  2:29822/100‑1 (MQ=255)
   gAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATTCATTATca                                                                                                 <  2:25759/100‑1 (MQ=255)
      gTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATcacca                                                                                              >  1:1790419/1‑100 (MQ=255)
       ttCAGGCCATCAACCAGACCAAAGAAGTGGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTcaccgc                                                                                             <  2:1862364/100‑3 (MQ=38)
        tCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCaccacc                                                                                            <  1:1611392/100‑1 (MQ=255)
         cAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCAccacca                                                                                           <  2:1743569/100‑1 (MQ=255)
            gCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCaaac                                                                                        <  2:800424/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)
              cATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTc                                                                                      <  2:2542813/100‑1 (MQ=255)
               aTCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCt                                                                                     <  2:1546912/100‑1 (MQ=255)
                tCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTg                                                                                    <  1:628115/100‑1 (MQ=255)
                tCAACCAGACCAAAGAAGTTGGAGTTACGAGAGGTAGCAACGCCGCCAACGCGACCAACGAAGAAGTCAGCGCTGTATGCAGTATCACCACCAAATTCTg                                                                                    <  1:23240/100‑1 (MQ=255)
                 cAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTgg                                                                                   >  1:295386/1‑100 (MQ=255)
                  aaCCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGgc                                                                                  <  1:2385573/100‑1 (MQ=255)
                   aCCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGgca                                                                                 <  2:58797/100‑1 (MQ=255)
                    ccAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGgcag                                                                                >  1:2182259/1‑100 (MQ=255)
                     cAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGgcagc                                                                               <  1:2219079/100‑1 (MQ=255)
                      aGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGgcagca                                                                              <  2:2310366/100‑1 (MQ=255)
                      aGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGgcagca                                                                              >  2:668462/1‑100 (MQ=255)
                      aGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGgcagca                                                                              <  2:1371203/100‑1 (MQ=255)
                         ccAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCatat                                                                           >  1:1025916/1‑100 (MQ=255)
                         ccAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCatat                                                                           >  1:2254619/1‑100 (MQ=255)
                          cAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATc                                                                          >  2:1529939/1‑100 (MQ=255)
                           aaaGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCg                                                                         <  2:384969/100‑1 (MQ=255)
                           aaaGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCg                                                                         >  2:1024821/1‑100 (MQ=255)
                           aaaGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCg                                                                         >  1:617900/1‑100 (MQ=255)
                            acgaagTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGCAGTCATCGCTGTATGCCGTGCCACCGCCAAACTCAGGCAGCATATCgg                                                                        <  2:2583378/98‑1 (MQ=255)
                            aagaagTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCgg                                                                        <  1:871177/100‑1 (MQ=255)
                            aagaagTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCgg                                                                        <  2:1167459/100‑1 (MQ=255)
                            aagaagTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCgg                                                                        <  1:2376512/100‑1 (MQ=255)
                              gaagTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGgtg                                                                      <  2:375956/100‑1 (MQ=255)
                                agTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTa                                                                    <  2:919445/100‑1 (MQ=38)
                                 gTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTaa                                                                   >  1:1432119/1‑100 (MQ=38)
                                 gTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTaa                                                                   <  2:1644077/100‑1 (MQ=255)
                                   tGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAAcc                                                                 >  1:2074744/1‑100 (MQ=38)
                                       gTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGt                                                             >  1:512209/1‑100 (MQ=38)
                                        ttACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTg                                                            <  2:1663794/100‑1 (MQ=255)
                                         tACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGc                                                           >  1:415469/1‑100 (MQ=255)
                                           cGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGcat                                                         <  1:1881806/100‑1 (MQ=255)
                                            gATAGGTAGCAACGCCGCCAACACGACCAACGAAGTAGTCGTCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGcatc                                                        >  2:2532933/1‑100 (MQ=255)
                                            gATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGcatc                                                        <  2:1659208/100‑1 (MQ=255)
                                            gATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGCGcatc                                                        >  2:1485870/1‑100 (MQ=255)
                                             aTAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGcatca                                                       <  1:2123115/100‑1 (MQ=255)
                                                 gTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATaaa                                                   >  1:75202/1‑100 (MQ=38)
                                                   ggCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAAcc                                                 >  1:370123/2‑100 (MQ=38)
                                                   aGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGGTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAAcc                                                 <  2:2357353/100‑1 (MQ=37)
                                                   aGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAAcc                                                 <  2:207671/100‑1 (MQ=38)
                                                   aGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAAcc                                                 <  2:603457/100‑1 (MQ=38)
                                                   aGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAAcc                                                 <  2:232579/100‑1 (MQ=38)
                                                   aGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAAcc                                                 <  2:682629/100‑1 (MQ=38)
                                                     cAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACcac                                               >  1:1795486/1‑100 (MQ=255)
                                                      naCGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTAGCACCACCAAATTCTGGCAGCATATCGGGGTAACCCAGTGCATCAGAAACcaca                                              >  1:2113983/2‑100 (MQ=255)
                                                        cgccgcCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACAcc                                            >  1:1484546/1‑100 (MQ=38)
                                                        cgccgcCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACAcc                                            >  1:2386289/1‑100 (MQ=38)
                                                         gccgccAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCg                                           >  1:445487/1‑100 (MQ=38)
                                                         gccgccAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCg                                           >  1:449368/1‑100 (MQ=38)
                                                         gccgccAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCg                                           >  1:1702312/1‑100 (MQ=38)
                                                         gccgccAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCg                                           >  1:1955816/1‑100 (MQ=38)
                                                          ccgccAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGt                                          <  2:1791476/100‑1 (MQ=38)
                                                          ccgccAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGt                                          >  1:795293/1‑100 (MQ=38)
                                                          ccgccAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGt                                          <  1:2369696/100‑1 (MQ=255)
                                                           cgccAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTa                                         <  2:1412181/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:1487797/100‑1 (MQ=38)
                                                              cAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGtt                                      <  2:2226196/100‑1 (MQ=38)
                                                               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:64962/1‑100 (MQ=255)
                                                                    gACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGcc                                <  1:1686432/100‑1 (MQ=255)
                                                                      ccAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGt                              >  2:2424157/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)
                                                                             agaagTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGaa                       <  2:1598767/100‑1 (MQ=255)
                                                                              taaatCATCGCTGTAGGCAGTCTCCCCTCCAAATTCTAGCAGCAGGTCGGCGTAACCCAGTGCATCATAATCCCCACCGTAGTTACGGCCGTAATCGaaa                      <  2:2210630/96‑1 (MQ=255)
                                                                              gaagTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGaaa                      <  1:2241415/100‑1 (MQ=255)
                                                                              gaagTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGaaa                      <  2:1979841/100‑1 (MQ=255)
                                                                                   cATCGCTGTATGCCGTGTCACCACCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAAcc                 >  1:2207538/1‑100 (MQ=255)
                                                                                    aTCGCTGTATTCAGTATCACCACCAAAGGCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACGCCGTAGTTACGGCCGTAATCGAAAGAACCa                <  2:1417284/100‑1 (MQ=255)
                                                                                      cGCTGTATGCAGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAAc              >  1:510220/1‑100 (MQ=39)
                                                                                       gCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTCCGGCCGTAATCGAAAGAACCAACg             >  1:1169936/1‑100 (MQ=255)
                                                                                       gCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACg             >  2:877330/1‑100 (MQ=255)
                                                                                        cTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGt            <  1:1921459/100‑1 (MQ=255)
                                                                                          gTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCa          >  2:1847328/1‑100 (MQ=255)
                                                                                            aTGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGc        >  1:1988973/1‑100 (MQ=38)
                                                                                            aTGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGc        >  1:998172/1‑100 (MQ=38)
                                                                                            aTGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGc        >  1:683724/1‑100 (MQ=38)
                                                                                            aTGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGc        >  1:389236/1‑100 (MQ=255)
                                                                                             tGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCg       >  1:1960501/1‑100 (MQ=38)
                                                                                                 gTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTAtt   <  1:930731/100‑1 (MQ=255)
                                                                                                 gTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTAtt   >  1:1009060/1‑100 (MQ=255)
                                                                                                  tATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTAttt  <  1:1041095/100‑1 (MQ=255)
                                                                                                   |                                                                                                  
AGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTT  >  NC_000913/986455‑986652

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