Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913753,1250TA45.6% 126.4 / 50.0 92intergenic (‑330/+60)ybgD/gltAputative fimbrial‑like adhesin protein/citrate synthase
Reads supporting (aligned to +/- strand):  ref base T (43/7);  new base A (0/42);  total (43/49)
Fisher's exact test for biased strand distribution p-value = 3.45e-19
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATT  >  NC_000913/752976‑753275
                                                                                                                                                      |                                                                                                                                                      
cATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCAtata                                                                                                                                                       <  1:746025/151‑1 (MQ=255)
 aTTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATaa                                                                                                                                                      <  1:959905/151‑1 (MQ=255)
  ttttCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATaa                                                                                                                                                      >  1:15128/1‑150 (MQ=255)
   tttCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATaaa                                                                                                                                                     >  1:617067/1‑150 (MQ=255)
      cATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATc                                                                                                                                                 <  1:568550/151‑1 (MQ=255)
       aTAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATc                                                                                                                                                 >  2:21489/1‑150 (MQ=255)
          gggTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAAcc                                                                                                                                             <  2:447360/151‑1 (MQ=255)
            gTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCg                                                                                                                                           >  2:847755/1‑151 (MQ=255)
                taAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCat                                                                                                                                       >  2:1396295/1‑151 (MQ=255)
                taAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCat                                                                                                                                       >  2:790627/1‑151 (MQ=255)
                  aaTTCGACTACCATGAAGTTTTTAATGGAAATTATTGGTTTGCTGATAATTTGATCTTTTCTATTTTTTTTAAATATCTATATAGGTCTTTTAATGGATTTTATTTTTACAATTTTTTGTTTTTAGGCATAAAAAAATCAACCCGCCatat                                                                                                                                     <  2:1120611/151‑1 (MQ=255)
                   aTGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATg                                                                                                                                    >  2:256330/1‑151 (MQ=255)
                    tGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATg                                                                                                                                    >  1:747772/1‑150 (MQ=255)
                      cGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAAc                                                                                                                                 >  2:119024/1‑151 (MQ=255)
                       gACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAAcg                                                                                                                                >  2:321571/1‑151 (MQ=255)
                         cTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAAcggc                                                                                                                              >  1:1188874/1‑151 (MQ=255)
                            ccATGAATTTTTTAATTAAAATTTTTGGTTTGTTAATATTTTGAGTTTTTTTTTTTTTTTAAAATATTTATATGGGTTTGTAAATGGTTTTTTTTTTTAAATTTTTTTGTGTTTAGGAATAAAAAAACCAACCCGCCATATGAAcggcgg                                                                                                                            <  2:494657/150‑1 (MQ=255)
                              aTGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAAcggcg                                                                                                                             >  2:764359/1‑147 (MQ=255)
                              aTGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGtt                                                                                                                         >  2:805340/1‑151 (MQ=255)
                               tGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCggggt                                                                                                                         >  2:778226/1‑148 (MQ=255)
                               tGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGtt                                                                                                                         >  1:487354/1‑150 (MQ=255)
                               tGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGtt                                                                                                                         >  2:51807/1‑150 (MQ=255)
                                 aattttttAATTGAAAGTTTTGGGTTGCTGATATTTTGAGTTTTTCTTTTTTTTTAAAATATCTATATAGGTCTGTAAAGGGTTTTTTTTTTTACATTTTTTTGTTTTTAGGCATAAAAAAACCAACCCGCCATATGAACGGCGGGt                                                                                                                          <  2:1364256/144‑1 (MQ=255)
                                    ttttaATTTAAAAGTATGGGTTGCGTTATATTTTGAGTTTTTTTATTCTTTTAAAATATTTATATAGATCTATTAATGATTTTTATTTTAAAATTTTTTTGTTTTTAGGCATAAAAAATTAACCCCGCCATATGAACGGCGGGTTaaaa                                                                                                                     <  2:171192/149‑1 (MQ=255)
                                      tttAATTAAAAGTATTGGTTTGCTGATATTTTGAGTTTTTCTTTTTTTTTTAAATATTTATATAGGTCTGTTAATGTTTTTTTTTTTTACATTTTTTTGTTTTTAGACATAAAAAAATCACCCCGCCATATGAACGGCGGGTTAAAatat                                                                                                                  <  2:498293/150‑1 (MQ=255)
                                          aTTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATAttt                                                                                                                >  1:410914/1‑148 (MQ=255)
                                          aTTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATAttt                                                                                                                >  2:1147393/1‑148 (MQ=255)
                                             gAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTTAAATATTTACaaa                                                                                                           >  2:736659/1‑149 (MQ=255)
                                              aaaaTATTGGGTTGATGATTATTTGAACTGTTTTTTTTTTTTTAAATATTTAAATGGTTCTTTTAATGGATTTTATTTTCACATTTTTTTGTTTTTAGACAAAAAAAAAAAACCCCGCCATATGAACGGCGGGTTAAAATATTTACaa                                                                                                            <  2:322842/148‑1 (MQ=255)
                                              aaaGTTTTGGGTTGCTGATAATTTGAGCTTTTCTTTTTTTTTTAAATATCTATATAGGTCTGTTAATGGTTTTTTTTTTTACATTTTTTTGTGTTTAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAAc                                                                                                           <  1:1356674/149‑1 (MQ=255)
                                              aaaGTATTGGGTTTATAATAATTTAGGTTTTTTTTTTTTTTTTAAATTTCTATATAGGTTTGTTAATGAATTTTATTTTTAAAATTTTTTGTTTTTATACATAAAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTa                                                                                                               <  2:498974/145‑1 (MQ=255)
                                                 gttttGGTTTTCTAATATTTTGGGTTTTTCTATTCTTTTAAAATATTTATATGGGTCTGTTAATGGTTTTTTTTTTTACAATTTTTGTTGTTTAGGCATATAAAAATAACCCCGCCATATGAGCGGCGGGTTAAAATATTTACAACTTa                                                                                                        <  1:1122180/146‑1 (MQ=255)
                                                  tATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTaac                                                                                                      >  2:1006544/1‑148 (MQ=255)
                                                    ttGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTTAAATATTTACAACTTAGCaa                                                                                                    >  2:351977/1‑150 (MQ=255)
                                                      ggtttGCTGATAATTTGCCCTGTTCTATTCTTTTAAAATATCTATATAGTTCTTTTAATGGATTTTATTTTTAAAATTTTTTGTTTTCAGGCATAAAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAAtg                                                                                                  <  1:799022/147‑2 (MQ=255)
                                                      gggTTGTTAATATTTGGAGCTGTTCTATTTTTTTTAAATATCTATATAGTTTTTTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAAt                                                                                                   <  1:860499/149‑1 (MQ=255)
                                                      gggTTGCTGATAATTTGAGTTTTTTTTTTTTTTTAAAATATCTATATAGGTTTGTTAATGGTTTTTATTTTTACATTTTTTTGTTTTTAGGAATAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATc                                                                                                  <  1:1250955/150‑1 (MQ=255)
                                                       ggTTGTTAATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTTTGTTAATGGATTTTATTTTTACAATTTTTTGTTTTTAGGCATAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCaa                                                                                                <  1:795447/151‑1 (MQ=255)
                                                              tGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCAtt                                                                                           >  1:1369218/1‑150 (MQ=255)
                                                                aTAATTTAAGTTGTTTTTTTTTTTTAAAATATTTATATAGTGTTTTTAATGATTTTTATTTTTACAATTTTTTTTGTTTAGCAAAAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTa                                                                                          <  2:101033/149‑1 (MQ=255)
                                                                 tAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTaa                                                                                         >  2:1265408/1‑149 (MQ=255)
                                                                    tttGAGTTGTTTTTTTTTTTTTAAATATCTATATAGGTCTTTTAATGGATTTTATTTTTAAAATTTTTTGTGTTTAGGAAAATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCAt                                                                                            <  1:313042/143‑1 (MQ=255)
                                                                      tGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTg                                                                                   >  1:766080/1‑150 (MQ=255)
                                                                       gggCTGTTTTATTTTTTTTAAATATCTATATAGGTTTGTTAATGGATTTTATTTTTAAAATTTTTTGTGTTTAGGCATAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGa                                                                                  <  1:554045/148‑1 (MQ=255)
                                                                       gAGCTTTTCTATTTTTTTTAAATATCTATATAGGTCTTTTAATGGATTTTATTTTTAAAATTTTTTGTGTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTg                                                                                   <  1:1271997/149‑1 (MQ=255)
                                                                           ttttttttttttttaAAATTTCTATATAGGTCTTTTAATGTTTTTTTTTTTTACAATTTTTTGTTTTTAGGCATATAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATc                                                                              <  2:1150798/140‑1 (MQ=255)
                                                                           ttttcttttcttTTTAAATATCTATATAGGTTTGTTAATGGATTTTATTTTTACATTTTTTTGTGTTTAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATc                                                                              <  1:710610/148‑1 (MQ=255)
                                                                           ttttCTATTTTTTTTAAATATCTATATAGTTCTGTTAATGGTTTTTATTTTTAAAATTTTTTGTTTTTAGGCATAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGatat                                                                               <  2:410903/147‑1 (MQ=255)
                                                                             tttttttttttttAAATATTTATATAGGTCTGTTAATGTATTTTTTTTTTACATTTTTTTGTTTTTAGACATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCg                                                                             <  2:465992/141‑1 (MQ=255)
                                                                             tttttttttttttAAATATCTATATAGTTTTTTTAATGAATTTTATTTTTAAATTTTTTTGTTTTTAGACATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGatat                                                                               <  2:977536/139‑1 (MQ=255)
                                                                             ttCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCt                                                                           >  1:1319639/1‑151 (MQ=255)
                                                                                ttttttttttAAATATCTATATAGGTTTTTAAATGGTTTTTTTTTTTACATTTTTTTGTTTTTAGCCATATAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCttt                                                                         <  2:814277/145‑1 (MQ=255)
                                                                                tATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCtttt                                                                        >  1:609452/1‑151 (MQ=255)
                                                                                tATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAAATTAGCAATTAACCATTAACGCGTTGTTTCGCtttt                                                                        >  2:706039/1‑151 (MQ=255)
                                                                                     tttttAATTATTTATATAGTTTTGTTAATGGTTTTTATTTTTACAATTTTTTGTTTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAg                                                                    <  1:1195112/150‑1 (MQ=255)
                                                                                     tttttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTAAAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTaaa                                                                     >  1:1100004/1‑149 (MQ=255)
                                                                                          aaaTATTTATATGGGTCTGTAAATGTTTTTTATTTTTACAATTTTTTTTGTTTAGCCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTcgcg                                                               <  1:559170/150‑1 (MQ=255)
                                                                                          aaaTATCTATATAGGTTTGTTAATGGTTTTTATTTTTAAATTTTTTTGTGTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGt                                                              <  2:773536/151‑1 (MQ=255)
                                                                                          aaaTATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAAATTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTcgcg                                                               >  2:1238763/1‑150 (MQ=255)
                                                                                              atCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAAAGCTTGATATCGCTTTTAAAGTCGCGgtttt                                                          >  1:605757/1‑151 (MQ=255)
                                                                                               tCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTTAAGTCGCGttttt                                                          >  2:1387216/1‑150 (MQ=255)
                                                                                               tCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGtttt                                                           <  1:1396321/149‑1 (MQ=255)
                                                                                                cTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTTAAGTCGCGttttt                                                          >  2:273790/1‑149 (MQ=255)
                                                                                                   tataGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCatat                                                     >  1:306230/1‑151 (MQ=255)
                                                                                                   tataGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCata                                                      >  1:883332/1‑150 (MQ=255)
                                                                                                    ataGGTTTTTTAAGGGATTTTTTTTTTACATTTTTTTGTTTTTAGGCAAAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTc                                                         <  1:778246/146‑1 (MQ=255)
                                                                                                      aGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGAGATCGCTTTTAAAAGCGCGTTTTTTATATCCt                                                  >  2:457282/1‑151 (MQ=255)
                                                                                                       ggTCTGTTAATGGTTTTTTTTTTTAAATTTTTTTTTTTTTAGGCATAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTg                                                 <  2:1376415/151‑1 (MQ=255)
                                                                                                         ttttttAATGTTTTTTTTTTTTACATTTTTTTGTTTTTAGGAATAAAAAAACCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGt                                                <  2:951672/146‑1 (MQ=255)
                                                                                                         tCTGTTAATGGTTTTTATTTTTACAATTTTTTGTGTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGt                                                <  1:691301/150‑1 (MQ=255)
                                                                                                           tGTTAATGGATTTTTTTTTTACATTTTTTTGTTTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGtat                                              <  1:455371/150‑1 (MQ=255)
                                                                                                           tGTTAATGGATTTTATTTTTAAAATTTTTTGTTTTTAGGAATAAAAAAACCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGta                                               <  1:119026/149‑1 (MQ=255)
                                                                                                             ttaaggttttttttttttAAATTTTTTTGTTTTTAGAAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGtata                                             <  2:1080855/141‑1 (MQ=255)
                                                                                                             ttAATGTATTTTTTTTTTACATTTTTTTGTTTTTAGGAATAAAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGt                                                <  1:447235/146‑1 (MQ=255)
                                                                                                             ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGAAATCCACCATTAACGCTTGATATCGCTTTTTAAGTCGCGTTTTTCATATCCTGtat                                              >  2:978673/1‑148 (MQ=255)
                                                                                                             ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAAATTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGtata                                             >  2:1091381/1‑149 (MQ=255)
                                                                                                               aaTGGTTTTTTTTTTTAAATTTTTTTGTTTTTAGGAATAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACa                                           <  1:1217422/149‑1 (MQ=255)
                                                                                                               aaTGGTTTTTATTTTTAAATTTTTTTGTGTTTAGGCATAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTCAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAg                                          <  1:776647/150‑1 (MQ=255)
                                                                                                                aTGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAg                                          >  1:820018/1‑149 (MQ=255)
                                                                                                                aTGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAATCGCGTTTTTCATATCCTGTATACAGCt                                        >  2:826672/1‑151 (MQ=255)
                                                                                                                  ggATTTTATTTTTACAATTTTTTGTTTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCt                                        <  1:745008/149‑1 (MQ=255)
                                                                                                                   gATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAAACAACCATTAACGCTTGATATCGCTTTTAAAGACGCGGTTTTCATATCCCGTATACAGCTGa                                      >  1:1400714/1‑150 (MQ=255)
                                                                                                                     ttttttttttACATTTTTTTGTGTTAAGGCATAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCt                                        <  2:700271/146‑1 (MQ=255)
                                                                                                                     ttttttttttACAATTTTTTGTTTTTAGGAAAAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTg                                       <  1:420443/147‑1 (MQ=255)
                                                                                                                     ttttattttTACAATTTTTTGTGTTTAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGAcg                                    <  1:351986/150‑1 (MQ=255)
                                                                                                                      tttattttTAAAATTTTTTGTGTTTAGGCATAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGAcgc                                   <  1:68836/150‑1 (MQ=255)
                                                                                                                          tttttATAATTTTTTTTTTTTAGGAATAAAAAAAAAAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGAc                               <  2:173382/150‑1 (MQ=255)
                                                                                                                                    tttttGTGTTTAGGCATAAAAAAACCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCtt                     <  1:273797/150‑1 (MQ=255)
                                                                                                                                       ttGTTTTTAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCAt                  <  1:1371856/150‑1 (MQ=255)
                                                                                                                                        tgtgTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGGTAAAATATTTACAACTTAGCAATTAACCATTAAAGCTTGATATCGCTTTTAAAGGAGGGTTTTTCATATCCTGTATAAAGCTGATGCGGACGGGGAAACTTCATAc                >  2:1201270/1‑151 (MQ=255)
                                                                                                                                           gtTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCg              >  1:469932/1‑150 (MQ=255)
                                                                                                                                            tttAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCg              <  1:1183241/149‑1 (MQ=255)
                                                                                                                                               aGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCAc          >  1:583227/1‑150 (MQ=255)
                                                                                                                                                    tataAAAATCAACCCGCCATATGAACGGCGGGGTTAAATATTTACAAATTAGCAATTAACCATTAACCCTTGGTATCGCTTTTAAAGGCGCGTTTTTCATATCCTGTATTCAGCTGAAGCGGGCGGGGAAACTTCATACCGTc            >  1:1028066/1‑143 (MQ=255)
                                                                                                                                                     ataAAAATCAACCCGCCATATGAACGGCGGGGTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCAt   >  1:706094/1‑151 (MQ=255)
                                                                                                                                                      taAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGGCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCAtt  >  1:140030/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                      
CATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATT  >  NC_000913/752976‑753275

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