Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913986,4130AG27.7% 463.9 / 216.8 311R190R (CGT→CGC) ompFouter membrane porin 1a (Ia;b;F)
Reads supporting (aligned to +/- strand):  ref base A (106/118);  new base G (17/69);  total (123/188)
Fisher's exact test for biased strand distribution p-value = 8.28e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.54e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 318 total aligned reads displayed.

TTCTTGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACG  >  NC_000913/986314‑986512
                                                                                                   |                                                                                                   
tgcttGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATGCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGaa                                                                                                     <  1:2571616/98‑1 (MQ=255)
  cttGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAAca                                                                                                   >  2:1507456/1‑99 (MQ=255)
    tGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGc                                                                                                 >  1:3658563/1‑100 (MQ=255)
     gCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAgcggcg                                                                                                >  1:618383/1‑100 (MQ=255)
     gCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCg                                                                                                >  1:2255436/1‑100 (MQ=255)
        ggTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGc                                                                                             >  2:1242539/1‑100 (MQ=255)
           tGGTACGGTGCGCTGCACCGTAAGCACCAACGATACCAAAGACTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAgt                                                                                          <  2:4533517/100‑1 (MQ=255)
            ggTACGGTCAGCTGCACCATAAGCACCAACGATACCAACGCCTTCGTATTCGTAGCTGATAGAACCGCCACCACCGTCGCCGTAAGAGCGGCGTGCAgtg                                                                                         <  2:3056009/100‑1 (MQ=255)
            ggTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAGCGGCGTGCAgtg                                                                                         <  2:1843782/100‑1 (MQ=255)
            ggTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAgtg                                                                                         <  2:2067887/100‑1 (MQ=255)
            gctACGGTTAGCTGCACCATAAGTGCCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCTCCAACCCCGTCGCCGTCGGAACCGCGTCCAgtg                                                                                         <  2:1165823/98‑1 (MQ=255)
             gTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAgtgt                                                                                        >  1:4277936/1‑100 (MQ=255)
                  gTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGc                                                                                   <  1:4983188/100‑1 (MQ=255)
                   tCAGCTGCACCATAAGCACCACCGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTGGCCGGTAGAACGGCGTGCAGTGTCACGCt                                                                                  <  1:4809362/100‑1 (MQ=255)
                      gCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGt                                                                               >  2:3814247/1‑100 (MQ=255)
                      gCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGt                                                                               >  1:2733895/1‑100 (MQ=255)
                        tGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAGCGGCGTGCAGTGTCACGCTCGttt                                                                             <  2:2736676/100‑1 (MQ=255)
                        tGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAGCGGCGTGCAGTGTCACGCTCGttt                                                                             <  2:1931440/100‑1 (MQ=255)
                        ggCACTATAAGCACCAACGATATCACAGCCTGCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGCGCAGTGTCACGCTCGttt                                                                             <  2:1551054/99‑1 (MQ=255)
                         gCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGtttt                                                                            >  1:1147979/1‑100 (MQ=255)
                         gCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGtttt                                                                            <  2:4220300/100‑1 (MQ=255)
                            ccATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTAc                                                                         >  1:2457787/1‑100 (MQ=255)
                             cATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAGCGGCGTGCAGTGTCACGCTCGTTTTTAcc                                                                        <  2:1639118/100‑1 (MQ=255)
                              aTAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAGCGGCGTGCAGTGTCACGCTCGTTTTTTccc                                                                       <  2:717256/100‑1 (MQ=255)
                              aTAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTAccc                                                                       >  2:2862451/1‑100 (MQ=255)
                               tAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAGCGGCGTGCAGTGTCACGCTCGTTTTTACCCa                                                                      <  2:3371754/100‑1 (MQ=255)
                                 aGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCATCACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAgg                                                                    <  2:5003749/100‑1 (MQ=255)
                                 aGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAGCGGCGTGCAGTGTCACGCTCGTTTTTACCCAgg                                                                    <  2:3757643/100‑1 (MQ=255)
                                 aGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAgg                                                                    <  2:777975/100‑1 (MQ=255)
                                  gCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAGCGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGt                                                                   >  1:2212392/1‑100 (MQ=255)
                                       aaCGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGa                                                              >  2:3233000/1‑100 (MQ=255)
                                        aCGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGaa                                                             >  2:714215/1‑100 (MQ=255)
                                         cGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAGCGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAAc                                                            <  2:3746804/100‑1 (MQ=255)
                                              ccAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGTCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGa                                                       >  2:2513101/1‑100 (MQ=255)
                                              ccAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGa                                                       >  2:3467271/1‑100 (MQ=255)
                                              ccAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGa                                                       >  2:1440936/1‑100 (MQ=255)
                                               cAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGaa                                                      >  2:362651/1‑100 (MQ=255)
                                                aaaGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAGCGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAg                                                     >  1:3138606/1‑100 (MQ=255)
                                                aaaGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAg                                                     <  1:4045206/100‑1 (MQ=255)
                                                 aaGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGt                                                    <  1:3242697/100‑1 (MQ=255)
                                                   gCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTc                                                  <  1:725331/100‑1 (MQ=255)
                                                    ccTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCa                                                 >  2:4233388/1‑100 (MQ=255)
                                                    ccTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCa                                                 >  1:2428140/1‑100 (MQ=255)
                                                     ctccgtattcgtaGCTGATAGAACCGCCAACACCGTCGCCGTTAGAGCGGCGTGCAGTGTCACGCTCGTATTTACCCAGGTACTGAACAGCGAAGTTCAg                                                >  1:1964860/4‑100 (MQ=255)
                                                     cTTTGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAGCGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAg                                                >  1:1349526/1‑100 (MQ=255)
                                                     cTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAg                                                <  1:2894798/100‑1 (MQ=255)
                                                      ttcgtattcgtaGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAgg                                               <  2:481936/100‑1 (MQ=255)
                                                       tcgtattcgtaGCTGATAGAACCGCCAACATCGTCGCCGTTAGGGCGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGc                                              <  2:816980/100‑1 (MQ=255)
                                                       tcgtattcgtaGCTGATAGAACCGCCAACACCGTCGCCGTTAGAGCGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGc                                              <  2:3551146/100‑1 (MQ=255)
                                                       tcgtattcgtaGCTGATAGAACCGCCAACACCGTCGCCGTTAGAGCGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGc                                              <  2:3303785/100‑1 (MQ=255)
                                                       tcgtattcgtaGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGc                                              <  2:433204/100‑1 (MQ=255)
                                                       tcgtattcgtaGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGc                                              <  2:3185473/100‑1 (MQ=255)
                                                       tcggattcgtaGCGGATAGAACCGCCAACACCGTCGCCGTTAGAGCGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGc                                              <  2:195399/100‑1 (MQ=255)
                                                        cgtattcgtaGCTGATAGAACCGCCAACACCGTCGCCGTTAGAGCGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGcc                                             <  2:4977835/100‑1 (MQ=255)
                                                        cgtattcgtaGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACTCAGGTACTGAACAGCGAAGTTCAGGcc                                             <  2:4116423/100‑1 (MQ=255)
                                                           attcgtaGCTGATAGAACCGCCAACACCGTCGCCGTTAGAGCGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATc                                          <  2:2884555/100‑1 (MQ=255)
                                                           attcgtaGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATc                                          >  1:1881759/1‑100 (MQ=255)
                                                             tcgtaGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCaa                                        <  2:1341757/100‑1 (MQ=255)
                                                              cgtaGCTGATAGAACCGCCAACACCGTCGCCGTTAGAGCGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAAc                                       <  2:2423543/100‑1 (MQ=255)
                                                              cgtaGCTGATAGAACCGCCAACACCGTCGCCGTTAGAGCGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAAc                                       <  2:4904181/100‑1 (MQ=255)
                                                               gtaGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAAcc                                      >  2:3907002/1‑100 (MQ=255)
                                                                 aGCTGATAGAACCGCCAACACCGTCGCCGTTAGAGCGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAg                                    <  2:2297136/100‑1 (MQ=255)
                                                                 aGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAg                                    >  1:3378700/1‑100 (MQ=255)
                                                                  gCTGATAGAACCGCCAACACCGTCGCCGTTAGAGCGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGa                                   >  1:5022856/1‑100 (MQ=255)
                                                                  gCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGa                                   >  2:192102/1‑100 (MQ=255)
                                                                   cTGATAGAACCGCCAACACCGTCGCCGTTAGAGCGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAgcc                                  <  2:3731469/100‑3 (MQ=255)
                                                                   cTGATAGAACCGCCAACACCGTCGCCGTTAGAGCGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGAc                                  <  2:310866/100‑1 (MQ=255)
                                                                   cTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAAcagatc                                  >  1:5073894/1‑95 (MQ=255)
                                                                    tGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGAcc                                 >  1:1380451/1‑100 (MQ=255)
                                                                    tGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGAcc                                 <  2:2256712/100‑1 (MQ=255)
                                                                    tGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGAcc                                 >  1:3831705/1‑100 (MQ=255)
                                                                     gATAGAACCGCCAACACCGTCGACGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCa                                >  1:4311042/1‑100 (MQ=255)
                                                                       tAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCaaa                              >  1:1251407/1‑100 (MQ=255)
                                                                        aGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAaag                             >  1:415439/1‑100 (MQ=255)
                                                                         gAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACGCAGGTACTGAAGCGCGAAGTTCAGGCTATCAACCAGACCAaaga                            >  1:4729017/1‑100 (MQ=255)
                                                                         gAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAaaga                            >  2:2573945/1‑100 (MQ=255)
                                                                         gAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAaaga                            >  1:3070798/1‑100 (MQ=255)
                                                                            ccGCCAACACCGTCGCCGTTAGAGCGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGt                         >  1:44314/1‑100 (MQ=255)
                                                                              gCCAACACCGTCGCCGTTAGAGCGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTg                       <  2:4020263/100‑1 (MQ=255)
                                                                              gCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTg                       <  2:1131003/100‑1 (MQ=255)
                                                                               ccAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTgg                      <  2:4615434/100‑1 (MQ=255)
                                                                               ccAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTgg                      <  2:4795894/100‑1 (MQ=255)
                                                                                cAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGa                     >  2:1798266/1‑100 (MQ=255)
                                                                                  acatCGTCGCCATGAGGGCGTCCTCCAGTCTCACGCTCGTTTTCACCCAGGTACTGAACAGCGAAGTTCAGGCCATTAACCAGACCAAAGAAGTTGTAGt                   <  2:4254792/100‑1 (MQ=18)
                                                                                   cacCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGtt                  <  1:956949/100‑1 (MQ=255)
                                                                                   cacCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGtt                  >  1:2119824/1‑100 (MQ=255)
                                                                                   cacCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGtt                  >  2:1733998/1‑100 (MQ=255)
                                                                                     ccGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTan                <  1:3539657/100‑2 (MQ=255)
                                                                                     ccGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTAc                <  1:3036744/100‑1 (MQ=255)
                                                                                      cGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACg               <  1:1541206/100‑1 (MQ=255)
                                                                                        tCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGAt             <  2:4951391/100‑1 (MQ=255)
                                                                                        tCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACACAGGTACTGAACAGCGAAGTTCCGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGAt             <  2:438579/100‑1 (MQ=255)
                                                                                          gCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAg           >  1:2683176/1‑100 (MQ=255)
                                                                                           ccGTTAGAGCGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAgg          <  2:2185820/100‑1 (MQ=255)
                                                                                           ccGTTAGAACGGCGTGCAGTGTCACGCTCGTGTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAgg          <  2:3983490/100‑1 (MQ=255)
                                                                                            cGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGt         <  2:885931/100‑1 (MQ=255)
                                                                                             gTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTa        >  1:3902759/1‑100 (MQ=255)
                                                                                               tAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAGAGAAGTTGGAGTTACGATAGGTAGc      <  2:4046293/100‑1 (MQ=255)
                                                                                               tAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGc      <  1:1044461/100‑1 (MQ=255)
                                                                                                aGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCa     >  1:2274883/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
TTCTTGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACG  >  NC_000913/986314‑986512

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