Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913986,5690TA38.0% 147.2 / 165.5 150P138P (CCA→CCT) ompFouter membrane porin 1a (Ia;b;F)
Reads supporting (aligned to +/- strand):  ref base T (41/52);  new base A (29/28);  total (70/80)
Fisher's exact test for biased strand distribution p-value = 5.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 153 total aligned reads displayed.

ATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCC  >  NC_000913/986470‑986668
                                                                                                   |                                                                                                   
  cAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCTGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTgg                                                                                                   <  1:1984325/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)
    aCCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAAAGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGgca                                                                                                 >  1:476450/1‑100 (MQ=255)
     ccAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGgcag                                                                                                >  2:2047758/1‑100 (MQ=255)
     ccAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGgcag                                                                                                >  1:2182259/1‑100 (MQ=255)
      cAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGgcagc                                                                                               <  1:947140/100‑1 (MQ=255)
      cAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGgcagc                                                                                               <  1:2219079/100‑1 (MQ=255)
       aGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGgcagca                                                                                              <  2:184754/100‑1 (MQ=255)
       aGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCCCCACCAAATTCTGgccgca                                                                                              >  2:2547403/1‑100 (MQ=255)
       aGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGgcagca                                                                                              >  2:668462/1‑100 (MQ=255)
       aGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGgcagca                                                                                              <  2:2310366/100‑1 (MQ=255)
       aGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGgcagca                                                                                              <  2:1371203/100‑1 (MQ=255)
          ccAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCatat                                                                                           >  1:1025916/1‑100 (MQ=255)
           cAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATc                                                                                          >  2:1529939/1‑100 (MQ=255)
            aaaGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCg                                                                                         <  2:384969/100‑1 (MQ=255)
            aaaGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCg                                                                                         >  1:617900/1‑100 (MQ=255)
            aaaGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCg                                                                                         >  1:60616/1‑100 (MQ=255)
            aaaGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCg                                                                                         >  2:1024821/1‑100 (MQ=255)
             acgaagTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGCAGTCATCGCTGTATGCCGTGCCACCGCCAAACTCAGGCAGCATATCgg                                                                                        <  2:2583378/98‑1 (MQ=255)
             aagaagTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCgg                                                                                        <  2:1167459/100‑1 (MQ=255)
             aagaagTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCgg                                                                                        <  1:871177/100‑1 (MQ=255)
             aagaagTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCgg                                                                                        <  1:2376512/100‑1 (MQ=255)
                  gTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTaa                                                                                   >  1:1432119/1‑100 (MQ=38)
                  gTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTaa                                                                                   <  2:1644077/100‑1 (MQ=255)
                    tGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAAcc                                                                                 >  1:2074744/1‑100 (MQ=38)
                    cggAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAAcc                                                                                 >  1:1710895/2‑100 (MQ=38)
                           aCGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGca                                                                          <  2:1465256/100‑1 (MQ=38)
                             gATAGGTAGCAACGCCGCCAACACGACCAACGAAGTAGTCGTCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGcatc                                                                        >  2:2532933/1‑100 (MQ=255)
                             gATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGcatc                                                                        <  2:1659208/100‑1 (MQ=255)
                             gATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGCGcatc                                                                        >  2:1485870/1‑100 (MQ=255)
                              aTAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGcatca                                                                       >  1:1167522/1‑100 (MQ=38)
                              aTAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGcatca                                                                       <  1:2123115/100‑1 (MQ=255)
                                 ggTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATaa                                                                    >  1:1865114/1‑100 (MQ=255)
                                 ggTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATaa                                                                    <  2:544030/100‑1 (MQ=255)
                                   tAGCAACGCCGCCAACACGACCAATGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATaaan                                                                  <  1:711453/100‑2 (MQ=255)
                                   tAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATaaan                                                                  <  1:2113864/100‑2 (MQ=255)
                                    aGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAAcc                                                                 <  2:207671/100‑1 (MQ=38)
                                    aGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAAcc                                                                 <  2:603457/100‑1 (MQ=38)
                                    aGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAAcc                                                                 <  2:682629/100‑1 (MQ=38)
                                     gCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACca                                                                >  1:1944292/1‑100 (MQ=255)
                                      cAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACcac                                                               <  2:918968/100‑1 (MQ=255)
                                       aaCGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACcaca                                                              >  2:1786649/1‑100 (MQ=255)
                                         cgccgcCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACAcc                                                            >  1:2149892/1‑100 (MQ=38)
                                         cgccgcCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACAcc                                                            >  1:2386289/1‑100 (MQ=38)
                                          gccgccAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCg                                                           >  1:445487/1‑100 (MQ=38)
                                          gccgccAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCg                                                           >  1:1702312/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)
                                           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:2311402/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)
                                               cAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGtt                                                      <  2:2226196/100‑1 (MQ=38)
                                                aaCACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTa                                                     <  2:2502603/100‑1 (MQ=255)
                                                    cGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGc                                                 >  2:1861456/1‑100 (MQ=255)
                                                     gACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGcc                                                <  1:1686432/100‑1 (MQ=255)
                                                     gACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGcc                                                >  1:64962/1‑100 (MQ=255)
                                                       ccAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGt                                              >  2:2424157/1‑100 (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)
                                                              agaagTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGaa                                       <  2:1598767/100‑1 (MQ=255)
                                                              agaagTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGaa                                       >  2:2034569/1‑100 (MQ=255)
                                                               taaatCATCGCTGTAGGCAGTCTCCCCTCCAAATTCTAGCAGCAGGTCGGCGTAACCCAGTGCATCATAATCCCCACCGTAGTTACGGCCGTAATCGaaa                                      <  2:2210630/96‑1 (MQ=255)
                                                               gaagTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGaaa                                      <  1:2241415/100‑1 (MQ=255)
                                                               gaagTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGaaa                                      <  2:1979841/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)
                                                                             aTGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGc                        >  1:998172/1‑100 (MQ=38)
                                                                             aTGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGc                        >  1:683724/1‑100 (MQ=38)
                                                                             aTGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGc                        >  1:421804/1‑100 (MQ=255)
                                                                             aTGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGc                        >  1:389236/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)
                                                                                      caccgcCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAg               >  1:691436/1‑100 (MQ=255)
                                                                                      caccgcCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAg               >  1:2273872/1‑100 (MQ=255)
                                                                                        ccaccaAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGAc             >  1:200388/1‑100 (MQ=255)
                                                                                        ccaccaAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGAc             >  1:904505/1‑100 (MQ=255)
                                                                                         ngccaAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGAcc            >  1:302758/3‑100 (MQ=255)
                                                                                         cgccaAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGAcc            >  1:78284/3‑100 (MQ=255)
                                                                                         caccaAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGAcc            >  1:728836/1‑100 (MQ=255)
                                                                                          accaAATTCTGGCAGCATATCCGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGCCAGCGTATTTAAGaaaa           >  2:248614/1‑97 (MQ=255)
                                                                                           ccaAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCcg          <  2:2147981/100‑1 (MQ=255)
                                                                                            caAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCcgc         <  1:1708767/100‑1 (MQ=255)
                                                                                             aaaTTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCcgcg        <  2:936475/100‑1 (MQ=255)
                                                                                                 tCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATg    <  2:2059787/100‑1 (MQ=255)
                                                                                                 tCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATg    <  1:2200595/100‑1 (MQ=255)
                                                                                                  cTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGc   <  1:1806954/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
ATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCC  >  NC_000913/986470‑986668

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