Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,312,309 G→A 11.4% P228L (CCG→CTG)  alsR ← d‑allose‑inducible als operon transcriptional repressor; autorepressor; repressor of rpiR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,312,3090GA11.4% 91.6 / 3.6 35P228L (CCG→CTG) alsRd‑allose‑inducible als operon transcriptional repressor; autorepressor; repressor of rpiR
Reads supporting (aligned to +/- strand):  ref base G (13/18);  new base A (2/2);  total (15/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.24e-01

TTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGC‑TGGGCGACAGAGACAAAAAAAGCGTCCAGCAAA‑GTTAGTTGTAATATTCTTGCCGAGGCATTACG‑ACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATAATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGCGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGAAAAAGCGGGCGGCACGGTGGATCTC  >  NC_000913/4312025‑4312604
                                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                                       
tttCAGTTATTTATTCATTATAATGTTCTAACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGC‑TGGGCGACAGAGACAAAAAAAGCGTCCAGCAAA‑GTTAGTTGTAATATTCTTGCCGAGGCATTACG‑ACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATAATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATg                                                                                                                                                                                                                                                                                         <  1:295573/301‑1 (MQ=255)
tttCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGC‑TGGGCGACAGAGACAAAAAAAGCGTCCAGCAAA‑GTTAGTTGTAATATTCTTGCCGAGGCATTACG‑ACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATAATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATg                                                                                                                                                                                                                                                                                         <  2:171852/301‑1 (MQ=255)
 ttCAGATAGTTATTCATAATAAGGGTGTGACGGTAAGGGGACAGGTAACGCCGCGATGGCGGAATTACCTTTTGCTTAGTTCAGCGCGCCTGGTGAGAAGAAATCAACAAGTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGC‑TGGGCGACAGAGAC‑AAAAAAGCGTCCAGCAAA‑GTTAGTTGTAATATTCTTGCCGAGGCATTACG‑ACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATAAGCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATggg                                                                                                                                                                                                                                                                                       <  1:388086/301‑1 (MQ=255)
                                               aaCTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGC‑TGGGCGACAGAGACAAAAAAAGCGTCCAGCAAA‑GTTAGTTGTAATATTCTTGCCGAGGCATTACG‑ACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATAATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACgg                                                                                                                                                                                                                                          >  1:266269/1‑301 (MQ=255)
                                                 cTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGC‑TGGGCGACAGAGACAAAAAAAGCGTCCAGCAAA‑GTTAGTTGTAATATTCTTGCCGAGGCATTACG‑ACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATAATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGcc                                                                                                                                                                                                                                        >  2:500314/1‑301 (MQ=255)
                                                                tttATCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGC‑TGGGCGACAGAGACAAAAAAAGCGTCCAGCAAA‑GTTAGTTGTAATATTCTTGCCGAGGCATTACG‑ACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATAATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTg                                                                                                                                                                                                                         <  2:378623/300‑1 (MQ=255)
                                                                           ttgagttCAGCGCGCTTGGTGAGAAGAAATCAATAATTGGGCCGGTTTTTTTCATATTAATATTAGCTTGTTCAATGTTGAGC‑TGGGCGACAGAGACAAAAAAAGCGTCCAGCAAA‑GTTAGTTGTAAGATTCTTGCCGAGGCATTACG‑ACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATAATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGa                                                                                                                                                                                                              <  1:463366/296‑1 (MQ=255)
                                                                             tATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGC‑TGGGCGACAGAGACAAAAAAAGCGTCCAGCAAA‑GTTAGTTGTAATATTCTTGCCGAGGCATTACG‑ACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAACTAATATAATCGGCCAGTTTCGCTATCAGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAAt                                                                                                                                                                                                            >  2:513025/1‑301 (MQ=255)
                                                                              aTTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGC‑TGGGCGACAGAGACAAAAAAAGCGTCCAGCAAA‑GTTAGTTGTAATATTCTTGCCGAGGCATTACG‑ACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATAATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGaaga                                                                                                                                                                                                           >  2:463930/1‑299 (MQ=255)
                                                                                   aGCGCGGCTGGTGAGAAGAAATCAACAATTGCGTCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGC‑TGGGCGACAGAGACAAAAAAAGCGTCCAGCAAA‑GTTAGTTGTAATATTCTTGCCGAGGCATTACG‑ACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATAATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCa                                                                                                                                                                                                      <  1:499667/301‑1 (MQ=255)
                                                                                        accgacggtttttcgattttaaaTTGGGCTGTTTTTTTG‑ATATGTATAATAGCGTGTTCAATGTTTAGCGTGGGGGA‑AGAGAAAAAAAAAGCCTTCCGCAAAGGTTAGTTGAATTATT‑TTGGCGGGGCAATAAGAACTTATTAACGGGTTTTCCGGGCCT‑GTGAGCAAAATAAATAATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCa                                                                                                                                                                                                <  1:490170/280‑1 (MQ=255)
                                                                                           tGGTGAGAAGAAATCTACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGC‑TGGGCGACAGAGACAAAAAAAGCGTCCAGCAAA‑GTTAGTTGTAATATTCTTGCCGAGGCATTACG‑ACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATAATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGc                                                                                                                                                                                              <  2:550930/301‑1 (MQ=255)
                                                                                                    aaaaTCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGC‑TGGGCGACAGAGACAAAAAAAGCGTCCAGCAAA‑GTTAGTTGTAATATTCTTGCCGAGGCATTACG‑ACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATAATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACAtct                                                                                                                                                                                     <  2:167758/300‑1 (MQ=255)
                                                                                                       attaacaaTTGTGGCGGTTTTTTGCATATTAATATTAGGTTGTTTAAGGTTGAGC‑TGGGTGACAGAGACAAAAAAAGCGTCCAGCAAA‑GTTAGTTGTAATATTCTTGCCGAGGCATTACG‑ACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATAATCGGCCAGTTTCGCTATCAGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCtcct                                                                                                                                                                                  <  1:512353/298‑1 (MQ=255)
                                                                                                          aacaaTTGTGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGC‑TGGGCGACAGAGACAAAAAAAGCGTCCAGCAAA‑GTTAGTTGTAATATTCTTGCCGAGGCATTACG‑ACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATAATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTACGCCCGGAATGGGTCACTACCAGCACAACATCtccttcc                                                                                                                                                                               <  1:298301/301‑1 (MQ=255)
                                                                                                            caaTTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGC‑TGGGCGACAGAGACAAAAAAAGCGTCCAGCAAA‑GTTAGTTGTAATATTCTTGCCGAGGCATTACG‑ACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATAATCGGCCAGTTTCGCTATCGGTGAATGTTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGTTTCGCCC‑GAATTGGTTCCTACCaaaaaccgaaaaacaaaactt                                                                                                                                                                            >  2:490816/1‑281 (MQ=255)
                                                                                                                ggcgcGGTTTTTTTGCATATTATTATTTGCTTGTGCAATTTTGAGC‑TGGGCGACAGAGACAAAAAAAGCGGCCAGCAAA‑GTTAGTTGTAATATTCTTGCCGAGGCATTACG‑ACCTAAGAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATAATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTaac                                                                                                                                                                         <  1:416543/300‑1 (MQ=255)
                                                                                                                               cATATTAATATTAGCTTGTTCAATGTTGAGC‑TGGGCGACAGAGACAAAAAAAGCGTCCAGCAAA‑GTTAGTTGTAATATTCTTGCCGAGGCATTACG‑ACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATAATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGAcatc                                                                                                                                                          <  2:266537/301‑1 (MQ=255)
                                                                                                                                 tattaatattaGCTTGTTCAATGTTGAGC‑TGGGCGACAGAGACAAAAAAAGCGTCCAGCAAA‑GTTAGTTGTAATATTCTTGCCGAGGCATTACG‑ACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATAATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGAcatcaa                                                                                                                                                        >  2:294314/1‑300 (MQ=255)
                                                                                                                                                           ttgtgctcttcctttctcAAAAAAAGCGTCCAGCAAA‑GTTAGTTGTAATATTCTTGCCGAGGCATTACG‑ACCTAATAACGGCGTTTCCGGGTCTGGTGAGCAAATAATATAATCGGCCAGTTTCGCTATCAGTGAATGGTAGCTATGGGTTATACAAATATTCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGcc                                                                                                                               <  2:503218/284‑1 (MQ=255)
                                                                                                                                                                 ggCGACAGAGACAAAAAAAGCGTCCAGCAAA‑GTTAGTTGTAATATTCTTGCCGAGGCATTACG‑ACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATAATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGAATAGGCCTGACac                                                                                                                         >  2:547763/1‑300 (MQ=255)
                                                                                                                                                                   ttctgttatcAAAAAAAGCGTCCAGCAAA‑GTTAGTTGTAATATTCTTGCCGAGGCATTACG‑ACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATAATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTGCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCg                                                                                                                       <  2:230133/292‑1 (MQ=255)
                                                                                                                                                                    gACAGAGACAAAAAAAGCGTCCAGCAAA‑GTTAGTTGTAATATTCTTGCCGAGGCATTACG‑ACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATAATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCACCATGAAGCTAGCATCAGGAAAGCCCTGAccacgg                                                                                                                      >  1:303060/1‑296 (MQ=255)
                                                                                                                                                                            cAAAAAAAGCGTCCAGCAAA‑GTTAGTTGTAATATTCTTGCCGAGGCATTACG‑ACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATAATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATACTCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTAACTCACTGGTTCGCCCGGAATGGGTCACTACCAGAACAACATCTCCTTCCTGTAACACCGAAGCGGAAATCATCATGATGAGAGAATCAGGCACGGCCAGACAGCgagatcggaa                                                                                                              >  1:229890/1‑292 (MQ=255)
                                                                                                                                                                            cAAAAAAAGCGTCCAGCAAA‑GTTAGTTGTAATATTCTTGCCGAGGCATTACG‑ACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATAATCGGCCAGTTTCGCTATCAGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCCCTGGTTCGCCCGGCATGGGTCACTACCAGCAAAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGAAGTGAACATCATGATAAGcaaaaacggaagaacgtcg                                                                                                              >  1:502546/1‑283 (MQ=255)
                                                                                                                                                                              aaaaaaTCGTCCAGCAAA‑GTTAGTTGTAATATTCTTGCCGAGGCATTACG‑ACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATAATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATgc                                                                                                            <  1:293910/301‑1 (MQ=255)
                                                                                                                                                                              aaaaaaGCGTCCAGCAAA‑GTTAGTTGTAATATTCTTGCCGAGGCATTACG‑ACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATAATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAACGTACGCCAATgc                                                                                                            >  2:301192/1‑301 (MQ=255)
                                                                                                                                                                                         cAGCAAA‑GTTAGTTGTAATATTCTTGCCGAGGCATTACG‑ACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATAATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGCGCAAGAACTtg                                                                                                 >  2:585370/1‑301 (MQ=255)
                                                                                                                                                                                                           cctttccctacccgacgctcttccgatcttaaCGGCGTTTCCGGGGCTGGTGAGCAAATAATATAATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGCGCAAGAACTTGTGCTGTACATCAGcac                                                                                 <  1:259046/272‑1 (MQ=255)
                                                                                                                                                                                                                 ttgctggtggaTTTCG‑TCCTAAGAGTGGGGTTTCCGGGGCTGGTGAGCAAATAATATAATCGGCCAGTTTCGTTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCTCCGTTCTTTTTTGCCAGTTCTACGGCCGCTGTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGCGCAAGAACCTGTGCTgatgtacagcacaagttct                                                                          <  1:415153/291‑19 (MQ=255)
                                                                                                                                                                                                                     ggAGGCATTACG‑ACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATAATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGCGCAAGAACTTGTGCTGTACATCAGCACAGATAGCAttt                                                                      <  1:584646/300‑1 (MQ=255)
                                                                                                                                                                                                                      gAGGCATTACG‑ACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATAATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGAACGGCCTGACAGCGTACGCCAATGCCCAAGAAATTGTGCTGTACATCAGCAAAGAAAACAtttt                                                                     >  1:465767/1‑300 (MQ=255)
                                                                                                                                                                                                                        ggCATTACG‑ACCGAATAAAGGCGTTTCCGGGGCTGGTGAGCAAATAATATAATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTTCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGCGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGAt                                                                   <  1:547090/301‑1 (MQ=255)
                                                                                                                                                                                                                                        taaCGGCGTTTCCGGGGCTGGTGAGCAAATAATATAATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGCGCAAGAACTTGTGCTGTACATCAGCACAGAtcggaagagcacacgtctgaactcca                                                    >  2:259314/1‑276 (MQ=255)
                                                                                                                                                                                                                                                                                          cTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGCGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGAAAAAGCGGGCGGCACGGTGGAtctc  >  1:489234/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                                       
TTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGC‑TGGGCGACAGAGACAAAAAAAGCGTCCAGCAAA‑GTTAGTTGTAATATTCTTGCCGAGGCATTACG‑ACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATAATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGCGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGAAAAAGCGGGCGGCACGGTGGATCTC  >  NC_000913/4312025‑4312604

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑