Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913986,4130AG27.5% 248.3 / 110.6 167R190R (CGT→CGC) ompFouter membrane porin 1a (Ia;b;F)
Reads supporting (aligned to +/- strand):  ref base A (46/75);  new base G (17/29);  total (63/104)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.20e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 168 total aligned reads displayed.

TCTTGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACG  >  NC_000913/986315‑986512
                                                                                                  |                                                                                                   
tcttGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAAc                                                                                                    >  1:1037291/1‑100 (MQ=255)
 cttGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACg                                                                                                   >  1:13902/1‑100 (MQ=255)
   tGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGc                                                                                                 <  1:2184393/100‑1 (MQ=255)
     cAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGt                                                                                               >  1:1952383/1‑100 (MQ=255)
      aGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTCTAACGGCgac                                                                                              <  2:114890/100‑3 (MQ=255)
       ggTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGc                                                                                             <  2:1338038/100‑1 (MQ=255)
        gTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTCTAACGGCgacgg                                                                                            <  1:1250118/100‑5 (MQ=255)
             tACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTc                                                                                       <  1:1563166/100‑1 (MQ=255)
              aCGGTCAGCTGCACCAGAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGCCGCCGTTAGAACGGCGTGCAGTGTcc                                                                                      <  2:273907/100‑2 (MQ=255)
               cGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCAc                                                                                     >  2:856300/1‑100 (MQ=255)
               cGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCAc                                                                                     <  1:943192/100‑1 (MQ=255)
                    aGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAGCGGCGTGCAGTGTCACGCTCg                                                                                <  2:1966842/100‑1 (MQ=255)
                     gCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACAGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGt                                                                               >  2:499415/1‑100 (MQ=255)
                       tGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGttt                                                                             <  2:137955/100‑1 (MQ=255)
                        gCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGtttt                                                                            >  2:1856869/1‑100 (MQ=255)
                         cACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAGCGGCGTGCAGTGTCACGCTCGttttt                                                                           <  2:1607753/100‑1 (MQ=255)
                           ccATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTAc                                                                         <  2:409186/100‑1 (MQ=255)
                           ccATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTAc                                                                         >  2:1223543/1‑100 (MQ=255)
                           ccATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGCTAGCACGGCGTGCAGTGCCACGCCCGtctttcc                                                                         >  2:2112518/1‑98 (MQ=255)
                            cATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTAcc                                                                        <  1:2618103/100‑1 (MQ=255)
                             aTAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTAccc                                                                       >  1:1533872/1‑100 (MQ=255)
                              tAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAGCGGCGTGCAGTGTCACGCTCGTTTTTACCCa                                                                      <  2:2154616/100‑1 (MQ=255)
                              tAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCa                                                                      >  1:1653282/1‑100 (MQ=255)
                              tAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCa                                                                      <  1:452363/100‑1 (MQ=255)
                              tAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCa                                                                      >  2:915333/1‑100 (MQ=255)
                              tAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCa                                                                      <  1:1486117/100‑1 (MQ=255)
                              tAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCa                                                                      <  1:1856322/100‑1 (MQ=255)
                                aGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAGCGGCGTGCAGTGTCACGCTCGTTTTTACCCAgg                                                                    <  2:55619/100‑1 (MQ=255)
                                aGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAgg                                                                    <  1:79908/100‑1 (MQ=255)
                                aGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAgg                                                                    <  1:2208640/100‑1 (MQ=255)
                                aGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAgg                                                                    <  1:2570366/100‑1 (MQ=255)
                                 gCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGt                                                                   <  2:1347305/100‑1 (MQ=255)
                                 gCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGt                                                                   >  1:680477/1‑100 (MQ=255)
                                 gCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGt                                                                   >  1:2148323/1‑100 (MQ=255)
                                        ngATGCCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTAGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAAc                                                            >  1:242375/2‑100 (MQ=255)
                                          aTACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAg                                                          <  2:2492783/100‑1 (MQ=255)
                                          aTACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAg                                                          <  2:1963284/100‑1 (MQ=255)
                                           tACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGc                                                         >  1:2355198/1‑100 (MQ=255)
                                               taaGCGTTCGTATCCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAg                                                     <  2:2421581/99‑1 (MQ=255)
                                               aaaGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAGCGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAg                                                     <  2:1345415/100‑1 (MQ=255)
                                               aaaGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAg                                                     <  2:1411536/100‑1 (MQ=255)
                                               aaaGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAg                                                     <  1:851569/100‑1 (MQ=255)
                                                    cTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAGCGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAg                                                >  1:2536750/1‑100 (MQ=255)
                                                    cTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAg                                                <  1:1019418/100‑1 (MQ=255)
                                                     ttcgtattcgtaGCTGATAGAACCGCCAACACCGTCGCCGTTAGAGCGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAgg                                               <  2:713135/100‑1 (MQ=255)
                                                     ttcgtattcgtaGCTGATAGAACCGCCAACACCGTCGCCGTTAGAGCGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAgg                                               <  2:479331/100‑1 (MQ=255)
                                                     ttcgtattcgtaGCTGATAGAACCGCCAACACCGTCGCCGTTAGAGCGGCGTGCAGTGTCACGCTCGTTTTGACCCAGGTACTGAACAGCGAAGTTCAgg                                               >  1:2401259/1‑100 (MQ=255)
                                                      tcgtattcgtaGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGc                                              <  2:814047/100‑1 (MQ=255)
                                                      tcgtattcgtaGCTGAGAGAACCGCCAACACCGTCGCCGTGAGAGCGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGc                                              <  2:2130723/100‑1 (MQ=255)
                                                       cgtattcgtaGCTTATAGAACCGCCAACACCGTCGACGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGGACTTAACAGCGAAGTTCAGGcc                                             <  2:325699/100‑1 (MQ=255)
                                                        gtattcgtaGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCa                                            <  1:1573200/100‑1 (MQ=255)
                                                        gtattcgtaGCTGATAGAACCGCCAACACCGGCGCCGTTGGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCa                                            >  1:714129/1‑100 (MQ=255)
                                                             cgtaGCTGATAGAACCGCCAACACCGTCGCCGTTAGAGCGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAAc                                       <  2:2014395/100‑1 (MQ=255)
                                                             cgtaGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAAc                                       <  2:316831/100‑1 (MQ=255)
                                                             cgtaGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAAc                                       >  1:555351/1‑100 (MQ=255)
                                                              ttaGCTGATAGAACCGCCAACACCGTCGCCGTTAGAGCGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAAcc                                      <  2:727931/99‑1 (MQ=255)
                                                              gtaGCTGATAGAACCGCCACCACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAAcc                                      <  2:2318296/100‑1 (MQ=255)
                                                              ctaGCTGATAGAACCGCCAACACCGTCGCCGTTAGAGCGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAAcc                                      <  2:105847/99‑1 (MQ=255)
                                                               taGCTGATAGAACCGCCAACACCGTCGCCGTTAGAGCGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCa                                     >  1:741889/1‑100 (MQ=255)
                                                                tgCTGATAGAACCGCCAACACCGTCGCCGTTAGAGCGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAg                                    >  1:2155704/2‑100 (MQ=255)
                                                                ggCTGATAGAACCGCCAACACCGTCGCCGTTAGAGCGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAg                                    >  1:236611/2‑100 (MQ=255)
                                                                aGCTGATAGAACCGCCAACACCGTCGCCGTTAGAGCGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTGAGGCCATCAACCAg                                    >  1:1771821/1‑100 (MQ=255)
                                                                aGCTGATAGAACCGCCAACACCGTCGCCGTTAGAGCGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAg                                    <  2:570934/100‑1 (MQ=255)
                                                                aGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAg                                    >  1:422582/1‑100 (MQ=255)
                                                                aGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAg                                    >  1:2627924/1‑100 (MQ=255)
                                                                 gCTGATAGAACCGCCAACACCGTCGCCGTTAGAGCGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTCCTGAACAGCGAAGTTCAGGCCATCAACCAGa                                   >  1:2232348/1‑100 (MQ=255)
                                                                  cTGATAGAACCGCCTACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCACGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGAc                                  >  2:956875/1‑100 (MQ=255)
                                                                  cTGATAGAACCGCCAACACCGTCGCCGTTAGAGCGGCGTGCAGTGTCACGCTCGTTTTTACCCAGCTGCTGAGCAGCGAAGTGCACGCCATCAAcctacg                                  >  1:2063598/1‑96 (MQ=255)
                                                                   tGAGAGAACCGCCAACACCGTCGCCGTTAGAGCGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGAcc                                 >  1:2029942/1‑100 (MQ=255)
                                                                    gATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCa                                >  2:634531/1‑100 (MQ=255)
                                                                     aTAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCaa                               <  1:1127113/100‑1 (MQ=255)
                                                                        gAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAaaga                            >  1:1197667/1‑100 (MQ=255)
                                                                        gAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAaaga                            >  1:1390778/1‑100 (MQ=255)
                                                                         aaCCGCCAACACCGTCGCCGTTAGAGCGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAaagaa                           <  2:2061940/100‑1 (MQ=255)
                                                                         aaCCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAaagaa                           <  1:578492/100‑1 (MQ=255)
                                                                         aaCCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAaagaa                           <  1:2320910/100‑1 (MQ=255)
                                                                          cccGCCAACACCGTCGCCGTTAGAGCGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAaagaag                          >  1:2272055/2‑100 (MQ=255)
                                                                          aCCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAaagaag                          >  2:2313289/1‑100 (MQ=255)
                                                                           ccGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCTAAGAAGt                         <  2:1229660/100‑1 (MQ=255)
                                                                           ccGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGt                         >  2:695623/1‑100 (MQ=255)
                                                                              ccAACACCGTCGCCGTTAGAGCGGCGTGCAGTGTCACGCTCGTTTTTCACCAGGTCGTGGAAAGCGAAGTTCCGGCAATCGGTCAGGCCAACGACGTTgg                      <  2:1296282/100‑1 (MQ=255)
                                                                              ccAACACCGTCGCCGTTAGAGCGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTgg                      <  2:2210933/100‑1 (MQ=255)
                                                                              ccAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTgg                      <  1:1072308/100‑1 (MQ=255)
                                                                              ccAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTgg                      <  1:329372/100‑1 (MQ=255)
                                                                                aaCACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAg                    <  2:787856/100‑1 (MQ=255)
                                                                                aaCACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAg                    <  1:1302045/100‑1 (MQ=255)
                                                                                     cggcGCCGTTAGAGCGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTATTGAACAGCGAAGTTCAGGCCATCAACCAGACCAGAGAAGTTGGAGTTACg               <  2:2552446/97‑1 (MQ=255)
                                                                                     cGTCGCCGTTAGAGCGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACg               <  2:224188/100‑1 (MQ=255)
                                                                                     cGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACg               <  2:1023564/100‑1 (MQ=255)
                                                                                     cGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACg               <  1:2337592/100‑1 (MQ=255)
                                                                                        cGCCGTTAGAGCGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATa            <  2:2504920/100‑1 (MQ=255)
                                                                                        cGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATa            <  1:1362888/100‑1 (MQ=255)
                                                                                        cGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATa            <  2:386648/100‑1 (MQ=255)
                                                                                        cGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATa            <  2:386638/100‑1 (MQ=255)
                                                                                          ccGTTAGAGCGGCGTGCAGTGTCACGCTCGGTTTTACCCAGGTTCTGAACAGCGAAGTTCAGGCCATCAACCAGTCCAAAGAAGTTGGAGTTACGATAgg          <  2:1071635/100‑1 (MQ=255)
                                                                                           cGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGt         <  2:2073126/100‑1 (MQ=255)
                                                                                             ttAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAg       <  2:1967424/100‑1 (MQ=255)
                                                                                              tAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGc      >  2:1685641/1‑100 (MQ=255)
                                                                                                gAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCaa    >  1:1077310/1‑100 (MQ=255)
                                                                                                  aCGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTGACGATAGGTAGCAAcg  >  1:59709/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
TCTTGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACG  >  NC_000913/986315‑986512

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