Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913986,5600AG25.7% 436.8 / 165.7 263G141G (GGT→GGC) ompFouter membrane porin 1a (Ia;b;F)
Reads supporting (aligned to +/- strand):  ref base A (104/90);  new base G (31/36);  total (136/127)
Fisher's exact test for biased strand distribution p-value = 3.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.60e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 295 total aligned reads displayed.

GTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCC  >  NC_000913/986461‑986659
                                                                                                   |                                                                                                   
gTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATcacca                                                                                                     >  1:3460875/1‑100 (MQ=255)
   cAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCAccgcca                                                                                                  >  1:4840316/1‑100 (MQ=38)
   cAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCAccacca                                                                                                  <  1:1262362/100‑1 (MQ=255)
   cAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCAccacca                                                                                                  <  2:615732/100‑1 (MQ=255)
      tccATTAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCaaac                                                                                               <  2:4729379/99‑2 (MQ=255)
      gCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCaaac                                                                                               >  1:2012209/1‑99 (MQ=39)
        cATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGTCGTATCACCACCAAATTc                                                                                             >  2:4272764/1‑100 (MQ=255)
          tCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTcag                                                                                           <  2:3408221/100‑3 (MQ=38)
           cAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTgg                                                                                          <  2:262405/100‑1 (MQ=255)
            aaCCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGgc                                                                                         <  1:367203/100‑1 (MQ=255)
            aaCCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGgc                                                                                         >  1:2981211/1‑100 (MQ=255)
              ccAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGACGCAAACCCGACCAACGAAGAAGTCATAGCTGTATGAAGTAACACCCCCAAATCCTAgcag                                                                                       >  1:2573308/1‑100 (MQ=255)
               cAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGgcagc                                                                                      <  2:4835860/100‑1 (MQ=255)
               cAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGgcagc                                                                                      <  2:1591660/100‑1 (MQ=255)
               cAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGgcagc                                                                                      <  2:2894259/100‑1 (MQ=255)
               cAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGgcagc                                                                                      <  1:2080502/100‑1 (MQ=255)
                aGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGgcagca                                                                                     <  2:3599751/100‑1 (MQ=255)
                 gACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCat                                                                                    >  1:4538096/1‑100 (MQ=255)
                    cAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATc                                                                                 >  1:3350114/1‑100 (MQ=255)
                    cAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATc                                                                                 >  2:4916879/1‑100 (MQ=255)
                    cAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATc                                                                                 >  1:2918225/1‑100 (MQ=255)
                      aagaagTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCgg                                                                               >  1:2182441/1‑100 (MQ=255)
                      aagaagTTGGAGTTACGATAGGGAGCAACGCCGCCAACACGACCGACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCgg                                                                               <  1:3086033/100‑1 (MQ=255)
                          agTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTa                                                                           >  1:2805213/1‑100 (MQ=255)
                           gTTGGAGTTACGATAGGTAGCAACACCGCGAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAAATCAGGCAGCATATCGGTGTaa                                                                          >  1:4809796/1‑100 (MQ=37)
                              ggATTTACGATAGGTAGCAACGCCTCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAAccc                                                                       <  1:3880727/100‑1 (MQ=255)
                              ggAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAAccc                                                                       >  1:1089139/1‑100 (MQ=38)
                              ggAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAAccc                                                                       >  1:470240/1‑100 (MQ=255)
                               gAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCa                                                                      >  1:2326240/1‑100 (MQ=255)
                                 gTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGt                                                                    >  1:1810530/1‑100 (MQ=38)
                                 gTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGt                                                                    >  1:4327857/1‑100 (MQ=255)
                                 gTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGCATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGt                                                                    >  1:317241/1‑100 (MQ=255)
                                  ttACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTg                                                                   <  2:846744/100‑1 (MQ=255)
                                     cGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACAACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGcat                                                                <  2:1029732/100‑1 (MQ=255)
                                      gATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGcatc                                                               <  2:3854328/100‑1 (MQ=38)
                                      gATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGcatc                                                               <  2:4326143/100‑1 (MQ=255)
                                       aTAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGcatca                                                              <  1:3156801/100‑1 (MQ=255)
                                          ggTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATaa                                                           >  2:1875311/1‑100 (MQ=255)
                                           gTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATaaa                                                          <  2:2404935/100‑1 (MQ=38)
                                           gTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATaaa                                                          <  1:399793/100‑1 (MQ=255)
                                           gTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATaaa                                                          >  1:4071460/1‑100 (MQ=255)
                                             aGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAAcc                                                        >  1:833846/1‑100 (MQ=38)
                                             aGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAAcc                                                        <  2:2961752/100‑1 (MQ=255)
                                             aGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAAcc                                                        >  2:4524215/1‑100 (MQ=255)
                                              gCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACca                                                       >  1:5127651/1‑100 (MQ=255)
                                               cAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGGCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCAGCATAAACcac                                                      <  2:3660210/100‑1 (MQ=255)
                                                 ccgccgcCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCacac                                                    >  1:2412240/2‑100 (MQ=38)
                                                 aCGCCGCCAACACGACTAACAAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTTACCCAGTGCATCATAAACCacac                                                    <  2:4250417/100‑1 (MQ=255)
                                                 aCGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCacac                                                    <  2:313378/100‑1 (MQ=255)
                                                 aCGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCacac                                                    <  1:4598733/100‑1 (MQ=255)
                                                 aCGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCacac                                                    >  1:495999/1‑100 (MQ=255)
                                                 aCGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCacac                                                    >  2:1313238/1‑100 (MQ=255)
                                                  cgccgcCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACAcc                                                   >  1:757181/1‑100 (MQ=38)
                                                  cgccgcCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACAcc                                                   >  1:3758925/1‑100 (MQ=38)
                                                   gccgccAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCg                                                  >  1:1299892/1‑100 (MQ=38)
                                                   gccgccAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCg                                                  <  2:926070/100‑1 (MQ=255)
                                                      gccAACACGACCATCGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAg                                               <  2:3039927/100‑1 (MQ=38)
                                                      gccAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAg                                               <  2:1279324/100‑1 (MQ=38)
                                                      gccAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAg                                               <  2:2076007/100‑1 (MQ=38)
                                                      gccAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAg                                               <  2:73296/100‑1 (MQ=38)
                                                      gccAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAg                                               <  2:1802908/100‑1 (MQ=255)
                                                       ccAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGt                                              >  2:1702929/1‑100 (MQ=255)
                                                         aaCACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTa                                            >  1:965246/1‑100 (MQ=255)
                                                          atctgACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTAc                                           <  1:2333391/96‑1 (MQ=255)
                                                           cacGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACg                                          <  2:3947217/100‑1 (MQ=38)
                                                           cacGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACg                                          >  2:95473/1‑100 (MQ=255)
                                                             cGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCTCCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGc                                        <  1:3375668/100‑1 (MQ=255)
                                                             cGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGc                                        <  1:1350777/100‑1 (MQ=255)
                                                             cGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGc                                        <  1:2196404/100‑1 (MQ=255)
                                                             cGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACAACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGc                                        <  2:4613366/100‑1 (MQ=255)
                                                              gACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTAcagat                                       >  2:2305715/1‑96 (MQ=255)
                                                               aCCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCg                                      <  2:2205118/100‑1 (MQ=255)
                                                               aCCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCg                                      >  1:1928114/1‑100 (MQ=255)
                                                                ccAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGt                                     >  1:1519765/1‑100 (MQ=255)
                                                                ccAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGt                                     >  2:3192971/1‑100 (MQ=255)
                                                                  aaCGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTaa                                   <  2:833460/100‑1 (MQ=255)
                                                                  aaCGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTaa                                   >  2:1599753/1‑100 (MQ=255)
                                                                       agaagTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGaa                              <  2:2394013/100‑1 (MQ=38)
                                                                        gaagTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGaaa                             <  2:1168085/100‑1 (MQ=255)
                                                                         aagTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAg                            <  2:561907/100‑1 (MQ=38)
                                                                         aagTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAg                            <  2:4365734/100‑1 (MQ=255)
                                                                         aagCCATCGCTGTATGCCGTGTCCCCGCCAAACTCAGGCAGCATATCGGTGTAAACCAGTGCACCATAAACCACACCGTAGTTACGGCCGTAAACCAAAg                            <  2:494083/100‑1 (MQ=255)
                                                                             cATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAAcc                        >  2:4065888/1‑100 (MQ=255)
                                                                             cATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAAcc                        <  1:1042725/100‑1 (MQ=255)
                                                                              aTCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCa                       <  2:899087/100‑1 (MQ=38)
                                                                                 gCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACg                    <  2:331774/100‑1 (MQ=255)
                                                                                   tGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTc                  <  2:4155203/100‑1 (MQ=255)
                                                                                      gtGCAGTATCACCACCAAATTCTGGCAGCATATCGTTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGc               >  1:1900629/2‑100 (MQ=255)
                                                                                      gtGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGc               >  1:2428849/2‑100 (MQ=255)
                                                                                       tGCAGTATCACCACCAAATCCTGGCAGCATATCGGTGTAACCGAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCg              >  1:4392824/1‑100 (MQ=255)
                                                                                        gCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGGACCAACGTCAGCGt             >  2:2125084/1‑100 (MQ=255)
                                                                                        gCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGt             >  1:3512374/1‑100 (MQ=255)
                                                                                              tCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTaa       <  1:1757136/100‑1 (MQ=255)
                                                                                               caccacCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAg      >  1:3324572/1‑100 (MQ=255)
                                                                                               caccacCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGATACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAg      >  1:4550249/1‑100 (MQ=255)
                                                                                                 ccaccaAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGAc    <  1:2030837/100‑1 (MQ=255)
                                                                                                  caccaAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGAcc   <  1:58475/100‑1 (MQ=255)
                                                                                                   accaAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGAccc  >  1:4153287/1‑100 (MQ=255)
                                                                                                   accaAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGAccc  >  1:3271020/1‑100 (MQ=255)
                                                                                                   accaAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGAccc  >  1:4675928/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
GTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCC  >  NC_000913/986461‑986659

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