Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,558,291 C→T G8G (GGC→GGT)  rtcR → DNA‑binding transcriptional activator RtcR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,558,2910CT98.2% 175.4 / ‑6.5 55G8G (GGC→GGT) rtcRDNA‑binding transcriptional activator RtcR
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (43/10);  total (45/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

GCATTTTCAGTGGTCAGTAATTCGTAATTCATTTTGTTTTTCTCTTTTCGTTGTTTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGCTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTCATCGATCGACTGGAATTGTTGCACGACGCCCGCTCGCGCTCGCTATTTGAAACGCTTAAACGCGATATCGCCAGCGTTTCGCCAGAAACAGAAGTGGTGAGCGTTGAGATTGAACTGCATAACCCGTGGGATTTCGAAGAGGTCTACGCC  >  NC_000913/3558048‑3558537
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                      
gCATTTTCAGTGGTCAGTAATTCGTAATTCATTTTGTTTTTCTCTTTTCGTTGTTTGCTGTCCTGCTAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTACTTTTTTAATTCCTTGTTTTAACCAAGATTATATCTTTACGTCCGTAAACGGCGATTTCCAGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATATACGATGCGTACAACAGTGGCTTTTGGTTTTGTCg                                                                                                                                                                                                                                                 >  1:526082/1‑251 (MQ=255)
 cATTTTCAGTGGTCAGTAATTCGTAATTCATTTTGTTTTTCTATTTACGTTGTTTGCTGTCCTGATAATTGAACCCGTCGTGCCAGACACTTAAATAAGCAGCTTAATTTTATAATTCATTGTTTTAAAAAAGCTTATATCTTTACGTCCGTAACCGGAGATTTCACGCCACGCCCATTTACAGATAATGAAATATCGTCTTTTATAAGGATATATAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCgg                                                                                                                                                                                                                                                >  1:242214/1‑251 (MQ=255)
  aTTTTCAGTGGTCAGTAATTCGTAATTCATTTTGTTTTTCTCTTTTCGTTGTTTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGt                                                                                                                                                                                                                                               >  1:440072/1‑251 (MQ=255)
                     tcgtaattcATTTTGTTTTTCTCTTTTCGTTGTTTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCg                                                                                                                                                                                                                            >  1:90021/1‑251 (MQ=255)
                     tcgtaattcATTTTGTTTTCCTCTTTTCCTTGTTTGCTGTCCTGATAATCGCAACCGTCGTGCCAGAAAATTACATAAGCAGCTTCATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCACGCAAAGCCAATTTAACGATAATGAAATATCGTCTTTTATAAGGATATCTACGATGCGTAAAACAGTGGATTTTGGTTTTGTCGGTACCGTACTGGATTATGCCg                                                                                                                                                                                                                            >  2:575369/1‑251 (MQ=255)
                       gtaattcATTTTGTTTTTCTCTTTTCGTTGTTTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCggg                                                                                                                                                                                                                          >  1:272654/1‑251 (MQ=255)
                                      tttCTCTTTTCGTTGTTTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAgcgc                                                                                                                                                                                                           >  1:493729/1‑251 (MQ=255)
                                                 gttgtctccggTCCTGATCATTGCAACGGTCGTCCCAGACACTTAAATACGCACCTTAATTTTTTACTTCATTACCTTACAAACGCATACATCTTTACGACCGTAAACGGCGATTTCACGCACAGCCAATTTACCGATGATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGAATTTTGGATTTGTCGATACCGTACTGGATTATACCGTGAGCGACCGCCAGCGCTGGTaaaa                                                                                                                                                                                                   >  2:394849/11‑248 (MQ=255)
                                                       tgctgtggtgataacagccaCCGTCGCGCCTGAATATTACATGAGCAGCTTATATTTTGAATTCATTGTCTTAACAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGCTAATGAAATATCGTCTTTTCTAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCCACGCTGGTCAAAATggagtac                                                                                                                                                                                          >  2:90002/20‑246 (MQ=255)
                                                       tGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTcc                                                                                                                                                                                          >  1:261420/1‑251 (MQ=255)
                                                         cTGTCCTGGTAATTGCAAACGTCGTGTCCGAAAATTAAATACGCAGCTTAATTTTTTCATTCCTTGTTTTAAAAAAGATTATATCTTTACGTCCGTACCCGGAGATTTCCCGCAAAGCCAATTTACCGATTATGAAATATCGCCTTTTTTAAGGCTATCTAAGATGCGTACAACAGTGGCTTTTGGTTTATTCGGTACCGTACTGGATTATGCCGGGCGAGGCAGTCAGCGCTGGTCAAAATGGCGTCCg                                                                                                                                                                                         >  2:563554/1‑250 (MQ=255)
                                                                 aTAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCACGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGt                                                                                                                                                                                <  1:404296/251‑1 (MQ=255)
                                                                    aTTGCAAACGTCGTGCCCGAAAATTAAATAAGAAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGCTATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGACGGTACCGTACTGGATTATGCCGGGCGCGGCCGTCAGCGCTGGTCACACTGGCGTCCGACCCTCTGttt                                                                                                                                                                              >  2:335349/1‑250 (MQ=255)
                                                                        cAACCGTCGTTCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGcc                                                                                                                                                                         <  1:203668/251‑1 (MQ=255)
                                                                            cgtcgtGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGc                                                                                                                                                                      >  1:306269/1‑250 (MQ=255)
                                                                                gtGCCAGAAAATTAAATACGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTCCGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACAGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGAGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGa                                                                                                                                                                  >  1:381775/1‑250 (MQ=255)
                                                                                       aaaaTTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTgg                                                                                                                                                          <  1:523064/251‑1 (MQ=255)
                                                                                                                ttAATTCATTGTTTTAAAAAAGATTATATATTTACGTCCGTAACCGGAGATTTCAAGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTCATCGATCGACTGGCATTGTTGc                                                                                                                                  >  2:391412/1‑250 (MQ=255)
                                                                                                                      cATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGCTTTCCAGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGCCGTCCGACACTCTGTTTATGCCCGCAAGAATCGTTGGTCATCGATTGACTGGAATTGTTGCACGACGc                                                                                                                           >  2:83001/1‑251 (MQ=255)
                                                                                                                      cATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTCATCGATCGACTGGAATTGTTGCACGACGc                                                                                                                           >  2:67244/1‑251 (MQ=255)
                                                                                                                                 aaaaGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCAAATTTACCGATATTGAAATCTCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTCTGGCTTTGTCGATACCGTAATGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAAAGGCGTCCGCCACTCTGTTTATGCCAGCAAGAATCGTTGGTCATCGCTCGACTGGAATTTTAAAACGACGCCcgctagcgct                                                                                                                >  1:57054/1‑251 (MQ=255)
                                                                                                                                          tatCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGcc                                                                                                                                                                                                                             >  1:311980/1‑133 (MQ=255)
                                                                                                                                          tatCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGcc                                                                                                                                                                                                                             <  2:311980/133‑1 (MQ=255)
                                                                                                                                          tatCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTCATCGATCGACTGGAATTGTTGCACGACGCCCGCTCGCGCTCGCTATTTg                                                                                                       >  1:557602/1‑251 (MQ=255)
                                                                                                                                          tatCTTTACGTCCGGAAACGGAGATTTCCCGCAAAGCCAATTTACAGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTCATCGATCGACTGGAATTGTTGCACGACGCCCGCTCGCGCTCGCTATTTg                                                                                                       >  2:262885/1‑251 (MQ=255)
                                                                                                                                               ttACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTCATCGATCGCCTGGAATTGTTGCACGACGCCCGCTCGCGCTCGCTATTTGAAACg                                                                                                  >  1:486775/1‑251 (MQ=255)
                                                                                                                                                    tCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTCATCGATCGACTGGAATTGTTGCACGACGCCCGCTCGCGCTCGCTATTTGAAACGCTTa                                                                                              >  2:480646/1‑250 (MQ=255)
                                                                                                                                                        tAACCGGCGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTCATCGATCGACTGGAATTGTTGCACGACGCCCGCTCGCGCTCGCTATTTGAAACGCTTAAAcgc                                                                                         >  1:168293/1‑251 (MQ=255)
                                                                                                                                                        tAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTCATCGATCGACTGGAATTGTTGCACGACGCCCGCTCGCGCTCGCTATTTGAAACGCTTAAAcgc                                                                                         >  1:551538/1‑251 (MQ=255)
                                                                                                                                                              gagaTTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTCATCGATCGACTGGAATTGTTGCACGACGCCCGCTCGCGCTCGCTATTTGAAACGCTTAAACGCGATATc                                                                                   <  2:253396/251‑1 (MQ=255)
                                                                                                                                                               agaTTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTCATCGATCGACTGGAATTGTTGCACGACGCCCGCTCGCGCTCGCTATTTGAAACGCTTAAACGCGATATCg                                                                                  >  1:232926/1‑251 (MQ=255)
                                                                                                                                                               agaTTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTCAGCGATCGACTGGAATTGTTGCACGACGCCCGCGCGCGCTCGCTATTTGAAACGCTTAAACTCGATATCg                                                                                  <  2:304966/251‑1 (MQ=255)
                                                                                                                                                                              ccAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTCATCGATCGACTGGAATTGTTGCACGACGCCCGCTCGCGCTCGCTATTTGAAACGCTTAAACGCGATATCGCCAGCGTTTCGcc                                                                     >  1:137465/1‑249 (MQ=255)
                                                                                                                                                                               cAATTTACCGATACTGACATATCGTCTTTCATACGGATCTCTAAGCTGCGAAAAACAGAGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGACGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCCAGAATCGTTGGTCATCGATCGACTGGAATTGTTGCACGACGCCCGCTCGCGCTCGCTATTTGAAACGCTTAAACGCGATATCGCCAGCGTTTCGcc                                                                     >  2:550904/1‑248 (MQ=255)
                                                                                                                                                                                    tACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTAACGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTCATCGATCAACTGGAATTGTTGCACGACGCCCGCTCGAGCTCGCTATTTCAAACGCTTAAACGCGATATCGCCAGAGTTTCGCCAGAAAc                                                               >  1:383665/1‑249 (MQ=255)
                                                                                                                                                                                    tACCGATAATGAAATATAGTCTTTTATAAGGTTATCAAAGCTGCGTAAAACAGTGAATTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGACGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTCATCGATCGACTGGAATTGTTGCACGACGCCCGCTCCCGCTCGCTCTTTGAAACGCTTAAACGCGAACTAGCCAGCGTTTCGcccgaccc                                                               >  2:343944/1‑243 (MQ=255)
                                                                                                                                                                                     aCCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGGTTATGCCAGCAAGAATCGTTGGTCATCGATCGACTGGAATTGTTGCACGACGCCCGCTCGCGCTCGCTATTTGAAACGCTTAAACGCGATATCTCCAGCGTTTCGCCAGAAACAGa                                                            <  1:438849/251‑1 (MQ=255)
                                                                                                                                                                                      ccGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTCATCGATCGACTGGAATTGTTGCACGACGCCCGCTCGCGCTCGCTATTTGAAACGCTTAAACGCGATATCGCCAGCGTTTCGCCAGAAACAGa                                                            >  2:204866/1‑250 (MQ=255)
                                                                                                                                                                                         aTAATGACATATCGTCTTTTATAAGGATCTCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTCATCGATCGACTGGAATTGTTGCACGACACCCGCTCGCGCTCGCTATTTGAAACGCTTACACGCGATATCGCCAGCGTTTCGCCAGAAACAGAAgtg                                                        >  2:118222/1‑251 (MQ=255)
                                                                                                                                                                                               caaTATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTCATCGATCGACTGGAATTGTTGCACGACGCCCGCTCGCGCTCGCTATTTGAAACGCTTAAACGCGATATCGCCAGCGTTTCGCCAGAAACAGAAGTGGTGAGc                                                  >  2:448177/2‑251 (MQ=255)
                                                                                                                                                                                               aaaTATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTCATCGATCGACTGGAATTGTTGCACGACGCCCGCTCGCGCTCGCTATTTGAAACGCTTAAACGCGATATCGCCAGCGTTTCGCCAGAAACAGAAGTGGTGAGc                                                  >  1:184916/1‑251 (MQ=255)
                                                                                                                                                                                                  tatCGTCTTTTCTAAGGATATCTAAGATGTGTAAAACAGAGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTCATCGCTAGGCTGGAATTGTTGCACCACGCCCGCTCGCCCTCGCTATTAGAAACGCTTAAACGCGATATCCCCAGCGTTTCGCCAGATACAGAAGGGGTGAGAGtt                                               >  2:139913/1‑251 (MQ=255)
                                                                                                                                                                                                  tatCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTCATCGATCGACTGGAATTGTTGCACGACGCCCGCTCGCGCTCGCTATTTGAAACGCTTAAACGCGATATCGCCAGCGTTTCGCCAGAAACAGAAGTGGTGAGCGtt                                               >  1:478446/1‑251 (MQ=255)
                                                                                                                                                                                                   atCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGCTTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTCATCGATCGACTGGAATTGTTGCACGACGCCCGCTCGCGCTCGCTATTTGAAACGCTTAAACGCGATATCGCCAGCGTTTCGCCAGAAACAGAAGTGGTGAGCGTTg                                              >  1:555415/1‑251 (MQ=255)
                                                                                                                                                                                                      gTCTTTTATAAGGATATCTAAGATGCGGAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTCATCGATCGACTGGAATTGTTGCACGACGCCCGCTCGCGCTCGCTATTTGAAACGCTTAAACGCGATATCGCCAGCGTTTCGCCAGAAACAGAAGTGGTGAGCGTTg                                              >  1:114385/1‑248 (MQ=255)
                                                                                                                                                                                                         ttttATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTCATCGATCGACTGGAATTGTTGCACGACGCCCGCTCGCGCTCGCTATTTGAAACGCTTAAACGCGATATCGCCAGCGTTTCGCCAGAAACAGAAGTGGTGAGCGTTGAGATTg                                        >  1:252635/1‑251 (MQ=255)
                                                                                                                                                                                                               aaGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTCATCGATCGACTGGAATTGTTGCACGACGCCCGCTCGCGCTCGCTATTTGAAACGCTTAAACGCGATATCGCCAGCGTTTCGCCAGAAACAGAAGTGGTGAGCGTTGAGATTGAACTGc                                  <  2:227316/251‑1 (MQ=255)
                                                                                                                                                                                                                       cTAAGATGCGTAACACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTCATCGATCGACTGGAATTGTTGCACGACGCCCGCTCGCGCTCGCTCTTTGAAACGCTTAAACGCGATATCGCCAGCGTTTCGCCAGAAACAGAAGTGGTGAGCGTTGAGATTGAACTGCATAAACCg                          >  1:481722/1‑251 (MQ=255)
                                                                                                                                                                                                                         aaGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTCATCGATCGACTGGAATTGTTGCACGACGCCCGCTCGCGCTCGCTATTTGAAACGCTTAAACGCGATATCGCCAGCGTTTCGCCAGAAACAGAAGTGGTGAGCGTTGAGATTGAACTGCATAAccc                           <  2:440072/248‑1 (MQ=255)
                                                                                                                                                                                                                         aaGATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTCATCGATCGACTGGAATTGTTGCACGACGCCCGCTCGCGCTCGCTATTTGAAACGCTTAAACGCGATATCGCCAGCGTTTCGCCAGAAACAGAAGTGGTGAGCGTTGAGATTGAACTGCATAACCCGTg                        >  2:551005/1‑251 (MQ=255)
                                                                                                                                                                                                                           gATGCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTCATCGATCGACTGGAATTGTTGCACGACGCCCGCTCGCGCTCGCTATTTGAAACGCTTAAACGCGATATCGCCAGCGTTTCGCCAGAAACAGAAGTGGTGAGCGTTGAGATTGAACTGCATAACCCGTggg                      <  2:461930/251‑1 (MQ=255)
                                                                                                                                                                                                                              gCGTAAACCAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTAAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTCATCGATCGACTGGCATTGTTGCACGACGCCCGCTCGCGCTCGCTATTTGAAACGCTTAAACGCGATATCGCCAGCGTTTCGCCAGAAACAGACGAGGTGCGCGATGAGATTGAACTGCATAACCCGTGGGAtt                   >  2:188791/1‑251 (MQ=255)
                                                                                                                                                                                                                              gCGTAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTCATCGATCGACTGGAATTGTTGCACGACGCCCGCTCGCGCTCGCTATTTGAAACGCTTAAACGCGATATCGCCAGCGTTTCGCCAGAAACAGAAGTGGTGAGCGTTGAGATTGAACTGCATAACCCGTGGGAtt                   >  1:314846/1‑251 (MQ=255)
                                                                                                                                                                                                                                 tAAAACAGTGGCTTTTGGTTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTCATCGATCGACTGGAATTGTTGCACGACGCCCGCTCGCGCTCGCTATTTGAAACGCTTAAACGCGATATCGCCAGCGTTTCGCCAGAAACAGAAGTGGTGAGCGTTGAGATTGAACTGCATAACCCGTGGGATTTCg                >  1:19702/1‑251 (MQ=255)
                                                                                                                                                                                                                                      cAGTGGCTTTTGGTTTTATCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTAGGTCATCGATCGACTGGAATTGTTGCa                                                                                                                                 >  2:349322/1‑133 (MQ=255)
                                                                                                                                                                                                                                           gcttttggttttGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTCATCGATCGACTGGAATTGTTGCACGACGCCCGCTCGCGCTCGCTATTTGAAACGCTTAAACGCGATATCGCCAGCGTTTCGCCAGAAACAGAAGTGGTGAGCGTTGAGATTGAACTGCATAACCCGTGGGATTTCGAAGAGGTCTa      <  1:448177/251‑1 (MQ=255)
                                                                                                                                                                                                                                               ttggttttGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGCCCCTCTGTTTATGCCAGCAAGAATCGTTGGTCATCGATCGACTGGACTTGTTGCACGACGCCCGCTCGCGCTCGCTATTTGAAACGCTTAAACGCGATATCGCCAGCGTTTCGCCAGAAACAGAAGTGGTGAGCGTTGAGATTGACATGCATAACCCGTGGGATTTCGAAGAGGTCTACgcc  >  1:65188/1‑251 (MQ=255)
                                                                                                                                                                                                                                               ttggttttGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTCATCGATCGCCTGGAATTGTTGCACGACGCCCGCTCGCGCTCGCTATTTGAAACGCTTAAACGCGATATCGCCAGCGTTTCGCCAGAAACAGAAGTGGTGAGCGTTGAGATTGAACTGCATAACCCGTGGGATTTCGAAGAGGTCTACgcc  >  2:193260/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                      
GCATTTTCAGTGGTCAGTAATTCGTAATTCATTTTGTTTTTCTCTTTTCGTTGTTTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGCTTTGTCGGTACCGTACTGGATTATGCCGGGCGCGGCAGTCAGCGCTGGTCAAAATGGCGTCCGACACTCTGTTTATGCCAGCAAGAATCGTTGGTCATCGATCGACTGGAATTGTTGCACGACGCCCGCTCGCGCTCGCTATTTGAAACGCTTAAACGCGATATCGCCAGCGTTTCGCCAGAAACAGAAGTGGTGAGCGTTGAGATTGAACTGCATAACCCGTGGGATTTCGAAGAGGTCTACGCC  >  NC_000913/3558048‑3558537

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