Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091315,4100TG39.8% 168.7 / 23.1 89intergenic (+112/‑35)dnaJ/insL1chaperone Hsp40, DnaK co‑chaperone/IS186 transposase
Reads supporting (aligned to +/- strand):  ref base T (3/50);  new base G (35/0);  total (39/50)
Fisher's exact test for biased strand distribution p-value = 2.01e-21
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCT  >  NC_000913/15261‑15534
                                                                                                                                                        |                                                                                                                             
aTGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGtt                                                                                                                              <  2:826349/151‑1 (MQ=255)
  ggTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTtgt                                                                                                                            <  4:433968/151‑1 (MQ=255)
  ggTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTtg                                                                                                                             <  2:227565/150‑1 (MQ=255)
    tgtgAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGCACCCCCCCAAAAGCTTAATTAAAGGGTGTgg                                                                                                                          >  1:1393641/1‑151 (MQ=255)
     gtgAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCACTCCCCCAAACGCCCACTTAAGGGGTGTgg                                                                                                                          >  3:7125/1‑150 (MQ=255)
     gtgAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTgg                                                                                                                          <  2:740992/150‑1 (MQ=255)
     gtgAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGt                                                                                                                         <  2:748916/151‑1 (MQ=255)
     ctgAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTGACCCCCCTGAACGCCGACTTTAGGGGTgggg                                                                                                                          >  1:59586/2‑150 (MQ=255)
       gaagaaGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTAt                                                                                                                      <  2:310808/151‑1 (MQ=255)
        aagaagTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTAtt                                                                                                                     <  1:184213/151‑1 (MQ=255)
        aagaagTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGAGTACTCCCCCAAACGCCAACTTAAGGGGtgtg                                                                                                                           >  2:851479/1‑146 (MQ=255)
         agaagTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCACCCCCCTAAGCGCGAACTTTAGGGGTGTGGttttt                                                                                                                    >  1:790372/1‑150 (MQ=255)
          gaagTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTAc                                                                                                                   <  3:536973/151‑1 (MQ=255)
            agTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGc                                                                                                                 <  1:772287/151‑1 (MQ=255)
                ttttGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGTACTCCCCTAAAAGCCAACTTTATGGGTGGGGGAttttg                                                                                                                  >  4:54815/1‑143 (MQ=255)
                 tttGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGa                                                                                                             <  1:75323/150‑1 (MQ=255)
                 tttGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGAAAACTCCCCCAAGCGGTAACATAAAGGGTGTGGTATTACGcccggt                                                                                                            >  2:1131683/1‑147 (MQ=255)
                   tGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGata                                                                                                           <  4:186590/150‑1 (MQ=255)
                    gacgacCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATg                                                                                                         <  2:837355/151‑1 (MQ=255)
                       gacCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCACTCCCCCAAACGCTAACTTTAGGGGTGTGGGATTTAGCCCGATattatt                                                                                                      >  1:1417399/1‑151 (MQ=255)
                        acCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGAGCACCCCCCTAAACGCTAACTTTAAGGGTGTGGGATTAAGCCTGGTATGAttt                                                                                                     >  2:1419102/1‑151 (MQ=255)
                         ccTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCTCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTa                                                                                                    <  1:1174638/151‑1 (MQ=255)
                          cTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGCACCCCCCCAAGAGCGAACTTAAGGGGTGTGGGATTTTGCCTGATATGAttt                                                                                                     >  1:543263/1‑149 (MQ=255)
                            gACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACg                                                                                                 <  3:108990/151‑1 (MQ=255)
                              cccGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTg                                                                                               <  2:1089350/151‑1 (MQ=255)
                              cccGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTg                                                                                               <  1:796463/151‑1 (MQ=255)
                                 gCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCg                                                                                            <  3:54816/151‑1 (MQ=255)
                                   tAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATTGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCACTCCCCTAAGCGCTAACTTAAGGGGTGTGGTATTAAGCCTGATATGATTTAACGTGCCg                                                                                            >  2:862164/1‑149 (MQ=255)
                                    aaCCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGTGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCATGATATGATTTAACGTGCCGAt                                                                                          <  1:1370731/150‑1 (MQ=255)
                                    aaCCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGTACTCCCATAAGCGCTAACTTAAGGGGTGTGGGATTTAGCCTGATATTAttttaaaag                                                                                               >  1:394026/1‑139 (MQ=37)
                                    aaCCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGa                                                                                           <  1:1358460/149‑1 (MQ=255)
                                    aaCCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGAt                                                                                          <  1:139707/150‑1 (MQ=255)
                                     aCCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGa                                                                                        <  2:18388/151‑1 (MQ=255)
                                       cTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCACTCCCATAAACGCTAACTTAAAGGGTGTGGGATTTAACCCGATATGATTTAAAGTGCCGATGaa                                                                                       >  2:243166/1‑150 (MQ=255)
                                         ccccAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAAtt                                                                                     <  3:341845/150‑1 (MQ=255)
                                         ccccAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAAtt                                                                                     <  1:1353172/150‑1 (MQ=255)
                                         ccccAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTa                                                                                    <  1:1285692/151‑1 (MQ=255)
                                          cccAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAAt                                                                                      <  2:803789/148‑1 (MQ=255)
                                                aGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCAc                                                                             <  2:481734/151‑1 (MQ=255)
                                                 gCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACg                                                                            <  2:1417466/151‑1 (MQ=255)
                                                 gCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACg                                                                            <  1:924714/151‑1 (MQ=255)
                                                   cTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATGTCTCCCCATAGGGGTTACTTAAGGGGTGTGGGGTTAAGCCTGGTATGATTTTAAGCGCCGACGGAATATTCCCACGAt                                                                          >  1:557347/1‑151 (MQ=255)
                                                           gggCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTc                                                                  <  1:1455863/151‑1 (MQ=255)
                                                            agCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGAGATTAATCCCCCAAGAGGTAAACATAAAGGTGCGGGGTTACGCGTGTAAATATTTAAAAAGTCGATGGATTACACTCACGCTaaactgca                                                                  >  4:106284/2‑144 (MQ=255)
                                                               aGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTGACTGCCATAAAAGCGAACTTAAGGGGTGTGGGGTTACGCCCGGGTTGGTTTTAAGTGCCGATGAATTACTCTTACGGTAAATGGTCAGc                                                               >  2:404734/1‑150 (MQ=255)
                                                               aGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGGAACCCCATAAGGGCCAACCTAAGGGGTGGGGGATTTACGCCGGTATGGTTTAAACTGGCGGAGAAAATCTTCCCCGATAACTGGGCAGc                                                               >  2:669239/1‑150 (MQ=255)
                                                               aGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGAACACTCCCATAAACGCCCAATTAAGGGGTGTGGTGTTAAACCTGATAAGATTTTACGGGGCGATGAAATAATCTCACggaaaattgtcgg                                                                >  2:948458/1‑137 (MQ=255)
                                                                ggCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACCCCCCTAAACGCTAACTTAAGGGGTGTGGGATTAAGCCTGGTATGATTTAAAGTGCCGATGAATTAATCTCACGATAAATgg                                                                    >  3:193646/1‑144 (MQ=255)
                                                                     gggTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTc                                                          <  2:334826/149‑1 (MQ=255)
                                                                      ggTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTgg                                                       <  2:1543421/151‑1 (MQ=255)
                                                                      ggTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTgg                                                       <  1:409683/151‑1 (MQ=255)
                                                                          aaaaTAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTGACTCCCATAAACGCTAACTTTAAGGGTGGGGGATTACGCCTGATATGATTTAAAGTGGCGATGAATTATTCTCACGATAAATGGTTAGGAATTCTGGcccc                                                   >  2:924768/1‑150 (MQ=255)
                                                                           aaaTAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACATAAAGGGTGTGGTATTACGCCTGGTATTATTTAAAGTGCCGATGAAATACTCTCACGATAAATGGTCAGCAATTTTGGccc                                                    >  4:243918/1‑149 (MQ=255)
                                                                             aTAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGCACTCCCATAAACGCTAACTTAAGGGGTGTGGGGTTAAGCCCGGTATTATTTTAACTGCCGATGAATTAATTTCACGATAAATGGGCAGGAATTCttggcc                                                    >  2:1256521/1‑142 (MQ=255)
                                                                               aGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTg                                              <  2:135452/151‑1 (MQ=255)
                                                                                tggTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCACTCCCATAAGCGGCAAACTTAGGGGTGTGGGATTAAGCCCGATATGATTTAAAGTGCCGATGAATTACTTTCACGATAAATGGTTAGCAATTTTGGCCCCTAtt                                               >  4:76461/2‑149 (MQ=255)
                                                                                gggTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTgg                                             <  1:433312/151‑1 (MQ=255)
                                                                                 ggTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGt                                            <  1:1694154/151‑1 (MQ=255)
                                                                                  gTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGAGATTCCCTCCCAAAACGCCAAATTTACAGGTGTGGGATTTCGTACGGTATGATTTAACGTGACGGCGAAATAATCTTCCGATAAATTGTAAGCAATTCTgtaccctcttgg                                             >  1:757748/1‑138 (MQ=255)
                                                                                   tGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGACCATCCCCTAAACGCTAAATTAAAGGGTTTGGGATTTACCCTGGTATTATTTTAAGTGCCGGTGAATTACTCTCACCCTAAATGGGCCGCAAATCTAGCCCCTATTTGt                                            >  1:1589066/1‑149 (MQ=255)
                                                                                     cGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGc                                        <  2:184134/151‑1 (MQ=255)
                                                                                         gAAGATATGCGAGCACCTGTAAAGTGGCGGGGGACACCCCCCTAAAAGCCAACTTAAGGGGTGTGGTATTACGCCTGATATTATTTAAAGTGCCGATGAATTACTCTCACGAAAAATAGTCAGCAAATTTTGCCCATTTTTGTAAGacc                                      >  2:1543482/1‑149 (MQ=17)
                                                                                         gAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCCTAAACGCTAACTTAAAGGGTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAAATGGTCAGCAATTTTGGCCCCTATTTGTAAGccc                                      >  2:365086/1‑149 (MQ=17)
                                                                                          aaGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCACAAACGCTAACTTAAAGGGTGTGGGGTTTCGCCTGATATGATTTAACGTGCCGCCGAAATAATCTCCCCATAAATGGTCAGCAATTCTTGCCCCTATTTGTAAAccc                                      >  4:593940/1‑148 (MQ=17)
                                                                                          aaGATATGCGAGCACCTGTAAAGTGGCGGGGATCAATCCCCTAAGCCCCAAATTTAGGGGTGTGGGATTTTGCCTGATAAGATTTAAAGTGCCGGTGAATTAATCTTACGATAAATAGTCCGCAATTCTTGGCCATATTGGTaaagcc                                      >  3:530016/1‑144 (MQ=255)
                                                                                           aGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTTCCCTAAGCACGCACATAAGGGGGGTGGGATTAAGGCTGTTATGATTTTACAAGCCGGTGAATTACTCTCCCTATAAAAGGTGAGCAATTTTTTCCCCTAATTGTAAGacc                                      >  3:121798/1‑147 (MQ=255)
                                                                                           aGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCgaag                                  <  4:334766/151‑1 (MQ=255)
                                                                                            gATATGCGAGCACCTGTAAAGTGGCGGGGAGTACACCCCCAAGCGCCTAATTCAGGGGTGGGGGGTTATCCCTTGTATGAGTTAACGCGCCGAAGAATTTATCTCACTACAACACGTCTGCAAGTCTGGTCCCTATTgcgtagtgcacca                                  >  3:581843/1‑138 (MQ=255)
                                                                                             atatGCGAGCACCTGTAAAGTGGCGGGGGTTACTCCCCTAAACGCCAACTTAAAGGGTGTGGGGTTTCCCCTGATATTATTTTAAGTGCCGGTGAAATAATCTCACGATAACTTGTCAGCAAATCTTGCCCCTATTTGGAAGACCgaa                                   >  2:1686576/1‑148 (MQ=255)
                                                                                             atatGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCgaagaa                                <  1:964629/151‑1 (MQ=255)
                                                                                               atGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAAc                               <  4:384172/150‑1 (MQ=255)
                                                                                               atGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACt                              <  1:1561499/151‑1 (MQ=255)
                                                                                                tGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCCTAAGCGCTTAATTAAGGGGTGTGGGATTATGCCTGGTATTATTTTACATGCCGATGAATTACTTTCACGATAAATTGTCAGCAATTCTTGCCCATATTTGTAAGCCCggagaa                                >  2:1452242/1‑148 (MQ=255)
                                                                                                tGCGAGCACCTGTAAAGTGGCGGGGATCACCCCCCTAAGCGCTAACTTTAAGGGTGGGGGGTTAAGCCTGATATTAATTTAAGTGCCGGTGAATTAATCTCACGATAAATTGGCAGCAAATTTGGGCCCTATTTGTAAAACCgaagaa                                >  4:451512/1‑148 (MQ=255)
                                                                                                 gCGAGCACCTGTAAAGTGGCGGGGATCACTCCCCTAAACGCTAACTTAAAGGGTGTGGGATTACGCCTGATATGATTTTACGTGCCGATGAATTACTCTCACGATAAATGGTCAGCAATTCTTGCCCCTATTTGTAAGCCCGAAGAACTgg                            >  2:1564771/1‑151 (MQ=17)
                                                                                                 gCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTgg                            <  2:170319/151‑1 (MQ=255)
                                                                                                     gCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATAc                        <  1:1609414/151‑1 (MQ=255)
                                                                                                        ccTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACtt                      <  2:459354/150‑1 (MQ=255)
                                                                                                        ccTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTc                     <  3:209243/151‑1 (MQ=255)
                                                                                                             aaaGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCAc                <  2:1296883/151‑1 (MQ=39)
                                                                                                                gTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTa             <  1:1527185/151‑1 (MQ=37)
                                                                                                                gTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTa             <  2:1273469/151‑1 (MQ=37)
                                                                                                                   gCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATg          <  2:1208264/151‑1 (MQ=34)
                                                                                                                   gCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATg          <  4:517447/151‑1 (MQ=34)
                                                                                                                    cGGTGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTAAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGc         >  1:66559/1‑151 (MQ=18)
                                                                                                                      gggATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCg       <  4:521751/151‑1 (MQ=33)
                                                                                                                       ggATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCgg      >  1:1453868/1‑151 (MQ=33)
                                                                                                                           tCACTCCCATAAGCGCTAACTTAAGGGGTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCgggg    >  1:908042/1‑150 (MQ=17)
                                                                                                                            cACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGct  <  1:243126/151‑1 (MQ=31)
                                                                                                                            cACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGct  >  2:1105355/1‑151 (MQ=31)
                                                                                                                                                        |                                                                                                                             
ATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCT  >  NC_000913/15261‑15534

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: