Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913986,5690TA28.1% 569.6 / 235.6 307P138P (CCA→CCT) ompFouter membrane porin 1a (Ia;b;F)
Reads supporting (aligned to +/- strand):  ref base T (118/102);  new base A (49/37);  total (168/139)
Fisher's exact test for biased strand distribution p-value = 6.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 313 total aligned reads displayed.

ATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCC  >  NC_000913/986470‑986668
                                                                                                   |                                                                                                   
aTCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCt                                                                                                     >  2:844602/1‑100 (MQ=255)
aTCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCt                                                                                                     >  1:4911958/1‑100 (MQ=255)
 tCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTg                                                                                                    <  1:3906900/100‑1 (MQ=255)
   aaCCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGgc                                                                                                  >  1:1670103/1‑100 (MQ=255)
   aaCCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGgc                                                                                                  >  1:2981211/1‑100 (MQ=255)
   aaCCAGACAAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGgc                                                                                                  <  2:3830481/100‑1 (MQ=255)
      cAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGgcagc                                                                                               <  2:2894259/100‑1 (MQ=255)
      cAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGgcagc                                                                                               <  1:2080502/100‑1 (MQ=255)
           cAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATc                                                                                          >  1:3350114/1‑100 (MQ=255)
           cAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATc                                                                                          <  1:803950/100‑1 (MQ=255)
           cAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATc                                                                                          >  1:2918225/1‑100 (MQ=255)
            aaaGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCg                                                                                         >  2:4927616/1‑100 (MQ=255)
             aagaagTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGGCCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTATGGCAGCATATCgg                                                                                        <  2:954520/100‑1 (MQ=255)
             aagaagTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCgg                                                                                        >  1:2182441/1‑100 (MQ=255)
               gaagTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGgtg                                                                                      >  2:4233343/1‑100 (MQ=255)
                  gTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTaa                                                                                   >  1:4974085/1‑100 (MQ=38)
                  gTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTaa                                                                                   >  1:3172081/1‑100 (MQ=255)
                  gTTGGAGTTACGATAGGTAGCAACACCGCGAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAAATCAGGCAGCATATCGGTGTaa                                                                                   >  1:4809796/1‑100 (MQ=37)
                     ggAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAAccc                                                                                >  1:470240/1‑100 (MQ=255)
                       aGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAg                                                                              <  2:3618632/100‑1 (MQ=255)
                        gTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGt                                                                             >  1:4327857/1‑100 (MQ=255)
                          tACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGc                                                                           <  2:600434/100‑1 (MQ=255)
                          tACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAATGc                                                                           >  1:2329687/1‑100 (MQ=255)
                            cGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGcat                                                                         >  2:1517642/1‑100 (MQ=255)
                             gATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGcatc                                                                        <  2:4326143/100‑1 (MQ=255)
                              atgggTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGcatca                                                                       >  1:4049859/4‑100 (MQ=37)
                                 ggTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATaa                                                                    <  2:1889541/100‑1 (MQ=38)
                                 ggTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATaa                                                                    >  1:3138319/1‑100 (MQ=255)
                                 ggTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATaa                                                                    <  2:2190880/100‑1 (MQ=255)
                                 ggTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATaa                                                                    >  1:4437577/1‑100 (MQ=255)
                                  gTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATaaa                                                                   >  1:3314785/1‑100 (MQ=255)
                                   tAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAAc                                                                  <  2:3439114/100‑1 (MQ=38)
                                   tAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAAc                                                                  <  2:1172098/100‑1 (MQ=255)
                                    aGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAAcc                                                                 <  2:2961752/100‑1 (MQ=255)
                                    aGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAAcc                                                                 >  2:4524215/1‑100 (MQ=255)
                                     gCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACca                                                                >  1:5127651/1‑100 (MQ=255)
                                     gCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACca                                                                >  1:2321083/1‑100 (MQ=255)
                                        ccgccgcCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCacac                                                             >  1:3627571/2‑100 (MQ=38)
                                         cgccgcCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACAcc                                                            >  1:757181/1‑100 (MQ=38)
                                         cgccgcCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACAcc                                                            >  1:3481622/1‑100 (MQ=255)
                                          gccgccAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCg                                                           >  1:1299892/1‑100 (MQ=38)
                                          gccgccAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCg                                                           >  1:2422386/1‑100 (MQ=38)
                                            cgccAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTa                                                         >  1:3248977/1‑100 (MQ=38)
                                            cgccAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTa                                                         >  2:890508/1‑100 (MQ=255)
                                             gccAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAg                                                        <  2:820721/100‑1 (MQ=38)
                                             gccAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAg                                                        <  2:1279324/100‑1 (MQ=38)
                                             gccAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAg                                                        <  2:1802908/100‑1 (MQ=255)
                                             gccAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAg                                                        >  1:655867/1‑100 (MQ=255)
                                              ccAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGt                                                       >  2:1702929/1‑100 (MQ=255)
                                                  caggACGAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTTCATCATAAACCACACCGTAGTTACg                                                   <  2:1115343/97‑1 (MQ=255)
                                                  cacGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACg                                                   <  2:3947217/100‑1 (MQ=38)
                                                    cGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCTCCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGc                                                 <  1:3375668/100‑1 (MQ=255)
                                                    cGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGc                                                 <  1:1350777/100‑1 (MQ=255)
                                                    cGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGc                                                 >  1:4163298/1‑100 (MQ=255)
                                                    cGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGc                                                 <  2:712102/100‑1 (MQ=255)
                                                    cGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACAACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGc                                                 <  2:4613366/100‑1 (MQ=255)
                                                      aCCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCATCATATCGGTGTAAGCCAGTGCATCATAAACCACACCGTAGTTACGGCCg                                               >  2:5003971/1‑100 (MQ=255)
                                                      aCCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCg                                               >  1:1928114/1‑100 (MQ=255)
                                                         aaCGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTaa                                            >  2:1599753/1‑100 (MQ=255)
                                                         aaCGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTaa                                            >  1:1246088/1‑100 (MQ=255)
                                                          aCGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAAt                                           <  2:925711/100‑1 (MQ=255)
                                                             aagaagTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGa                                        <  1:2405283/100‑1 (MQ=255)
                                                              agaagTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGaa                                       <  2:2394013/100‑1 (MQ=38)
                                                              agaagTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGaa                                       <  1:2603860/100‑1 (MQ=255)
                                                              agaacTCATCGCTGTATGCAGTATCACCACCAAATTGTGGCAGCATATCAGTGTAACCCAGTGCATCATGAACCACACCGTAGTTACGGCCGTGATCGaa                                       >  2:4773298/1‑100 (MQ=255)
                                                                aagTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAg                                     <  2:4365734/100‑1 (MQ=255)
                                                                   tCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTGACGGCCGTAATCGAAAGAAc                                  <  2:2302065/100‑1 (MQ=255)
                                                                    cATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAAcc                                 <  2:3081520/100‑1 (MQ=255)
                                                                    cATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAAcc                                 <  1:3065433/100‑1 (MQ=255)
                                                                    cATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAAcc                                 <  1:1042725/100‑1 (MQ=255)
                                                                       cGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATGGCGGTGTAACCCAGTGCATCATAAACCACCCCGCAGTTACGGCCGTAATCGAAAGAACCAAc                              <  2:4896134/100‑1 (MQ=255)
                                                                        gCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACg                             <  2:331774/100‑1 (MQ=255)
                                                                            gttGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAATCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAg                         >  1:4388345/3‑100 (MQ=255)
                                                                             gtGCAGTATCACCACCAAATTCTGGCAGCATATCGTTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGc                        >  1:1900629/2‑100 (MQ=255)
                                                                              tGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCg                       >  2:2108471/1‑100 (MQ=255)
                                                                               gCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGt                      >  1:1768811/1‑100 (MQ=255)
                                                                                cAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTa                     >  2:3160147/1‑100 (MQ=255)
                                                                                 cgtgtCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTAt                    >  1:253367/5‑100 (MQ=38)
                                                                                 cgtgtCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTAt                    >  1:3850325/5‑100 (MQ=38)
                                                                                     tCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCATACCGTAGTTACGGCCGTAATCCAAAGAACCAACGTCAGCGTATTTaa                <  1:4612970/100‑1 (MQ=255)
                                                                                      caccacCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGATACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAg               >  1:4550249/1‑100 (MQ=255)
                                                                                        ccgccaAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGAc             >  1:5066618/4‑100 (MQ=255)
                                                                                        ccaccaAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGAc             <  1:2030837/100‑1 (MQ=255)
                                                                                         cgccaAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGAcc            >  1:1118294/3‑100 (MQ=255)
                                                                                         cgccaAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGAcc            <  2:2622687/98‑1 (MQ=255)
                                                                                         cgccaAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGAcc            >  1:987254/3‑100 (MQ=255)
                                                                                         cgccaAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGAcc            >  1:3420375/3‑100 (MQ=255)
                                                                                         cgccaAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGAcc            >  1:411682/3‑100 (MQ=255)
                                                                                          accaAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGAccc           >  1:4675928/1‑100 (MQ=255)
                                                                                          accaAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGCATTTAAGAccc           >  1:375011/1‑100 (MQ=255)
                                                                                            caAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCcgc         <  2:2410833/100‑1 (MQ=255)
                                                                                            caAATTCTGGCAGCATATAGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGAAATCGAAAGAACCAACGTCAGCGTATTTAAGACCcgc         <  1:937802/100‑1 (MQ=255)
                                                                                             aaaTTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAAAATATTTAAGAcaagcg        <  2:1380654/100‑6 (MQ=39)
                                                                                               aTTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGaa      <  2:3026423/100‑1 (MQ=255)
                                                                                                ttCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAAt     <  2:4114257/100‑1 (MQ=255)
                                                                                                  cTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGc   >  1:2005699/1‑100 (MQ=255)
                                                                                                  cTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGc   >  1:3328621/1‑100 (MQ=255)
                                                                                                  cTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGCCAGCGTATTCAAGACCCGCGCATGc   >  2:1017199/1‑100 (MQ=255)
                                                                                                   tGGCAGCATGTCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGGCAGCGTATTTAAGACCCGCGAATGcc  <  1:1539449/100‑1 (MQ=255)
                                                                                                   tGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGcc  <  2:4266970/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: ‑