Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913986,5690TA26.6% 475.4 / 190.5 259P138P (CCA→CCT) ompFouter membrane porin 1a (Ia;b;F)
Reads supporting (aligned to +/- strand):  ref base T (105/83);  new base A (38/30);  total (144/115)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.34e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 271 total aligned reads displayed.

ATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCC  >  NC_000913/986470‑986668
                                                                                                   |                                                                                                   
 tCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTg                                                                                                    <  1:1249750/100‑1 (MQ=255)
   aaCCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGgc                                                                                                  <  2:4305062/100‑1 (MQ=255)
   aaCCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGgc                                                                                                  >  1:803484/1‑100 (MQ=255)
   aaCCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGgc                                                                                                  <  2:894364/100‑1 (MQ=255)
   aaCCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGgc                                                                                                  >  2:3770511/1‑100 (MQ=255)
    aCCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGgca                                                                                                 >  1:3353403/1‑100 (MQ=255)
    aCCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGgca                                                                                                 <  1:91356/100‑1 (MQ=255)
    aCCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGgca                                                                                                 >  2:1225594/1‑100 (MQ=255)
     ccAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGgcag                                                                                                <  2:1966777/100‑1 (MQ=255)
       aGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACTAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGgcagca                                                                                              <  2:351108/100‑1 (MQ=255)
        gACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCat                                                                                             >  1:2775652/1‑100 (MQ=255)
           cAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATc                                                                                          <  2:2041740/100‑1 (MQ=255)
            aaaGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCg                                                                                         >  2:3732023/1‑100 (MQ=255)
             aagaagTTGGAGTTACGATAGGTAGCAACGCCGGCAACACGACGAACGAAGAAGTCATCGCTGTATGCAGTATAACCACCAAATTCTGGCAGCATATTgg                                                                                        >  2:3586332/1‑100 (MQ=255)
             aagaagTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCCGGCAGCATATCgg                                                                                        <  2:1788286/100‑1 (MQ=255)
             aagaagTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCgg                                                                                        <  2:1192116/100‑1 (MQ=255)
             aagaagTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCgg                                                                                        >  2:2647548/1‑100 (MQ=255)
             aagaagTGGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCgg                                                                                        <  2:3042634/100‑1 (MQ=255)
               gaagTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGgtg                                                                                      <  1:3717891/100‑1 (MQ=255)
                aagTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGgtgt                                                                                     <  1:3458215/100‑1 (MQ=255)
                  gTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTaa                                                                                   >  2:3101815/1‑100 (MQ=255)
                   ttGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAAc                                                                                  <  2:2346220/100‑1 (MQ=38)
                   ttGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAAc                                                                                  >  1:2030508/1‑100 (MQ=255)
                   ttGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAAc                                                                                  <  2:2854778/100‑1 (MQ=255)
                     ggAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAAccc                                                                                >  1:1705010/1‑100 (MQ=38)
                     ggAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAAccc                                                                                >  1:3849396/1‑100 (MQ=38)
                     ggAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAAccc                                                                                <  1:689785/100‑1 (MQ=255)
                       aGTTACGTTAGGTAGCAGCGCCGCCAACACGACCAACTAAGAAGTCATCGCTTTAGTCAGTATCACCACCAAATTCTGTCAGCATATCGGTGTAACCCAg                                                                              <  2:4141976/100‑1 (MQ=255)
                       aGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAg                                                                              <  1:2049210/100‑1 (MQ=255)
                         ttACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTg                                                                            <  1:2492369/100‑1 (MQ=255)
                         ttACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTg                                                                            <  2:777262/100‑1 (MQ=255)
                            cGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGcat                                                                         >  1:2998672/1‑100 (MQ=255)
                               tAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCCAACTCCGGCAGCATATCGGTGTAACCCAGTGcatcat                                                                      >  1:3587080/1‑100 (MQ=255)
                                 ggTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATaa                                                                    >  2:57430/1‑100 (MQ=255)
                                 ggTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATaa                                                                    >  1:1403699/1‑100 (MQ=255)
                                  gTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATaaa                                                                   >  2:2535472/1‑100 (MQ=255)
                                   tAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAAc                                                                  <  2:2507921/100‑1 (MQ=38)
                                   tAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAAc                                                                  >  2:79397/1‑100 (MQ=255)
                                    aGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAAcc                                                                 <  2:2228267/100‑1 (MQ=38)
                                     gCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACca                                                                <  2:1638330/100‑1 (MQ=255)
                                         cgccgcCAACTCGACCAACGAAGAAGTCATCGCTGTATGCCGTGCCACCGCCAAACTCAGGCAGCCTATCGGTGTAACCCAGTGCATCATAAACCCCAcc                                                            >  1:3865817/1‑100 (MQ=255)
                                          gtcgccAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCg                                                           >  1:1494640/3‑100 (MQ=38)
                                          gccgccAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCg                                                           <  2:1812464/100‑1 (MQ=38)
                                          gccgccAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCg                                                           >  1:2418371/1‑100 (MQ=255)
                                           ccgccAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGt                                                          >  1:3144017/1‑100 (MQ=255)
                                            cgccAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTa                                                         <  1:3271786/100‑1 (MQ=38)
                                              ccAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGt                                                       >  1:259168/1‑100 (MQ=255)
                                                aaCACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTa                                                     <  1:1957861/100‑1 (MQ=255)
                                                  cacGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACg                                                   <  1:3646145/100‑1 (MQ=255)
                                                   acGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAAGTCTGGCAGCATAGCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACgg                                                  <  2:3232877/100‑1 (MQ=255)
                                                   acGACCAACGAAGAAGTCATCGCTATATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACgg                                                  <  2:3962164/100‑1 (MQ=37)
                                                    cGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGc                                                 >  1:2195086/1‑100 (MQ=255)
                                                     gACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGcc                                                <  2:2134676/100‑1 (MQ=255)
                                                      aCCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCg                                               <  1:1690529/100‑1 (MQ=255)
                                                       ncAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGt                                              >  1:1028911/2‑100 (MQ=255)
                                                       ccAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGt                                              >  1:3043974/1‑100 (MQ=255)
                                                        cAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTa                                             >  1:4000671/1‑100 (MQ=255)
                                                         aaCGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTaa                                            >  1:1895015/1‑100 (MQ=255)
                                                          aCGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAAt                                           <  1:847059/100‑1 (MQ=255)
                                                          aCGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAAt                                           <  2:1444087/100‑1 (MQ=255)
                                                           cGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATc                                          <  2:4335257/100‑1 (MQ=255)
                                                              agaagTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGaa                                       >  2:2070584/1‑100 (MQ=255)
                                                               gaagTCATCGCTGGATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTATTTACGGCCGTATTCGaaa                                      <  1:3675386/100‑1 (MQ=255)
                                                                   tCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAAc                                  <  2:1352169/100‑1 (MQ=255)
                                                                    cATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAAcc                                 >  1:1309452/1‑100 (MQ=255)
                                                                     aTCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCa                                >  1:2402084/1‑100 (MQ=38)
                                                                     aTCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCa                                >  1:1845839/1‑100 (MQ=255)
                                                                      tCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCaa                               >  1:1625675/1‑100 (MQ=255)
                                                                        gatGTATGCAGTATCACCACCAAATTCTGGCAGCTTATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACg                             >  2:4009594/3‑100 (MQ=255)
                                                                        gCTGTATGCAGTATCACCACCAAATTCTGGCTGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACg                             >  1:4252860/1‑100 (MQ=255)
                                                                        gCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACg                             <  1:651386/100‑1 (MQ=255)
                                                                        gCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACg                             <  2:2711496/100‑1 (MQ=255)
                                                                            tATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGTCGTAATCGAAAGAACCAACGTCAg                         >  2:3189999/1‑100 (MQ=255)
                                                                            tATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAg                         <  2:3841071/100‑1 (MQ=255)
                                                                              tGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACTCCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCg                       <  2:1516338/100‑1 (MQ=255)
                                                                              tGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCg                       >  2:4220896/1‑100 (MQ=255)
                                                                              cgatctATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCg                       <  2:1057624/95‑1 (MQ=255)
                                                                               gCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGt                      >  1:435770/1‑100 (MQ=255)
                                                                                ccgtgtCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTa                     >  1:442487/6‑100 (MQ=38)
                                                                                cAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTa                     <  2:1683629/100‑1 (MQ=255)
                                                                                cAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTa                     >  1:3229656/1‑100 (MQ=255)
                                                                                  ntgtCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGGAGTTACGGCCTTAATCGAAAGAACCAACGTCAGCGTAtt                   >  1:3531483/4‑100 (MQ=38)
                                                                                    aTCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTa                 <  1:2665604/100‑1 (MQ=255)
                                                                                    aTCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTa                 >  2:2587587/1‑100 (MQ=255)
                                                                                     ccaccgcCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTaa                >  1:3911745/2‑100 (MQ=255)
                                                                                      caccacCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAg               <  2:3806995/100‑1 (MQ=255)
                                                                                       accgccAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGa              >  1:2444704/1‑100 (MQ=255)
                                                                                        ccgccaAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGAc             >  1:1878420/4‑100 (MQ=255)
                                                                                        ccaccaAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGAc             >  2:544662/1‑100 (MQ=255)
                                                                                         ngccaAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGAcc            >  1:557707/3‑100 (MQ=255)
                                                                                          gccaAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGAccc           >  1:2571693/2‑100 (MQ=255)
                                                                                          accaAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGAccc           >  2:2453600/1‑100 (MQ=255)
                                                                                           ccaAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCcg          >  1:1157384/1‑100 (MQ=255)
                                                                                           ccaAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGAccca          <  2:817390/100‑2 (MQ=255)
                                                                                            caAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCcgc         >  1:3738331/1‑100 (MQ=255)
                                                                                                ttCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAAt     >  2:483266/1‑100 (MQ=255)
                                                                                                 tCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATg    <  2:3668302/100‑1 (MQ=255)
                                                                                                  cTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGc   >  2:2009411/1‑100 (MQ=255)
                                                                                                  cTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGACAGCGAAAATAAGACCCACGAATGc   >  2:865735/1‑100 (MQ=255)
                                                                                                   tGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGcc  >  2:2551278/1‑100 (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: ‑