Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913753,1250TA37.8% 104.6 / 43.6 75intergenic (‑330/+60)ybgD/gltAputative fimbrial‑like adhesin protein/citrate synthase
Reads supporting (aligned to +/- strand):  ref base T (36/10);  new base A (0/28);  total (36/39)
Fisher's exact test for biased strand distribution p-value = 2.55e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

ATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATT  >  NC_000913/752977‑753275
                                                                                                                                                         |                                                                                                                                                       
aTTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATaa                                                                                                                                                       <  1:729585/151‑1 (MQ=255)
 ttttCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGATAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATaa                                                                                                                                                       >  1:941809/1‑150 (MQ=255)
   ttCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATaaaa                                                                                                                                                     <  1:811168/150‑1 (MQ=255)
     cATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATc                                                                                                                                                 >  2:91598/1‑151 (MQ=255)
      aTAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCa                                                                                                                                                <  2:1258151/151‑1 (MQ=255)
                    gCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGa                                                                                                                                   >  2:93805/1‑150 (MQ=255)
                     ctaCTACAATTAAGTTTTTAATTGAAATTATGGGTTTGCTAATATTTTGAGCTTTTTTTTTTTTTTTAAATTTTTATATAGTTTTGTTAATGGTTTTTTTTTTTACAATTTTTTGTGTTTAGGCATAAAAAAACAAACCCGCCATATGaa                                                                                                                                  <  1:1419027/148‑1 (MQ=255)
                         tACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAAcggcg                                                                                                                             >  2:741500/1‑151 (MQ=255)
                            cATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCggg                                                                                                                           >  2:933864/1‑150 (MQ=255)
                             aTGAATTTTTTATTTGAAGGTATTGGTTTGGTAAAATTTTGAGCTTTTTTATTCTTTTAAAATATCTTTATATCTTTATAAATGGATTTTATTTTTAAATTTTTTGTTTTTAAGGCATATAAAACCAACCCCGCCATATGAACGGCGGGtt                                                                                                                         <  1:600570/151‑1 (MQ=255)
                              tGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGtt                                                                                                                         >  2:1095888/1‑150 (MQ=255)
                                atttttttAATTAAAAGTTTTGGTTTGTTAATATTTTGAGTTTTTTTTTTTTTTTAAAATATCTATATAGTTTTTTAAATGGTTTTTTTTTTTACAATTTTTTGTGTTTAGCAATAAAAAAACCACCCCGCCATATGAAcggc                                                                                                                              <  1:137152/140‑1 (MQ=255)
                                atgTTTTTATTTAAAATTTTTGGGTGGTTGTTAATTTGCGTTGTTTTTTTCTTTTTAAATTTCTAAAAAGTTTTTTAATGATATTTTTTTTTTAAATTTTTTTGTGTTTGATCATAAAAAAAACAACCCGCCATATGAACGGCGGGTTa                                                                                                                        <  1:691887/147‑1 (MQ=255)
                                   tttttATTTAAAGTTATTGGGTGGTTGATATTTTGGGTTTTTTTATTCTTTTTAAATATTTATATAGGTCTGTAAATGTATTTTTTTTTTAAAATTTTTTGTGTTTAGGCATATAAAAAACACCCCGCCATATGAACGGCGGGtt                                                                                                                         <  1:558618/145‑1 (MQ=255)
                                   tttttAATTGAAAGTATTGGTGTTCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTATTGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCggg                                                                                                                           >  1:560795/1‑143 (MQ=255)
                                   tttttAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGGTaaaa                                                                                                                     >  1:1188812/1‑149 (MQ=255)
                                   tttttAATTAAAAGTTTTGGGTTGCTAATAATTTGAGCTTTTCTTTTTTTTTAAAATATTTATATAGTTCTTTTAATGGTTTTTATTTTTACAATTTTTTGTGTTAAGGAAAAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAat                                                                                                                    <  1:1407726/150‑1 (MQ=255)
                                             aaaGTATTGGGTTGCTGATATTTTGAGCTTTTCTTTTTTTTTTAAATATCTATATAGGTCTTTTAATGGATTTTATTTTTACAATTTTTTGTTTTTAGGAATAAAAAAACCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACtt                                                                                                         <  2:949374/151‑1 (MQ=255)
                                              aaGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTTAAATAttttaaaatt                                                                                                         >  1:343205/1‑143 (MQ=255)
                                               aatATGGGTTTTCTTATATTTTTTGCTTTTCTTTTCTTTTTAAATTTATTTAAAGGTTGTTAAATGATTTTTTTTTTTTTATTTTTTTTTTTTGATAAATATAAAAATAACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACt                                                                                                          <  1:737133/146‑1 (MQ=255)
                                                   ttGGGTTGTTGATATTTTGGGTTGTTCTATTTTTTTTAAATTTCTATATAGTTCTGTTAATGGATTTTATTTTTACATTTTTTTGTGTTTAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAg                                                                                                       <  2:1162920/147‑1 (MQ=255)
                                                   ttGGGTTGCTGATATTTTGAGTTTTTCTATTTTTTTTAAATATTTATATAGGTTTTTTAATGGATTTTATTTTTACATTTTTTTGTGTTTAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCaa                                                                                                    <  2:93549/150‑1 (MQ=255)
                                                     gggTGGATGTTAATTTGTGTTTTTCTATTTTTTTTAAATATTTATATAGTTTGTTAAAAGGTTTTTATTTTTACAATTTTTTTTTTTTAGACAAAAAAAAACCACACCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCa                                                                                                 <  1:1303658/151‑1 (MQ=255)
                                                              aTAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAACAATCAACCATTaa                                                                                         >  1:1199037/1‑150 (MQ=255)
                                                               tAATTTGAGTTTTTTTTTTTTTTTAAAATATCTATATAGTTCTTTTAATGGTTTTTTTTTTTAAAATTTTTTGTGTTTAGGCAAAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAAc                                                                                        <  1:114151/150‑1 (MQ=255)
                                                                   ttGAGCTTTTTTTTTTTTTTTAAATATCTATATAGTTCTTTAAATGGTTTTTTTTTTTAAATTTTTTTGTTTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCtt                                                                                    <  2:578791/150‑1 (MQ=255)
                                                                     gttctttttttttCTTTTAAAATTTTTTTATAGTTTTGTTAATGGATTTTTTTTTTAAATTTTTTTGTTTTTAGACATAAAAAAAGCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCtt                                                                                    <  2:761478/137‑1 (MQ=255)
                                                                      agtttttCTTTTTTTTTTAAATATTTATATAGGTTTTTAAATGGTTTTTATTTTAAAAATTTTTTGTTTTTAGGAAAAAAAAAACCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGa                                                                                  <  2:1008309/144‑1 (MQ=255)
                                                                      aGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTg                                                                                   >  1:550491/1‑148 (MQ=255)
                                                                        ttGTTTTTTTTTTTTTAAATATTTATATAGTTTTTTAAATGGTTTTTTTTTTTAAATTTTTTTTTTTTTAGCAAAAAAAAAACCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGatat                                                                               <  1:1265477/149‑1 (MQ=255)
                                                                         ttttCTTTTTTTTTTAAATATTTATATAGGTCTGTTAATGGATTTTTTTTTTACAATTTTTTGTTTTTAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGatat                                                                               <  2:885036/147‑1 (MQ=255)
                                                                           ttCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAACAATCAACCATTAAACCCTGATATCg                                                                             >  1:460469/1‑149 (MQ=255)
                                                                              ttttttttttAAATTTCTATATTGTTTTGTTAATGTTTTTGTTTTTAAAATTTTTTTGTTTTAAGCAAAAAAAAATTCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCtt                                                                          <  1:893559/144‑1 (MQ=255)
                                                                              tcttCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAAGCATATAAAAATCAACCCACAATATGCACGGCGGGTTAAATTATTTTCCATTTTAAATTCAACCATAAACCCTTGGTATCCCttt                                                                         >  1:207335/3‑150 (MQ=255)
                                                                              tATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCtttt                                                                        >  2:1103860/1‑151 (MQ=255)
                                                                              tATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAAACATTAACGCTTGATATCGCtttt                                                                        >  1:586204/1‑151 (MQ=255)
                                                                               aTTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCtttt                                                                        >  2:1395801/1‑150 (MQ=255)
                                                                                 tgttttaAAATATTAATATAGGGTTTATAATGAATTTTATTTTTAAATTTTTTTGTTTTTAGCAAAAAAAAAACCAAACCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCtttt                                                                        <  1:466727/146‑1 (MQ=255)
                                                                                 tCTTTTTAAATATCTATATAGGTCTTTTAATGGTTTTTATTTTTAAATTTTTTTGTGTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTaa                                                                      <  2:705147/150‑1 (MQ=255)
                                                                                      taaaaTATTTATATAGGTTTGTAAATGGTTTTTATTTTTACATTTTTTTGTGTTTAGGAATAAAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTcg                                                                 <  2:153603/148‑1 (MQ=255)
                                                                                        aaaTATCTATATAGGTTTTTAAATGTTTTTTTTTTTTAAATTTTTTTGTTTTAGGTAATATAAAAACCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTcg                                                                 <  1:748471/148‑1 (MQ=255)
                                                                                            atttatataGGTCTTTTAAGGGATTTTTTTTTTAAATTTTTTGGTGTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGtttt                                                           <  2:1193913/147‑1 (MQ=255)
                                                                                               tatataGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCa                                                        >  2:1395480/1‑150 (MQ=255)
                                                                                                  ataGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTAAAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCata                                                      >  2:70103/1‑149 (MQ=255)
                                                                                                  ataGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAAAATTTACAACTTAGCAATCAAACATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTtata                                                      >  2:459596/1‑149 (MQ=255)
                                                                                                    agttgttttAATGGTTTTTTTTTTTACAATTTTTTTTTTTAAGGCATAAAAAAATAACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCa                                                        <  1:649852/138‑1 (MQ=255)
                                                                                                     ggTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTTAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTTAAGTCGCGTTTTTTATATcc                                                   >  2:982/1‑149 (MQ=255)
                                                                                                     ggTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTg                                                 >  2:668220/1‑151 (MQ=255)
                                                                                                     ggTCTGTTAATGAATTTTATTTTTAAATTTTTTTGTGTTTAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTg                                                 <  2:261536/151‑1 (MQ=255)
                                                                                                       ttttttaaaTGTATTTTATTTTTATATTTTTTTTTTTTAAGTTTTATAAAAAACACCCCGCCATATGAACGGCGGGTTAAAATATTTACATCTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGttttt                                                          <  1:1095818/134‑1 (MQ=255)
                                                                                                       gttctttaaatgatttttttttttACATTTTTTTGTTTTAAGCAATAAAAAATTAACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATcc                                                   <  1:203540/135‑1 (MQ=255)
                                                                                                        cTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGGTATCGCTTTTAAAATCGCGTTTTTCATATcccggat                                              >  2:75007/1‑146 (MQ=255)
                                                                                                         tGTTAATGGATTTTATTTTTACAATTTTTTGGGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGGGAAAATATTTAGAAATTAGCAAAGCAACATTAAAGGGTGGTATCGATTTTTAAAATGCGTTTTTCTGATCCTgga                                               >  2:898121/1‑147 (MQ=255)
                                                                                                         tGTTAATGGATTTTATTTTTAATTTTTTTTGTTTTTGCTTATAGAAAAGTTTTTCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGt                                                <  2:464255/148‑1 (MQ=255)
                                                                                                           ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAAACATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGt                                                >  2:1333282/1‑146 (MQ=255)
                                                                                                              aaTGGTTTTTTTTTTTATAATTTTTTTTGTTTAGACATAAAAAAAACAAGCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGc                                         <  1:421664/151‑1 (MQ=255)
                                                                                                              aaTGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAAACATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTTTATAAAGc                                         >  2:1212531/1‑151 (MQ=255)
                                                                                                              aaTGAATTTTATTTTTACAATTTTTTGTGTTTAGCCATATAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAg                                          <  2:298879/150‑1 (MQ=255)
                                                                                                               aaggaattttattttAACATTTTTTTTTTTTTGGCAATAAAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACa                                           <  2:1034305/144‑1 (MQ=255)
                                                                                                                tGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTTAAATATTTACAACTTAGCAATCAACCATTAACACTTGATATCGATTTTAAAGACGCGTTTTTTATATCCCGTATAAAGCt                                        >  2:863120/1‑150 (MQ=255)
                                                                                                                 ggATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTTAAATAATTACAACTTTGCAATCCACCCTTAACGCTTTATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGtat                                              >  2:1229416/1‑143 (MQ=255)
                                                                                                                      ttttttttTATAATTTTTATGTATTACAGAATAAAAAAATCATCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGAc                                     <  1:1237004/149‑1 (MQ=255)
                                                                                                                      ttttttttTAAAATTTTTTTTGTTTAGAAAAAAAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGAcgc                                   <  1:588662/150‑1 (MQ=255)
                                                                                                                         attttTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCAATAACGCTTGATATTGCTTTTTAAGTCGCCTTTTTTATATCCTGTATACAGCTGAcgcg                                  >  1:416917/1‑148 (MQ=255)
                                                                                                                          tttttACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACATAGCAATCAACCATTAAAGCCTGATATCGCTTTTAAAGTCGCGTTTTTCATATACTGTATACAGCTGACGCGGAc                               >  1:929335/1‑150 (MQ=255)
                                                                                                                          tttttACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAAATTAACAATCAACCATTAAAGCTTGATATCGCTTTTAAAATCGCGTTTTTTATATCCTGTATACAGCTGACGCGGAc                               >  1:270441/1‑150 (MQ=255)
                                                                                                                          tttttACAATTTTTTGTGTTAAGGCATATAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGAc                               <  2:96249/150‑1 (MQ=255)
                                                                                                                          tttttAAATTTTTTTTTGTTAAGGAAAAAAAAAATAAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCgg                                 <  1:977043/148‑1 (MQ=255)
                                                                                                                               aaatttttttttttttAGGAATAAAAAAAACAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGc                           <  1:93659/145‑1 (MQ=255)
                                                                                                                                     ttttGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCACTTTTTAAGACCCGTTTTTCATATTCTGTATACAGCTGACGCGGACGGGCAATCtt                     >  2:479524/1‑149 (MQ=255)
                                                                                                                                             tttAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTTAAATATTTACAACTTAGCAATCAAACATTAAAGCTTGATATCGCTTTTTAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGGAGGGCAAACTTCATACCGt             >  1:1389566/1‑150 (MQ=255)
                                                                                                                                               tAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTTAAATATTTAAAAATTAGCAATCAAACCCTAACGCTTGATATCGCCTTTAAAGTCGCGTTTTTTATTTTCTGTATAAAGCTGACTCCGCCGGGCAATCTTCATTCCGtta           >  1:662563/1‑148 (MQ=255)
                                                                                                                                               tAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCAc          >  2:606253/1‑151 (MQ=255)
                                                                                                                                               tAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGGTAAAATATTTACAAATTAGCAATCAACCATTAAAGCTTGATATCGCTTTTAAAATCGCGTTTTTCCTATCCTGTATACAGCGGGCGCGGACGGGCAATCTTCATACCGTCAc          >  2:229885/1‑151 (MQ=255)
                                                                                                                                                   gCATATAAAAATCAACCCGCCATATGAACGGCGGGTTTAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAATCGCGTTTTTCATATCCTTTATACAGCTGACGCGGACGGGCAATCTTTCTTCCGTCACcgt       >  1:20680/1‑151 (MQ=255)
                                                                                                                                                         taAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCAtt  >  2:741274/1‑151 (MQ=255)
                                                                                                                                                         |                                                                                                                                                       
ATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATT  >  NC_000913/752977‑753275

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