Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913986,5510AC25.1% 564.7 / 204.0 291T144T (ACT→ACG) ompFouter membrane porin 1a (Ia;b;F)
Reads supporting (aligned to +/- strand):  ref base A (108/110);  new base C (31/42);  total (139/152)
Fisher's exact test for biased strand distribution p-value = 3.44e-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 291 total aligned reads displayed.

ACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTA  >  NC_000913/986453‑986649
                                                                                                  |                                                                                                  
 cAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGt                                                                                                  <  2:851753/100‑1 (MQ=255)
  aGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGtg                                                                                                 <  2:4694358/100‑2 (MQ=255)
  aGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTa                                                                                                 >  1:1158965/1‑100 (MQ=255)
    cGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATc                                                                                               <  1:4292509/100‑1 (MQ=255)
    cGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATc                                                                                               <  1:3793181/100‑1 (MQ=255)
     gAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATca                                                                                              <  2:2429136/100‑1 (MQ=255)
     gAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATca                                                                                              <  2:2093890/100‑1 (MQ=255)
     gAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATca                                                                                              <  2:629564/100‑1 (MQ=255)
      aaGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATcac                                                                                             <  1:4633895/100‑1 (MQ=255)
       aGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATcacc                                                                                            <  2:3673043/100‑1 (MQ=255)
          tCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCaccacc                                                                                         <  2:2652691/100‑1 (MQ=255)
          tCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCaccacc                                                                                         <  2:4047244/100‑1 (MQ=255)
           cAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCAccacca                                                                                        <  1:1262362/100‑1 (MQ=255)
             ggCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCaaa                                                                                      >  1:1460832/1‑100 (MQ=255)
              gCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCaaac                                                                                     >  1:2012209/1‑99 (MQ=39)
                cATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTc                                                                                   >  2:4131562/1‑100 (MQ=255)
                 aTCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCt                                                                                  <  1:1541351/100‑1 (MQ=255)
                 aTCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCt                                                                                  >  1:4911958/1‑100 (MQ=255)
                    aaCCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGgc                                                                               <  1:1842910/100‑1 (MQ=255)
                      ccAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGgcag                                                                             >  1:3896564/1‑100 (MQ=255)
                            cAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATc                                                                       >  1:3350114/1‑100 (MQ=255)
                            cAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATc                                                                       >  2:4916879/1‑100 (MQ=255)
                            cAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATc                                                                       <  1:803950/100‑1 (MQ=255)
                              aagaagTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCgg                                                                     >  1:2182441/1‑100 (MQ=255)
                                gaagTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGgtg                                                                   <  1:5026902/100‑1 (MQ=255)
                                  agTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTa                                                                 <  2:3047526/100‑1 (MQ=255)
                                  agTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTa                                                                 >  1:2805213/1‑100 (MQ=255)
                                   gTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTaa                                                                >  1:998224/1‑100 (MQ=38)
                                   gTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTaa                                                                >  2:2973422/1‑100 (MQ=255)
                                      ggAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAAccc                                                             >  1:2955721/1‑100 (MQ=38)
                                      ggAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAAccc                                                             >  1:2101817/1‑100 (MQ=38)
                                       gAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCa                                                            >  1:2326240/1‑100 (MQ=255)
                                        aGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAg                                                           <  2:3618632/100‑1 (MQ=255)
                                         gTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGt                                                          >  1:4750379/1‑100 (MQ=255)
                                         gTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGCATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGt                                                          >  1:317241/1‑100 (MQ=255)
                                             cGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGcat                                                      >  2:1517642/1‑100 (MQ=255)
                                              gATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGGTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGcatc                                                     <  2:1733319/100‑1 (MQ=38)
                                              gATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGcatc                                                     <  2:3854328/100‑1 (MQ=38)
                                              gATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGcatc                                                     <  2:315095/100‑1 (MQ=255)
                                              gATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGcatc                                                     <  2:4326143/100‑1 (MQ=255)
                                               aTAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGcatca                                                    <  1:3156801/100‑1 (MQ=255)
                                                  ggTAGCAACGCCGCGAAGATGACCAACGTTGAAATCTTCGCTGTATGCAGTATCACTATCAAATTCTGCCAGCTTGTCGCCGTGCACCTTTGCATCATaa                                                 >  2:2565089/1‑100 (MQ=255)
                                                  ggTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATaa                                                 <  2:1889541/100‑1 (MQ=38)
                                                  ggTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATaa                                                 >  2:2377030/1‑100 (MQ=255)
                                                  ggTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATaa                                                 >  1:4437577/1‑100 (MQ=255)
                                                   gTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATaaa                                                >  1:4071460/1‑100 (MQ=255)
                                                   gTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATaaa                                                >  1:3314785/1‑100 (MQ=255)
                                                    tAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAAc                                               <  2:3439114/100‑1 (MQ=38)
                                                    tAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAAc                                               <  2:1183588/100‑1 (MQ=255)
                                                     aGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAAcn                                              <  2:3874317/100‑2 (MQ=255)
                                                     aGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAAcc                                              <  1:53536/100‑1 (MQ=255)
                                                     aGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAAcc                                              >  2:4524215/1‑100 (MQ=255)
                                                        caCGCCGCCAACACGATCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACAAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACcaca                                           >  2:4583966/2‑100 (MQ=255)
                                                        caCGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACcaca                                           >  1:940589/2‑100 (MQ=255)
                                                        aaCGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACcaca                                           <  2:3206679/100‑1 (MQ=255)
                                                         ccgccgcCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCacac                                          >  1:2412240/2‑100 (MQ=38)
                                                           gccgccAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCg                                        >  1:456063/1‑100 (MQ=38)
                                                           gccgccAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCg                                        >  1:1299892/1‑100 (MQ=38)
                                                           gccgccAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCg                                        <  2:926070/100‑1 (MQ=255)
                                                           gccgccAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCg                                        >  2:4274339/1‑100 (MQ=255)
                                                           caccccAACACGACCAACGAATAAGACATCGCTGGAGGCAGAACCACCACCAAATTCTATCAGCATATCCGTGTAACCCATTGCATCATAAACAACACCg                                        <  2:292341/96‑1 (MQ=255)
                                                              gccAACACGACCAGCGAAGAAGGCATCTCGGTGTGCCGTGTCACCGCCAAACTCAGGCAGCAGAGCGGTGTAACACAGTGCATCATAAACCACACGGTAg                                     <  2:4252792/100‑1 (MQ=255)
                                                              gccAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAg                                     <  2:73296/100‑1 (MQ=38)
                                                              gccAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAg                                     <  2:1802908/100‑1 (MQ=255)
                                                                  atctgACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTAc                                 <  1:2333391/96‑1 (MQ=255)
                                                                  acacGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTAc                                 >  2:4190924/1‑100 (MQ=255)
                                                                   caggACGAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTTCATCATAAACCACACCGTAGTTACg                                <  2:1115343/97‑1 (MQ=255)
                                                                   cacGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACg                                >  2:95473/1‑100 (MQ=255)
                                                                     cGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCTCCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGc                              <  1:3375668/100‑1 (MQ=255)
                                                                     cGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGc                              <  2:1311175/100‑1 (MQ=255)
                                                                     cGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGc                              <  1:2196404/100‑1 (MQ=255)
                                                                     cGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGc                              >  1:4163298/1‑100 (MQ=255)
                                                                      gACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGcc                             <  2:336321/100‑1 (MQ=38)
                                                                       aCCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCATCATATCGGTGTAAGCCAGTGCATCATAAACCACACCGTAGTTACGGCCg                            >  2:5003971/1‑100 (MQ=255)
                                                                       aCCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCg                            <  2:2205118/100‑1 (MQ=255)
                                                                       aCCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCg                            >  1:1928114/1‑100 (MQ=255)
                                                                        ccAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGt                           >  1:205485/1‑100 (MQ=255)
                                                                          aaCGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTaa                         <  2:833460/100‑1 (MQ=255)
                                                                          aaCGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTaa                         >  2:1595300/1‑100 (MQ=255)
                                                                             gaagaaGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCg                      >  1:3703827/1‑100 (MQ=255)
                                                                               agaagTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGaa                    <  1:2603860/100‑1 (MQ=255)
                                                                               agaacTCATCGCTGTATGCAGTATCACCACCAAATTGTGGCAGCATATCAGTGTAACCCAGTGCATCATGAACCACACCGTAGTTACGGCCGTGATCGaa                    >  2:4773298/1‑100 (MQ=255)
                                                                                gaagTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGaaa                   >  2:579894/1‑100 (MQ=255)
                                                                                gaagTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGaaa                   <  1:3462058/100‑1 (MQ=255)
                                                                                 aagTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCGTCATAAACCACACCGTAGTTACGGCCGTAATCGAAAg                  >  1:4006201/1‑100 (MQ=255)
                                                                                 aagTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAg                  <  2:4365734/100‑1 (MQ=255)
                                                                                  agTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGa                 >  1:3093444/1‑100 (MQ=255)
                                                                                    tCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTGACGGCCGTAATCGAAAGAAc               <  2:2302065/100‑1 (MQ=255)
                                                                                     cATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAAcc              <  1:1042725/100‑1 (MQ=255)
                                                                                      aTCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCa             <  2:899087/100‑1 (MQ=38)
                                                                                        cGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCGCACCGTAGTTACGGCCGTAATCGAAAGAACCAAc           >  1:663093/1‑100 (MQ=255)
                                                                                         gCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACg          >  1:3032904/1‑100 (MQ=255)
                                                                                             tATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAg      <  1:3548993/100‑1 (MQ=255)
                                                                                             gttGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAATCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAg      >  1:4388345/3‑100 (MQ=255)
                                                                                              gtGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGc     >  1:2428849/2‑100 (MQ=255)
                                                                                              aTGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGc     >  1:4826894/1‑100 (MQ=255)
                                                                                               tGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCg    >  2:2404744/1‑100 (MQ=255)
                                                                                               tGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCg    <  2:1680362/100‑1 (MQ=255)
                                                                                                gCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGGACCAACGTCAGCGt   >  2:2125084/1‑100 (MQ=255)
                                                                                                gCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGt   >  1:1768811/1‑100 (MQ=255)
                                                                                                  |                                                                                                  
ACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTA  >  NC_000913/986453‑986649

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