Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913753,1250TA32.5% 171.6 / 43.7 83intergenic (‑330/+60)ybgD/gltAputative fimbrial‑like adhesin protein/citrate synthase
Reads supporting (aligned to +/- strand):  ref base T (52/4);  new base A (0/27);  total (52/31)
Fisher's exact test for biased strand distribution p-value = 6.17e-18
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TATCATGAATTTGTTGT‑ATGATGAATAAAATATAGGGGGGTATAGATAGACGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTAC‑GTGCCATTGCGAAAATGACGGTGAACATGGAAGACGGAATACCCATCGCTTTCAG  >  NC_000913/752923‑753361
                                                                                                                                                                                                           |                                                                                                                                                                                                                                             
gtacaatgaTTTGTTGTAATAAAGAAAAAAATAGGGGGGGGTATAAATAGATTTTTTTTTAATAGGGTAATAAATGCTATTACAGATAATTTTTAAATTAAAAGTATGGGTTTGCTGATAATTTGTCTTTTTTTTTTTTTTTTAAATGTCTATATGGGTTTTAAATTGTATTTTATTTTAACAATTTTTTGTGTTAAGCAATAAAAAATCCAACCCGCCATATGAACGGCGGGTTaaa                                                                                                                                                                                                             <  2:622304/230‑1 (MQ=255)
   cATTAATTTTTTGT‑ATGATAAAAAAAATATGGGGGGTTATAGATAGACGTTTTTTTCATGGGTTTATAAATGCGACTACCATAAATTTTTTAATTGAAATTTTTGGTTTGCTAATAATTTGGGTTTTTCTTTTTTTTTTAAATATCTATATAGGTTTTTTAATGTTTTTTTTTTTAACAATTTTTTGTTTTTAGGAAAAAAAAAACCACCCCGCCATATGAACGGCGGGTTAAAATAtt                                                                                                                                                                                                        <  2:2019713/239‑1 (MQ=255)
          ttgttgT‑ATGATGAATAAAATATAGGGGGGTATAGATAGACGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAata                                                                                                                                                                                                          >  1:860993/1‑230 (MQ=255)
                        aaTAAAATATAGGGGGGTATAGATAGACGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCGTATAAAAATCAACCCGCCATATGAACGGCGGTTTAAAATATTTAAAACTTAGAAATCAACCa                                                                                                                                                                                    >  1:2252479/1‑239 (MQ=255)
                          tAAAATATAGGGGGGTATAGATAGACGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTTCAACTTAGCAATCAACCAtt                                                                                                                                                                                  >  1:2396989/1‑239 (MQ=255)
                                    gggggTATAGATAGACGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGat                                                                                                                                                                        >  1:182405/1‑239 (MQ=255)
                                        gTATAGATAGACGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAACAATCAACCATTAACGCTTTATATcc                                                                                                                                                                    >  1:1529476/1‑238 (MQ=255)
                                         tataGATATACGTTATTTTAAGAGTGTTAAAATTGGCACTATGTTGAAGTTTTTAATTGAAATTGTTGGTTTGTTGATAATTTGATTTTTTCTATTTTTTTAAAATATTTTTATGGTTTTTTTAAGGTTTTTTTTTTTTAAAATTTTTTGTTTTTAGGAATAAAAAAAACACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTa                                                                                                                                                                                 <  2:2407521/225‑1 (MQ=255)
                                         tataGATAGACGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTAAAACTTAGCAATCAACCATTAACGCTTGATATc                                                                                                                                                                     >  1:2358855/1‑237 (MQ=255)
                                                  aCGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTAAAACTTAGCAATCAACCATTAACGCTTTATATAGCTTTTTAAg                                                                                                                                                           >  1:2419586/1‑238 (MQ=255)
                                                   cGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTc                                                                                                                                                         >  1:1193598/1‑239 (MQ=255)
                                                   cGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAAAGCTTGcta                                                                                                                                                                       >  1:2407521/1‑225 (MQ=255)
                                                      taTTTTCATAGGTTTAAAAATCGTACTACAATAAATTTTTTAATTAAAATTTTTGGTTTGTTATTAATTTGATTTTTTTTTTTTTTTTTAAATTTTTTTATAGTTTTGTTAATGGTTTTTTTTTTTAAAATTTTTTGTTTTTAAACATAAAAAAAACAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTcgcg                                                                                                                                                      <  2:1940889/238‑1 (MQ=255)
                                                        ttttCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTaa                                                                                                                                                                                              >  1:1786077/1‑196 (MQ=255)
                                                                   ttATAAATGCAACTAACATGAAGTTTTTAATTAAAATTATTGGGTTGTTGATTTTTTGATTTTTTTTTTTTTTTTTAAATATTTATATAGGTTTGTTAATGTATTTTTTTTTTAAAATTTTTTGTGTTTAGGAAAAAAAAAAACACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCaa                                                                                                                                                                                           <  2:389885/189‑1 (MQ=255)
                                                                           gCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTTAAATATTTTAAACTTATCAATCAACACTTAAAGCTTTATATAGCTTTTTAAATCACGTTTTTTATATcccttat                                                                                                                                     >  2:2090283/1‑230 (MQ=255)
                                                                             gACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTa                                                                                                                                                                                 >  1:389885/1‑189 (MQ=255)
                                                                                  ccATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTATCAATCAAACATTAACGCTTGATATTGCTTTTAAAATCGCGTTTTTCCTATTCTTTTTACAATTGAcgc                                                                                                                          >  1:1162624/1‑239 (MQ=255)
                                                                                          tttttAATTGAAAGTATTGGTTTGCTGATAATTTGAGCTTTTCTTTTCTTTTTAAATATCTATATAGGTCTTTTAATGGTTTTTTTTTTTACAATTTTTTGTGTTTAGGCATAAAAAAACCAACCCGCCATATGAACGGCGGGTTaaa                                                                                                                                                                                                             <  1:274702/148‑1 (MQ=255)
                                                                                          tttttAATTAAAATTATTGGGTGGTTAATAATTTGATTTTTTCTTTTATTTTTAATTATTTATATAGGTCTTTAAATGGATTTTATTTTTTCATTTTTTTTTTTTAAGGAATAAAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGc                                                                                                                  <  2:1162624/239‑1 (MQ=255)
                                                                                            tttAATTAAAATTATGGGTGTGTTGAGATTTGGAGTGTTTTTTTTCTTTTTAATTATCTAAATGGGTCTTAGAATGGATTTTTTTTTAACAATTTTTGTTTTTTAGGCATAAAAAAACAAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGGCGCGGACGGGCaa                                                                                                                <  2:1193598/239‑1 (MQ=255)
                                                                                                    aaaGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTTAAATATTTTCaa                                                                                                                                                                                                   >  2:274702/1‑148 (MQ=255)
                                                                                                    aaaGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTTAAGTCGCGTTTTTTATATACTGTATACAGCTGACGCGGATGGGCAATCTTCAt                                                                                                         >  1:388985/1‑238 (MQ=255)
                                                                                                     aaGTATTGAGTTGCTGATATTTTGATTTTTTTTTTTTTTTTAAAATATTTATATGGTTTTGTAAATGTATTTTTTTTTTACATTTTTTGGTTTTTAGAAATAAAAAAAACAACCCGCCATATGAACGGCGGGTTAGAATATTTACAACTTAGCAATCAACCATTAACGCTTGatat                                                                                                                                                                      <  2:102010/176‑1 (MQ=255)
                                                                                                        ttttGGGTTGCTGATAATTTGAGCTGTTTTTTTTTTTTTAAATTTCTATATAGGTCTGTTAATGGATTTTTTTTTTACAATTTTTTGTGTTTAGGCATAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTc                                                                                                                                                <  1:585179/193‑1 (MQ=255)
                                                                                                               ttGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAGAATATTTACAACTTAGCAATCAACCATTAAAGCTTGATATCGCTTTTaaa                                                                                                                                                            >  1:102010/1‑176 (MQ=255)
                                                                                                                 gCTGATAATTTGGGCTTTTTTTTTTTTTTTAAATATCTATATAGGTTTGTTAATGAATTTTATTTTTACAATTTTTTGTGTTTAGGAATAAAAAAATAAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCa                                                                                           <  1:1317052/239‑1 (MQ=255)
                                                                                                                  cTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTTCAACTTAACAATCAACCATTAACGCTTGATATTGCTTTTTAAGTCCCGTTTTTTATTTCCTGta                                                                                                                                      >  2:585179/1‑195 (MQ=255)
                                                                                                                  cTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTTAAGTCGCGGTTTTCATATCCTGTATACATCTGACGCGGACGGGCAATCTTCATACCGTCa                                                                                                  >  1:1501618/1‑231 (MQ=255)
                                                                                                                                  ttCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCg                                                                          >  1:1713763/1‑239 (MQ=255)
                                                                                                                                            tttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTTTATACAGCTGACGCGGCCGGGCAATCTTCATACCTTCACTTTGCa                                                                                           >  1:2689455/1‑212 (MQ=255)
                                                                                                                                            tttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCa                                                                 >  1:85690/1‑238 (MQ=255)
                                                                                                                                            ataagagaaaaatataGGTCTGTTAATGGTTTTTATTTTTACAATTTTTTGTGTTTAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCaa                                                                <  1:2144741/228‑1 (MQ=255)
                                                                                                                                                   atttatataGGTTTGTAAATGATTTTTATTTTAAAATTTTTTTGTTTTAAGCCATAAAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGAcg                                                                                                                           <  1:36613/170‑1 (MQ=255)
                                                                                                                                                   atttagataGTTCTTTTAATGTTTTTTTTTTTTAATTTTTTTTTTGTTTAAGCATAAAAAACACACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTAC‑g                                                        <  2:85690/233‑1 (MQ=255)
                                                                                                                                                   atCTATATAGTTCTGTTATGGGATTTTATTTTTACATTTTTTTTTGTTAAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTAC‑g                                                        <  1:989476/239‑1 (MQ=255)
                                                                                                                                                      tatatgGGTTTTTTAATGTTTTTTTTTTTTACATTTTTTTGTTTTTAGGAAAAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGc                                                                                                                  <  2:2168394/179‑1 (MQ=255)
                                                                                                                                                      tatataGTTTTTTTAATGGTTTTTTTTTTTACATTTTTTTGTTTTTAGGCATATAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                                                                                                                                              <  1:393959/151‑1 (MQ=255)
                                                                                                                                                      tatataGGTCTTTAAATGGATTTTTTTTTTACATTTTTTTGTGTTTAAGTATAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                                                                                                                                              <  1:2211588/151‑1 (MQ=255)
                                                                                                                                                      tatataGGTCTATTAATGGTTTTTTTTTTTACAATTTTTTTTTTTTAGACATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCtt                                                                                                                                                                 <  2:797667/132‑1 (MQ=255)
                                                                                                                                                      aatataGTTTTATTAATGGTTTTTATTTTTACATTTTTTTTTTTTAATGCAAAAAAAAATAACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAATTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCa                                                                                                  <  2:928396/194‑1 (MQ=255)
                                                                                                                                                      aatataGGGCGGTTAATGAATTTTATTTTTAAAATTTTTTGTGTTGAGAAAAAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                                                                                                                                              <  1:2558574/150‑1 (MQ=255)
                                                                                                                                                        tttaGGTCTGTTAATGTTTTTTATTTTTACATTTTTTTTGTTTTAGGCATAAAAAATCCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCa                                                                                                                                                                                        <  1:1679721/105‑1 (MQ=255)
                                                                                                                                                        tataGGTTTTTTAATGTTTTTTTTTTTTACAATTTTTTGTTTTTAGGCATAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGAcgcg                                                                                                                         <  1:2568214/170‑1 (MQ=255)
                                                                                                                                                        tataGGTTTGTTAAGGGTTTTTATTTTTAAATTTTTTTGTTTTTAGACATAAAAAAACCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCAt                                                                                                                                                                                   <  2:1275477/112‑1 (MQ=255)
                                                                                                                                                             ttCTTTTAATTGTTTTTATTTTTACATTTTTTTGTGTTAAGGAATAAAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGc                                                                                                                  <  2:1162457/171‑1 (MQ=255)
                                                                                                                                                             gTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTTAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATTGCTTTTAAAGTCGCGTTTTTTATATCCTTTATACATCTGACGCGGACGGGCa                                                                                                                 >  2:36613/1‑173 (MQ=255)
                                                                                                                                                             gTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCCACGGTAC‑GTGCCATTGCg                                              >  1:161326/1‑239 (MQ=255)
                                                                                                                                                             gTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGGTAAAATATTTACAACTTAGCAATCAACCATTAAAGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATTCTGTTTACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTTCCTCCAGTGGGCGATCCAGCCAACGGTAC‑ATGCCATTGCg                                              >  1:296702/1‑239 (MQ=255)
                                                                                                                                                              tCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTTAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCTTCCCGTCACTGTGCATTTCTCCCCACTGGGCGGTCCAACCCACCGTAC‑GTTCCCTTGCGa                                             >  2:209714/1‑239 (MQ=255)
                                                                                                                                                                tGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTTAAATTTTTAAAACTTAGCAATCAAACATTAACGCTTGATATCGCTTTTTAAATCGCGTTTTTCATATCCTGAATACAGCTGACGCGGAAGGGCAATCTTCATACCGTCACTGTGTAttt                                                                                        >  1:928396/1‑195 (MQ=255)
                                                                                                                                                                tGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTTAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCTCTCCCGTGGGCGATCCAGCCAACGGTAC‑TGGCCATTGCGaaa                                           >  2:241727/1‑239 (MQ=255)
                                                                                                                                                                tGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTTAAATATTTACAACTTAGCAAACAACCATTAACGCTTGATATCCCTTTTAAAGTCGCGTTTTTTATATCCTGTATAAAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTTTCTCCAGTGGGCGATCCAGCCAACGGTAC‑GTGCCATTgtgcaa                                           >  2:2699549/1‑234 (MQ=255)
                                                                                                                                                                tGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTTAAATATTTACAAATTAGCAATCAACCATTAAAGCTTGGTATCGCTTTTAAAAACGCGTTTTTTATATCCTGtata                                                                                                                                    >  2:393959/1‑151 (MQ=255)
                                                                                                                                                                tGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTcgc                                                                                                                                                       >  1:797667/1‑132 (MQ=255)
                                                                                                                                                                tGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTAC‑GTGCCATTGCGaaa                                           >  1:615886/1‑239 (MQ=255)
                                                                                                                                                                tGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAACCAACGGTAC‑GTGCCATTGCGaaa                                           >  1:338854/1‑239 (MQ=255)
                                                                                                                                                                tGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAAACTTCATa                                                                                                        >  1:2168394/1‑179 (MQ=255)
                                                                                                                                                                tGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAAACATTAACGCTTTATATCGCTTTTAAAGTCGCGTTTTTCCTATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCCCttgg                                                                                             >  1:1331300/1‑187 (MQ=255)
                                                                                                                                                                tGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAAATTAGCAATTAACCATTAACGCTTGATATCGCTTTTTAAGTGGCGGTTTTCTTAACCTTtata                                                                                                                                    >  2:2211588/1‑151 (MQ=255)
                                                                                                                                                                tGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGGTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTTAAATCGCGTTTTTTATTACCTGTATACCGCTGACGCGGACGGGCAATCTTCATACCTTCACTTTGCATTTCGCTCCATTGGGCGATCCAACCAACGGTACTGTACC‑TTGCGaaa                                           >  2:2045965/1‑239 (MQ=255)
                                                                                                                                                                tGTTAATGGATTTTATTTTTACAATTTTTTGGGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGGTTAAAAATTTACAAATTAACAATCAACCATTAACGCTTGATATAGCTTTTATAGTCGCGTTTTTCTTATCCTTTATAAACCTGAAGTGGACGGGGAATCTTCATACCGTCACTGTTCATTTCGCCTCAGCGGGCGACACCGCCACCGGGAC‑GTGCCATTTCGgaa                                           >  2:966225/1‑239 (MQ=255)
                                                                                                                                                                  ttAATGGTTTTTTTTTTTAAAATTTTTTGTTTTTAGAAAAAAAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCgg                                                                                                                        <  2:1804176/161‑1 (MQ=255)
                                                                                                                                                                  ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAat                                                                                                                                                                                                           >  1:2118040/1‑78 (MQ=255)
                                                                                                                                                                  ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTTAAACTTAGAAAACAACCATTAACGCTTGATATCCCTTTTTAAATCGCGTTTTTTATATCCTTTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCCCTGTGCATTTTGCTCCCGTGGGCGATCCAGCCAACGGGAC‑GTGCCATTGCTATAAt                                         >  2:1071840/1‑239 (MQ=255)
                                                                                                                                                                  ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTATAACTTAGCAATCAACCATTAACGCTTTATATTGCTTTTAAAGTCGCGTTTTTTATATCCTGTATACAGCTTACGCGGACGGGCAAt                                                                                                               >  2:2568214/1‑170 (MQ=255)
                                                                                                                                                                  ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGa                                                                                                                                                                         >  1:1275477/1‑112 (MQ=255)
                                                                                                                                                                  ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTTTGGCTTTCTCTCCCGTGGGCGGTCCAGCCAACCGGAC‑Gttc                                                     >  2:2144741/1‑225 (MQ=255)
                                                                                                                                                                  ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACg                                                              >  1:1475253/1‑219 (MQ=255)
                                                                                                                                                                  ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAAAGCTTGATATCGCTTTTTAAATCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTAGCTCCCGTTGGCGATCAAGCCCACGGTAC‑GTGcccatgtgcaaaa                                         >  2:1137821/1‑228 (MQ=255)
                                                                                                                                                                  ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGGTTTAAAATTTACAACTTTGCAATCAACCATTAACGCTTGATATCGCTTTTTAAATCGCGTTTTTTATATccctttta                                                                                                                                    >  2:2249528/1‑143 (MQ=255)
                                                                                                                                                                   tAATGTATTTTTTTTTTAAATTTTTTTGTGTTTAGGAATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                                                                                                                                              <  1:797628/138‑1 (MQ=255)
                                                                                                                                                                      tGGATTTTTTTTTTACATTTTTTTGTGTTTAGGCATATAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCt                                                                                    <  1:636912/193‑1 (MQ=255)
                                                                                                                                                                       ggATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTTAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATa                                                                                                        >  1:1162457/1‑172 (MQ=255)
                                                                                                                                                                        ttttttattttTACATTTTTTTGTTTTTAGGCATAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGc                                                                                                                  <  2:492248/159‑1 (MQ=255)
                                                                                                                                                                            ttattttTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTTAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATc                                                                                                              >  1:1804176/1‑161 (MQ=255)
                                                                                                                                                                            ttattttTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTTATACCGTCACTGGGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTAC‑GTGc                                                     >  1:747896/1‑217 (MQ=255)
                                                                                                                                                                            ttattttTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTAC‑GTGCCATTGCGAAAATGACGGTGAAc                               >  1:873053/1‑239 (MQ=255)
                                                                                                                                                                                ttttACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTTAAATATTTTCAAATTAACAATCAACCATTAAAGCTTGATATTGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAACTGAAGCGGAAGGGGAAACTTCATACTGTCACTTTGCATTTTGCTCCAGGGGGCg                                                                          >  2:636912/1‑193 (MQ=255)
                                                                                                                                                                                 tttAAAATTTTTGTTGTTTAGGAAAAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTAC‑GTGCCATTGCGAAAATGACGGTGAACATGGa                          <  2:182405/239‑1 (MQ=255)
                                                                                                                                                                                  ttACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTTAAAAATTTACAACTTAGCAAACAACCATTAAACCTTGATATTGCTTTTTAAGTTGCGTTTTTCATATCCTTTATAAAGCTGACGCGGACGGGAAAACTTCATACCGTCACTGTGCATTTCTCTCCAGCGGGCGAGCCCACCCACGGCAC‑GTGCCATTGCGAAAATAACGGGGAACATggga                         >  2:1514062/1‑237 (MQ=255)
                                                                                                                                                                                  ttACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATa                                                                                                        >  1:492248/1‑161 (MQ=255)
                                                                                                                                                                                               tttttAGAAATAAAAAAACCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTgg                                                                             <  2:134515/173‑1 (MQ=255)
                                                                                                                                                                                                        aaataAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTAC‑GTGCCATTGCGAAAATGACGGTGAACATGGAAGACGGAATACCCATCGCttt     <  1:1762389/235‑1 (MQ=255)
                                                                                                                                                                                                         tataAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGatat                                                                                                                                                                      <  2:1201641/76‑1 (MQ=255)
                                                                                                                                                                                                         tataAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGatat                                                                                                                                                                      >  1:1201641/1‑76 (MQ=255)
                                                                                                                                                                                                         tataAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGc                                                                   >  1:134515/1‑175 (MQ=255)
                                                                                                                                                                                                         tataAAAATCAACCCGCCATATGAACGGCGGGGTAAAATATTTACAAATTAACAACCAACCCTTAAAGCTTTGTATCGCATTTAAAGTCGGGTTTTACCTATTCTGTATAAACATGAAGAGGACGGGGCATATTGACACCCGCACCGTGCATTTCGCTCCCGTGGGGGAGCCAGCCAACGGGAA‑GTGCCACTGCGTACATAACCAGTAACCTGCAATACGGAAAACCCACCGCCTTCaa  >  2:60226/1‑238 (MQ=255)
                                                                                                                                                                                                           taAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTAC‑GTGCCATTGCGAAAATGACGGTGAACATGGAAGACGGAATAc               >  1:459441/1‑224 (MQ=255)
                                                                                                                                                                                                           taAAAATCAACCCGCCATATGAACGGCGGGGTaaaa                                                                                                                                                                                                            >  1:933208/1‑36 (MQ=255)
                                                                                                                                                                                                           |                                                                                                                                                                                                                                             
TATCATGAATTTGTTGT‑ATGATGAATAAAATATAGGGGGGTATAGATAGACGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTAC‑GTGCCATTGCGAAAATGACGGTGAACATGGAAGACGGAATACCCATCGCTTTCAG  >  NC_000913/752923‑753361

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑