Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,383,718 C→A 19.6% E82* (GAA→TAA)  ycjW ← LacI family putative transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,383,7180CA19.6% 133.1 / 5.3 73E82* (GAA→TAA) ycjWLacI family putative transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base C (49/7);  new base A (2/13);  total (52/21)
Fisher's exact test for biased strand distribution p-value = 1.10e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCATATTGACCTTCAACTTTGCCAATCACCACCACCGGAATATCATATTTATCGAGTTGGGCAAAAAATGACTCATCCGCTGGCGAACTGAGCATAATAATGCCTTTAATCATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGATTGTGCGACGTTTGTAATATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGAGATATTAGTTTGCTTATTCAGCACGCGTGATACGGTGGATTTTGATACGCCTGCAACCCTGGCAATATCATAAA  >  NC_000913/1383488‑1383948
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                            
ggCATATTGACCTTCAACTTTGCCAATCACCACCACCGGAATATCATATTTATCGAGTTGGTCAAAAAATGACTCATCCGCTGGCGAACTGAGCATAATAATGCCTTTAATCATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCTGGATTGTGCGCCGTTTGTAATATCACGTCACAACTTTCTTCTTCCGCTTTGGCGGTGTTTGCATCCACAACTTccatgaaaa                                                                                                                                                                                                                                      >  1:1071400/1‑231 (MQ=255)
ggCATATTGACCTTCAACTTTGCCAATCACCACCACCGGAATATCATATTTATCGAGTTGGGCAAAAAATGACTCATCCGCTGGCGAACTGAGCATAATAATGCCTTTAATCATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGATTGTGCGCCGTTTGTAATATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTGATGGCATGAAAAAATTCaaaaaaaa                                                                                                                                                                                                                                      >  1:199012/1‑239 (MQ=255)
ggCATATTGACCTTCAACTTTGCCAATCACCACCACCGGAATATCATATTTATCGAGTTGGGCAAAAAATGACTCATCCGCTGGCGAACTGAGCATAATAATGCCTTTAATCATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGATTGTGCGACGTTTGTACTATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTGTTGGCATGTAAAACTTACGAGaaaa                                                                                                                                                                                                                                      >  2:1007184/1‑239 (MQ=255)
ggCATATTGACCTTCAACTTTGCCAATCACCACCACCGGAATATCATATTTATCGAGTTGGGCAAAAAATGACTCATCCGCTGGCGAACTGAGCATAATAATGCCTTTAATCATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGATTGTGCGACGTTTGTAATATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTGATGGCATGCAAAACTTCAGAGaaaa                                                                                                                                                                                                                                      >  2:676171/1‑239 (MQ=255)
ggCATATTGACCTTCAACTTTGCCAATCACCACCACCGGAATATCATATTTATCGAGTTGGGCAAAAAATGACTCATCCGCTGGCGAACTGAGCATAATAATGCCTTTAATCATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGATTGTGCGACGTTTGTAATATAACGTCGAAACCTTCTTCTTCACCTTTGGCGGGGATGGCATGCAAAACTTTAGAGaaaa                                                                                                                                                                                                                                      >  1:851970/1‑239 (MQ=255)
    tatTGACCTTCAACTTTGCCAATCACCACCACCGGAATATCATATTTATCGAGTTGTGCAAAAAATGACTCATCCGCTGGCGAACTGAGCATAATAATGCCTTTAATCATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGATTGTGCGACGTTTGTAATATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTGATGGCATGCAAAAATTCTAGGAAAACCgg                                                                                                                                                                                                                                 >  2:405727/1‑239 (MQ=255)
           cTTCAACTTTGCCAATCACCACCACCGGAATATCATATTTATCGAGTTGGGCAAAAAATGACTCATCCGCTGGCGAACTGAGCATAATAATGCCTTTAATCATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGATTGTCCGACGTTTGTAATATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTGAGTCTATGCAAAACTTCAACAACAAACGGATAAcc                                                                                                                                                                                                                           >  1:902141/1‑238 (MQ=255)
                     gCCAATCACCACCACCGGAATATCATATTTATCGAGTTGGGCAAAAAATGACTCATCCGCTGGCGAACTGAGCATAATAATGCCTTTAATCATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGATTGTGCGACGTTTGTAATATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGttt                                                                                                                                                                                                               >  2:189979/1‑239 (MQ=255)
                                ccaccGGAATATCATATTTATCGAGTTGGGCAAAAAATGACTCATCCGCTGGCGAACTGAGCATAATAATGCCTTTAATCATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGATTGTGCGACGTTTGTAATATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACgg                                                                                                                                                                                                    >  2:1219566/1‑239 (MQ=255)
                                   ccGGAATATCATATTTATCGGGTTGGGAAAAAAATGATCACTCCGTGGGGCAACTGAGCAAAATAATGCCTTTAATCATTTTCTGCTTATTTTTGTTTTCGCTTTTTTGCAGTTCTTCTGCCGGATTGTGCGACGTTTGTAATATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTGATGGCATGCAAAACTTAAGAAAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTa                                                                                                                                                                                                 <  2:1057574/239‑1 (MQ=255)
                                           tCATATTTATCGAGTTGGGCAAAAAATGACTCATCCGCTGGCGAACTGAGCATAATAATGCCTTTAATCATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGATTGTGCGACGTTTGTAATATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCa                                                                                                                                                                                          >  2:1108545/1‑238 (MQ=255)
                                              tatTTATCGAGTTGGGCAAAAAATGACTCATCCGCTGGCGAACTGAGCATAATAATGCTTTTAATCATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGATTGTCCGACGTTTGTAATATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTGTTGGCATGCTAAACTTCAGAGAAAAACGGGTTAACCACCGTAGTTTTggggtaaaggggtgaaaacaacaa                                                                                                                                                                                      >  1:928876/1‑217 (MQ=255)
                                                 ttATCGAGTTGGGCAAAAAATGACTCATCCGCTGGCGAACTGAGCATAATAATGCCTTTAATCATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGATTGTGCGACGTTTGTAATATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATAACCATAAt                                                                                                                                                                                   >  2:552028/1‑239 (MQ=255)
                                                   aTCGAGTTGGGCAAAAAATGACTCATCCGCTGGCGAACTGATCATAATAATGCCTTTAATCATTTTCTGCTTAATTGTGCTTTCGCATTTTTGCAGATCTTCTGCCGGATTGTTCGACGTTTGTAATATCCCGTCGTAACCTTCTTCTTCAGCGTTGGAGGGGGTGGCATTCAAAACTTAAGGGCAAAACAGATAACCCGCCGTAGTGTTGGCGGAACAAGAAGAAAtacccaaaaggg                                                                                                                                                                                 >  1:1065707/1‑228 (MQ=255)
                                                   aTCGAGTTGGGCAAAAAATGACTCATCCGCTGGCGAACTGAGCATAATAATGCCTTTAATCATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGATTGTGCGACGTTTGTAATATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTGATGGGATGCAAACCTTCAGGGAAAAACGGATTCCCCCCCTGAATTTTGGTCGAACGGGtgaaaaaacacaaaaagg                                                                                                                                                                                 >  1:738582/1‑222 (MQ=255)
                                                     cGAGTTGGGCAAAAAATGACTCATCCGCTGGCGAACTGAGCATAATAATGCCTTTAATCATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGATTGTGCGACGTTTGTAATATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCa                                                                                                                                                                              >  2:167846/1‑239 (MQ=255)
                                                     cGAGTTGGGCAAAAAATGACTCATCCGCTGGCGAACTGAGCATAATAATGCCTTTAATCATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGATTGTGCGACGTTTGTAATATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTACTTTTGTCCGAACGGGAAAAAAAAACCAAAATGGCa                                                                                                                                                                              >  1:93901/1‑239 (MQ=255)
                                                     cGAGTTGGGCAAAAAATGACTCATCCGCTGGCGAACTGAGCATAATAATGCCTTTAATCATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGATTGTGCGACGTTTGTAATATCACGTCGAAACCTTCATCTTCAGCTTTGTCGGTGAGGGCGTGCAAAACTTCAGAGAAAAACGGATTACCCGCGGGAGTTTTGTTCGAACGGGTAGAAAtaaaaaaaaaggaa                                                                                                                                                                               >  1:294173/1‑226 (MQ=255)
                                                      gAGTTGGGCAAAAAATGACTCATCCGCTGGCGAACTGAGCATAATAATGCCTTTAATCATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGCTTGTGCGACGTTTGTACTATCACGTCGAAACCTTCTTCTTCCGCTTTGGCGGCGATGCTATGGAAAACTTCCGAGACAAACGGATTACCCACCGGTGTTTTTGGTCAACGGGGAGAAATCACCaaaaagagct                                                                                                                                                                              >  1:261425/1‑233 (MQ=255)
                                                         ttGGGCAAAAAATGACTCATCCGCTGGCGAACTGAGCATAATAATGCCTTTAATCATTTTCTGCTTAATTTTGCTTTCGCATTTTTTCAGATCTTCTGCCGGATTGTGCGACGTTTGTAATATCAAGTCGAAACCTTCTTCTTCAGCTTTGGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGATTACCCCCCGTAGTTTTGGGCGAAAGGGTAGaaaaccacaacaaagcaatca                                                                                                                                                                           >  1:843781/1‑221 (MQ=255)
                                                         ttGGGCAAAAAATGACTCATCCGCTGGCGAACTGAGCATAATAATGCCTTTAATCATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGTTTGTGCGACGTTTGTAATATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGAtt                                                                                                                                                                                                                              <  2:342781/189‑1 (MQ=255)
                                                         ttGGGCAAAAAATGACTCATCCGCTGGCGAACTGAGCATAATAATGCCTTTAATCATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGATTGTGCGACGTTTGTAATATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTGATGGCATGCCAAACTTCAAAGAAAAATGGATTACCCCCCGTTGTTTTGGGCCAAACGGTAGAAAACACCACaaaaccaaaca                                                                                                                                                                           >  1:544640/1‑230 (MQ=255)
                                                         ttGGGCAAAAAATGACTCATCCGCTGGCGAACTGAGCATAATAATGCCTTTAATCATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGATTGTGCGACGTTTGTAATATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTGATGGCATGCAAAACTTCAGGGAAAAACGGATTTCCCGCCGTAGTTTTGGTCGAAAGGGTAGAAAAAACCAAAATggaaaaaa                                                                                                                                                                          >  1:128759/1‑233 (MQ=255)
                                                         ttGGGCAAAAAATGACTCATCCGCTGGCGAACTGAGCATAATAATGCCTTTAATCATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGATTGTGCGACGTTTGTAATATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGAtt                                                                                                                                                                                                                              >  1:342781/1‑189 (MQ=255)
                                                         ttGGGCAAAAAATGACTCATCCGCTGGCGAACTGAGCATAATAATGCCTTTAATCATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGATTGTGCGACGTTTGTAATATCACGTAGAAACCTTCTTATTCGGCTTTGGCGGGGCTGGCTGACAAAACTTCATAGAAAAACGGATTACCCACCGCAGTTTTGGTCGAACGGTGAGaaaaaaacaaaaaggcaaaaa                                                                                                                                                                           >  1:1314287/1‑221 (MQ=255)
                                                           gggCAAAAAAGGCCTCATCGGCGGGCAAACTGAAAATAATAATGCCTTAAATCTTTTTCTTTTTAATTTTGTTTTCGCTTTTTTGCAGTTCTTCTGCCGGATTGTGGGACGTTTGAAATACCACGGCGAAACTTTCTTCTCCAGTTTTGGGGGTGATGGCAGGAAAAACTCAAGAAAAAAACGGATTCCCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCaaa                                                                                                                                                                         <  2:851970/238‑1 (MQ=255)
                                                                      gACTCATCCGCTGGCGAACTGAGCATAATAATGCCTTTAATCATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGATTGTGCGACGTTTGTAATATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAAAggc                                                                                                                                                             >  2:1035964/1‑238 (MQ=255)
                                                                                       aaCTGAGCATAATAATGCCTTTAATCATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGATTGTGCGACGTTTGTAATATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGATTTCCCGCCTTAGTTTTGGTCGACAGGGTAGAAATCACCAAAATGGCACAAAAAACCGAAAAGGGTAGGGCGCGGGCCAga                                                                                                                                           >  1:1003052/1‑238 (MQ=255)
                                                                                       aaCTGAGCATAATAATGCCTTTAATCATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGATTGTGCGACGTTTGTAATATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCa                                                                                                                                             >  2:115869/1‑237 (MQ=255)
                                                                                           gAGTATAATAATGCCTTTAATTATTTTCTTTTTAATTTTGTTTTCCCTTTTTTGCAGATCTTCTGCCGGATTGGGCGACGTTTGAAATATCACGTCAAACCTTTCTTCTCCGGTTTTGGGGGTGATGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGt                                                                                                                                       <  2:635890/239‑1 (MQ=255)
                                                                                                   taatGCTTTTAATCTTTTTCTGTTTATTTTTGCTTTCCTATTTTTGCAGATCTTCTGCGGGATTGTGCGACGTTTGAAATATCAGTGCGAACCTTTCTTCTTCAGTTTTGGGGGTGAGGGCTGGCAAAACTTCAAAAAAAAACGGTTTCCCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGa                                                                                                                               <  2:199012/239‑1 (MQ=255)
                                                                                                   taatCCCTTTAATCATTTTCTTTTTATTTTTGCTTTCGCTTTTTTGCAGTTCTTCTGCCGGATTGTGCGACTTTTTAAATTCCACGTCGAACCCTTCTTCTTCAGTTTTGGCGGTGATGGCATGAAAAACTTCAGAAAAAACCGTATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGa                                                                                                                               <  2:128759/239‑1 (MQ=255)
                                                                                                      tGCCTTTAATCATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGATTGTGCGACGTTTGTAATATCCCGTCGAAACCTTCTTCTTCAGCTTTGGCGGTGATGGCATGCAAAACTTCCAAGACAAACGGATTTCCCGACGTGGTTTTGGTCGAACCGGCAGAAATACCCCTAATGGCACAAAAACCCGAAGGGTGAGCGCAGCGGGTCAGCtgttatggatgaaca                                                                                                                            >  1:1183959/1‑225 (MQ=255)
                                                                                                      tGCCTTTAATCATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGATTGTGCGACGTTTGTAATATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTGGTGGCATTCAAAACTTCAGGGAAAAACGGATTACCAGCCGTAGTTTTGGCTGCACCGGTAGAAAACATCCACATGGCATTCAAACCCGAAGGGGTACACGCGGGGGCCacacttttgtgtgaaaac                                                                                                                            >  1:50231/1‑222 (MQ=255)
                                                                                                      tGCCTTTAATCATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGATTGTGCGACGTTTGTAATATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGAt                                                                                                                              >  2:1005064/1‑237 (MQ=255)
                                                                                                        ccTTTAATCATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGATTGTGCGACGTTTGTATTATAACTTCGAAACCTTCTTCTTCAGCTTTGGCGGTGTTGGCATGCTAAACTTCAGCAGAAATCGAATTACCCGCCGTAGCTTTGTTCGACTTGTTAGTAATCACCCTAATAGCATCACAATCCGAAGAGGTCATCACGCGAGCCACCCTGCTTGGCTggatcgt                                                                                                                          >  2:91143/1‑233 (MQ=255)
                                                                                                        ccTTTAATCATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGATTGTGCGACGTTTGTAATATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTCGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGAt                                                                                                                              >  2:150392/1‑235 (MQ=255)
                                                                                                        ccTTTAATCATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGATTGTGCGACGTTTGTAATATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGCAGAGGGCAGCCCGCGGGCCCGCTTGTTTGGCTGAAACTg                                                                                                                          >  2:149865/1‑239 (MQ=255)
                                                                                                        ccTTTAATCATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGATTGTGCGACGTTTGTAATATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATAACCATAATTGCATCAAAACCCGAAGAGGTCACCCGCGGGGCCAGCTTGTTTGGGTGGTaaag                                                                                                                          >  1:1291302/1‑236 (MQ=255)
                                                                                                        ccTTTAATCATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGATTGTGCGACGTTTGTAATATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGGAGTTTTGGTCGAAGGGGGAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGACCGGGCCAGCCTGTTTGGccgcaactg                                                                                                                          >  2:1286818/1‑231 (MQ=255)
                                                                                                           ttAATCATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGATTGTGCGACGTTTGTAATATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTaa                                                                                                                       >  2:489449/1‑239 (MQ=255)
                                                                                                           ttAATCATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGATTGTGCGACGTTTGTAATATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTaa                                                                                                                       >  2:1261537/1‑239 (MQ=255)
                                                                                                           ttAATCATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGATTGTGCGACGTTTGTAATATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCAGAATGGCATCAAAACCCGAAGAGGTCAGCGCCCGGGCAAGCTTGTTTGGCTGATACTGTaa                                                                                                                       >  2:765702/1‑239 (MQ=255)
                                                                                                           ttAATCATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGATTGTGCGACGTTTGTAATATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGGTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGCCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTaa                                                                                                                       >  2:1064123/1‑239 (MQ=255)
                                                                                                           taaattgttttCTGTTTAATTTTTGTTTCGTATTTGTTCAGTTCTTCGGCGGGTTTGTGGGACTTTTGAAATATCCCGTCAAAACCTTCTTCTCCAGTTTTGGGGGTGATGGCATGCAAAACTCAAGAAAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTaa                                                                                                                       <  2:693657/232‑1 (MQ=255)
                                                                                                             aaTCATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGATTGTGCGACGTTTGTAATATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCTTCAAAACCCGAAAAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGGAAtt                                                                                                                     >  2:1070691/1‑239 (MQ=255)
                                                                                                                                   cTTTCGCATTTTTGCAGATCTTCTGCCGGATTGTGCGACGTTTGTAATATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGCAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAAAACttt                                                                                               >  2:1174247/1‑239 (MQ=255)
                                                                                                                                      tCGCATTTTTGCAGATCTTCTGCCGGATTGTGCGACGTTTGTAATATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTTATTGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGATCGGGTAGAAATCACCATAATGGCATCAAAAACCGAAGAGGTCAGCGCGCGGGCCAGCCTGTTTTGCTGATACTGTAATTCTTCAATGGGCCGGAAACCTTTTTc                                                                                            >  1:811084/1‑239 (MQ=255)
                                                                                                                                       cGCATTTTTGCAGATCTTCTGCCGGATTGTGCGACGTTTGTAATATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGTCTGATACTGTAATTCTTAAATGGCCCGGAAAACTTTTTc                                                                                            >  1:928133/1‑238 (MQ=255)
                                                                                                                                           ttgttGCAGATCTTGTGGCGGTGTGTGGCACGTTTGAAAATTCCCGTCAAACCCTTCTTCTGCAGTTTTGGCGTGGAGGGAATGCAAAACTTAAGAAAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTcgcgcg                                                                                       <  2:843781/236‑1 (MQ=255)
                                                                                                                                                  aGATCTTCTGCCGGATTGTGCGACGTTTGTAATATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGt                                                                                                                                                             >  2:328534/1‑164 (MQ=255)
                                                                                                                                                    aTCTTCTGCCGGATTGTGCGCCGTTTGTAATATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTGATGGCATTCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCCAAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGGAACTGTAATTCTTCAAAGGCCCGCAAACCTTTTTTCGGGGCTTCCcggg                                                                              >  1:849119/1‑239 (MQ=255)
                                                                                                                                                              gattggggcggcgTTTGAAAACTCACGCAGAACCCTTTTTCTGCGGTTTGGGGGGTGATGGCATGAAAAATTTAAGAAAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGAGATATTAg                                                                     <  2:544640/231‑1 (MQ=255)
                                                                                                                                                                 tggggCGACTTTTGGAATACCACGTCGAAACCTTCTTCTTCAGCTTGGGCGGTGAGGGCATGAAAAACTTAAGAAAAAAACGGATTCCCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGAGATATTAGTTTg                                                                 <  2:93901/235‑1 (MQ=255)
                                                                                                                                                                    tgCGACGTTTGTAATATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGGTACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGAGATATTAGTTTGCtt                                                              >  2:651046/1‑239 (MQ=255)
                                                                                                                                                                           ttttTAATATCCCGCCAACACCTTTTTCTCCATTTTGGGCGTTGATGGAATGCAAACTTTAAAAAAAAAACGAATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGAGATATTAGTTTGCTTATTCAg                                                        <  2:738582/238‑1 (MQ=255)
                                                                                                                                                                              aaaaTACCACGGAGAAACCTTCTTTGTCAGGTTGGGCGGTGAGGGCATGAAAATCTTAAAAAAAAAGCGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGAGATATTAGTTTGCTTATTCAGCAcg                                                   <  2:1095974/237‑1 (MQ=255)
                                                                                                                                                                                     cACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGAGATATTAGTTTGCTTATTCAGCACGCGTGATa                                            >  2:491543/1‑239 (MQ=255)
                                                                                                                                                                                       cGTCGAAACCTTCTTCTTCAGCTTTGGCGGTGATGGAATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATAACCATAATGGCCTCAAAACCCGAAGAGGTCAGCGCGCGGGCAAGCTTGTTTGGCTTATCCTGTAATTCTTCAATGGCCCGGAACCCTTTTTCGCGCGCTTTCGGGGGGAAATTAGATTGGTTTTTCAGAAcgctgtaacgg                                          >  1:744682/1‑231 (MQ=255)
                                                                                                                                                                                            aaaCCTTCTTCTTCAGCTTTGGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAAAACCTTTTCGCGCGCTTCCGGGGAGCTATTAGTTTGCTTATTCAGCACCCGGGTTACgtcg                                       >  2:139133/1‑234 (MQ=255)
                                                                                                                                                                                                  tcttcttcAGCTTTGGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGGGCGCGGGCCAGCTTGGTTGGCTGATACTGTAATTCTTCAATGGCCCGGAACACTTTTTCGCGCGCTTTCGGGGTGATATTAGTTTGCTTATTCAGACcgcggggatatgggggttttgt                               >  1:339293/1‑221 (MQ=255)
                                                                                                                                                                                                  tcttcttcAGCTTTGGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGAGATATTAGTTTGCTTATTCAGCACGCCGGATACGGTGGAttttta                               >  2:196407/1‑237 (MQ=255)
                                                                                                                                                                                                   tttgttcAGCGTGGGGGGTGATGCAATACAAATCTTAAGAGAAAACCGGATTCCCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGAGATATTAGTTTGCTTATTCAGCACGCGTGATACGGTGGATTTTGAt                              <  2:744682/235‑1 (MQ=255)
                                                                                                                                                                                                   tttcttcCGTTTGGGGGGTGAGGGGATGAAAAACTTCGAAAAAAACCGGTTCCCCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGAGATATTAGTTTGCTTATTCAGCACGCGTGATACGGTGGATTTTGAt                              <  2:647738/238‑1 (MQ=255)
                                                                                                                                                                                                          aGCTTTGGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGAGATATTAGTTTGCTTATTCAGCACGCGTGATACGGTGGATTTTGATACg                           >  2:901573/1‑236 (MQ=255)
                                                                                                                                                                                                           gCTTTAGCGGTGATGGCATGCAAAACTTCAGAGAAAACTCGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGAATCAAAAACCGAAGGAGTCAGGGCGCGGGCCATCTTGTTTGGCTGATACTGTAATTATTCAATTGCCTGATCCACATTTTCTGGGGCTTCCGGGGGGGTCTTGATTTGCTTTTTAGCCAATGCGGTTACGAGGGATTTTGATAGGCCTGc                      >  1:278121/1‑239 (MQ=255)
                                                                                                                                                                                                                     ggtgtcctgaaaaacCTTAAGAAAAAAACTTATTCCCCGCCGTAGTGTTGGGCGAACGGGTAAAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGAGATATTAGTTTGCTTATTCAGCACGCGTGATACGGTGGATTTTGATACGCCTGCAACCCTGGCa            <  1:651046/228‑1 (MQ=255)
                                                                                                                                                                                                                      gATCGCATGCAAACCTTCAGAGAAAACAGAATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATGCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGAGATATTAGTTTGCTTATTCAGCACGCGTGATACGGTGGATTTTGATACGCCTGCAACCCTGGCa            <  2:261425/238‑1 (MQ=255)
                                                                                                                                                                                                                       tgtgggttgaaaaATTGAGAAAAAAATGGGTTTCCCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGAGATATTAGTTTGCTTATTCAGCACGCGTGATACGGTGGATTTTGATACGCCTGCAACCCTGGCa            <  1:552028/229‑1 (MQ=255)
                                                                                                                                                                                                                        gggCATACAAAGCTCAAGAGAAAACCGGTGTTCCCGCCGGAGTTTGGGTCGAACGGGTAGAAACCCCCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGAGATATTAGTTTGCTTATTCAGCACGCGTGATACGGTGGATTTTGATACGCCTGCAACCCTGGCAata         <  2:602586/238‑1 (MQ=255)
                                                                                                                                                                                                                         ggaaTGAAAACCTTAAGAAAAAACGGGATTCCCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGAGATATTAGTTTGCTTATTCAGCACGCGTGATACGGTGGATTTTGATACGCCTGCAACCCTGGCAata         <  2:339293/235‑1 (MQ=255)
                                                                                                                                                                                                                         ggaaTGAAAAACTTAAGAAAAAAACGAATCCCCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGAGATATTAGTTTGCTTATTCAGCACGCGTGATACGGTGGATTTTGATACGCCTGCAACCCTGGCAata         <  2:1291302/235‑1 (MQ=255)
                                                                                                                                                                                                                             tgataacgttCAAAGAAAAACGGTTTCTCCGCCGTACTTTTGGTACAACGGTGAGAAATCACGAAAATGGCTTCAAAACCCGAAAAGGTCAGCGCGCGGGACTGGTTGTTTGGCTGATGCTGTAATTATTCAATTGCCCGTAACACTTTTTCGTGAACTTCCGGGGAGATATTAGTTTGCTTATTCAGCACGCGTGATACGGTGGATTTTGATACGGCTGCAACCCTGGAAATAGCAt     <  1:58352/230‑1 (MQ=255)
                                                                                                                                                                                                                               cAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCTCGCGCTTCCGGGGAGATATTAGTTTGCTTATTCCGCACGCGCGATCCGGGGGCTTTTGATACGCCTGCAACCCTGGCAaaatagaaaa  >  1:1007842/1‑233 (MQ=255)
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                            
GGCATATTGACCTTCAACTTTGCCAATCACCACCACCGGAATATCATATTTATCGAGTTGGGCAAAAAATGACTCATCCGCTGGCGAACTGAGCATAATAATGCCTTTAATCATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGATTGTGCGACGTTTGTAATATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGAGATATTAGTTTGCTTATTCAGCACGCGTGATACGGTGGATTTTGATACGCCTGCAACCCTGGCAATATCATAAA  >  NC_000913/1383488‑1383948

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