Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,069,743 C→T A33T (GCA→ACA)  rutG ← pyrimidine:H(+) symporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,069,7430CT94.4% 218.1 / ‑4.5 74A33T (GCA→ACA) rutGpyrimidine:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base C (1/3);  new base T (45/22);  total (46/28)
Fisher's exact test for biased strand distribution p-value = 1.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

GCAATGCTGATATTCGGGTTAATGCCCTGACCGTTAAAGCCCGTCGCGGCGATCACCACGCCGACAAAAGCCGCGCTGGAGCCGAGATAACTGGGGACGCGCCCGCCGGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTCAC  >  NC_000913/1069505‑1069982
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                               
gcAATGCGGATATTGGGGTTAATGCCCTGACCGTGAAAGCCCGTCGCGGCGACCACCCCCCCGAAAAAACCCGCGCGGGATCCGAGATAACTGGGGACGCGCCCGCCGGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCaaa                                                                                                                                                                                                                                     <  1:464838/251‑1 (MQ=255)
          atatCGGGGTAATGCCCCTGCCCTTAAAAGCCCGCCGCGGCGACCCCCCCGCCGAAAAAAGCCGCGCTGAAGCCGAGAAAACGGGGAACGCGCCCGCCGGTAATAAAGAAAAACAGTAACGGGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCg                                                                                                                                                                                                                            <  1:515082/248‑1 (MQ=255)
               ggtttAATGCCCGAACCTTTAAACCCCGTCGCGGCGACCCCCCCCCCGACAAAAGCCGCGGTGGACCCGAGATATCGGGGGACGCGCCCGCCGGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATc                                                                                                                                                                                                                      <  2:468638/248‑1 (MQ=255)
                      ggCCCTGACCGTAAAAGCCCGTCGCGGCAATCCCCACGCCGACAAAAGCCGCGCGGGAGCCGAGATAACTGGGGACGCGCCCGCCGGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGcgcg                                                                                                                                                                                                               <  1:620940/250‑1 (MQ=255)
                               ccGTTAAAGCCCGTCGCGGCGATCACCACGCCGACACAAGCCGCGCTGGAGCCGAGATAACTGGGGACGCGCCCGCCGGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCAGTAACACCGTTGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAAAGGCAATCGTTCTTCCGGCGCGACCAcgcc                                                                                                                                                                                                       >  2:508687/1‑251 (MQ=255)
                                             cgcgGCGATCACCACGCCGACAAAAGCCGCGCTGGAGCCGAGATAACTGGGGACGCGCCCGCCGGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTCAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAgg                                                                                                                                                                                         >  1:383869/1‑251 (MQ=255)
                                                             ccGACAAAAGCCGCGCTGGAGCCGAGATAACTGGGGACGCGCCCGCCGGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGcc                                                                                                                                                                          >  1:643717/1‑250 (MQ=255)
                                                                    atccggcgtgtggcgaagaaatcggggggACGCCCACGACGGTAAAAAAGAAAACCAGTATCGTCCACATCCCCGACATTAAAATGGAAATATTGGGCTCCTGCCCAATTAAAATCGGAATTAACCCCGTCGCGCCAACCATCGCCACCGCGTGTTGAACGCCCTTTACGGTCGTTTGAGCAAGCGGCAATTTTTCATCCGGCGGGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGa                                                                                                                                                                  <  2:326017/226‑1 (MQ=255)
                                                                          cgcTGGAGCCGAGATAACTGGGGACGCGCCCGCCGGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGCAAACCg                                                                                                                                                            >  1:145243/1‑251 (MQ=255)
                                                                           ccggcgctgaaaaaactgggggagccgcccgccGGAAAAAAAGAAAAACAGTAACGCGCCGATCCCCGACATTAAAATGGAAAGATTGGGTTCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCCCCTGTTGTTGAACGCTCATTCCTGTCGTCTGAGCAAACGACAATCTTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGa                                                                                                                                                           <  1:291443/226‑1 (MQ=255)
                                                                              ggAGCCGAGATAACTGGGGACGCGCCCGCCGGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAAc                                                                                                                                                         >  1:42818/1‑250 (MQ=255)
                                                                                     ataaagtgggggccctgcccgccGGAAAAAAAGAAAAACAGTACCGGGCCGATCCCCGACATTAAAATGGAAAGATTGGTATCCAGCCCCATCAGAAGCGGCATTAACCCCGTCGCGCCAAACATCGCCACCGCGTTTTGAACGCCCATTTCTGTCGTCTGAGCAAACGGCAAGCGTTCATCCTGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACTGAACATTGCCAt                                                                                                                                                 <  2:42673/235‑1 (MQ=255)
                                                                                     agaTAACTGGGGACGCGCCCGCCGGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCAt                                                                                                                                                 >  1:607179/1‑251 (MQ=255)
                                                                                       aaaacggggcccgcccccgccGGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCATCCCCATCAGAATCGGCATTAACACCTTCGCGCCAAACATCGCCACCGCGTTTTGAACGCCCATTACGGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCa                                                                                                                                               <  1:710077/236‑1 (MQ=255)
                                                                                              gagagccacccgcaGGTAAAAAAGAAAAACAGAAACGTGCCGATCCCCGACATTAAAATGTAAATATGGGGATCCAGCCCCATCAGAATCGGCATTAACCCCTCCGCGCCAAACATCGCCTCCGCGTGTTGAACGCCCAGTACGGTCGGCTGAGCAAACGTCAATCGTTCATCCGGCGTAACCACGCCGCTCGCTGTAGAGGTCGATTTTAACTGCCAGTTAGGAAAACCGAACATTGCTATCAGCTGtct                                                                                                                                        <  2:466736/243‑1 (MQ=255)
                                                                                                      cccgccGGTAAAAAAAAAAAACAGGAACGTGCAGATCCCCGACATTATGATGGAAAGATGGGTATCGAGCCCCAGCAAAGTCGGCATTACCACCGGCCCGCTAAACATTGCCGCCGCTGTGTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAACCGTTCATCCGGCGCGTCCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAGACCGAACATTGCCATCAGCTGTCTCCTTAagt                                                                                                                                <  2:71861/251‑2 (MQ=255)
                                                                                                             gTAATAAAGAAAAACAGTACCGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTa                                                                                                                          >  2:662552/1‑250 (MQ=255)
                                                                                                              taataaAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGc                                                                                                                                                                           >  2:635795/1‑200 (MQ=255)
                                                                                                              taataaAGAAAAACAGTAACGTGCCGACCCCCGACATTAAAATGGAAAGATTGGGAGCCAGCCCCATCAGAAGCGTCATTAACACCGTCCCGCCAAACACCGCCACCGCTTGTTGAACGCCCATTACTGCCTTCTGAGCAAACGGCAATCGTTCACCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGc                                                                                                                                                                           <  1:635795/200‑1 (MQ=255)
                                                                                                                 taaAGAAAAACAGTAACGTGCCGATCCCCGCCATTAAACTGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGCGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACa                                                                                                                       >  2:82423/1‑249 (MQ=255)
                                                                                                                  aaaGAAAAACAGTAACGTGCCGAACCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCACCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTCCTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGc                                                                                                                    <  1:110124/251‑1 (MQ=255)
                                                                                                                      aaaaaCAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAggg                                                                                                                >  1:282228/1‑251 (MQ=255)
                                                                                                                        aaaCAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGCACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGgcg                                                                                                              >  2:126192/1‑251 (MQ=255)
                                                                                                                            aGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATc                                                                                                          >  2:189264/1‑251 (MQ=255)
                                                                                                                            aGTAACGGGCCGTTCCCCGACATAAAAAGGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACCCCGTCGCGCCAAACATCGCCCCCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATc                                                                                                          <  2:643717/251‑1 (MQ=255)
                                                                                                                               aaCGCGCCGATCCCCGACATTAAAAGGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACCCCTTCGCGCCAAACATCGCCGCCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGGCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAgc                                                                                                       <  2:231691/251‑1 (MQ=255)
                                                                                                                                    gCCGATCCCCGCCCTCAAAATGGAAAGATTGGTATCCAGCCACCTCAGAATCGGCATTAACACCGTCGCGCCAAACATCCCCACCGCGTGTTGAACGCCCCTTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGCGGTCGATTTTAACTGCCAGTGCGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTACCAAGCAGGGCGCATCAGCGCGTg                                                                                                  >  2:717795/1‑251 (MQ=255)
                                                                                                                                     ccaaaccccGACATTAAAATGGAAAAATTGGGATCCAGCCCCATCAGAACCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCTTCCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGa                                                                                                 <  2:115863/246‑1 (MQ=255)
                                                                                                                                      cGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCACACATCGCAACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGACAACCGAACATTGCCATCAGCTGTCTCCTTCAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGAt                                                                                                >  1:44761/1‑251 (MQ=255)
                                                                                                                                            ccGACATTAAAAGGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGGCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCTACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGc                                                                                          <  1:666895/251‑1 (MQ=255)
                                                                                                                                             cGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTCAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCg                                                                                         >  2:541913/1‑251 (MQ=255)
                                                                                                                                              gACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCg                                                                                         >  1:342507/1‑250 (MQ=255)
                                                                                                                                                   tAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGa                                                                                    >  2:252265/1‑250 (MQ=255)
                                                                                                                                                    aaaaTGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAAc                                                                                  >  2:498336/1‑251 (MQ=255)
                                                                                                                                                    aaaaTGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTAGCGCCACACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGAGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGCGGTTCACAAGCAGGGCCCCTCAGCGCGTGATAACTGCGATCGAAc                                                                                  >  2:207444/1‑251 (MQ=255)
                                                                                                                                                      aaTGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAAcc                                                                                 >  2:262033/1‑250 (MQ=255)
                                                                                                                                                           aaaGATTGGGATCCAGCCCCATCAGAATCGGCATTCACACCGTCGCGCCACACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGCAAACCGAACATTGCCCTCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGTTAACTGCGATCGAACCACAcc                                                                            >  2:659346/1‑250 (MQ=255)
                                                                                                                                                             aGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGt                                                                                                   >  1:395791/1‑225 (MQ=255)
                                                                                                                                                             aGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGACGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAACCCGAACATGGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGACCGAACCACACCAGc                                                                         >  1:713546/1‑251 (MQ=255)
                                                                                                                                                              gatggggagccgcccccATCATAAACGGCATTAACCCCGTCGCGCCAAACACCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGcc                                                                        <  2:261959/238‑1 (MQ=255)
                                                                                                                                                              gATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGcc                                                                        <  1:710389/251‑1 (MQ=255)
                                                                                                                                                                ttGGGATCCAACCCCTTCCGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTCAGCACACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGCACATTGCCATCAGCTGTCTCCTTACGGAGGTTACCAAGCAGGGCGCATCAGCGCGTGCTAACTGCGCTCGAACCACACCAGCCCg                                                                      >  2:106004/1‑251 (MQ=255)
                                                                                                                                                                 tGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGt                                                                     >  1:577416/1‑251 (MQ=255)
                                                                                                                                                                    gATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACCTCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCACTCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAgg                                                                  >  2:734426/1‑251 (MQ=255)
                                                                                                                                                                     acccAGCCCCATCAGAATCGGCATTACCCCCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCTTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAAcc                                                                                 <  1:704259/233‑1 (MQ=255)
                                                                                                                                                                     aTCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAAcc                                                                                 >  2:704259/1‑235 (MQ=255)
                                                                                                                                                                        cAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTATGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGggtgg                                                              >  2:722037/1‑250 (MQ=255)
                                                                                                                                                                             ccaacggAATCTGCATTAACACCGTCGCCCCAAACACCGCCCCCGCGTGTTGAACGCCCATGCCTGGCGTCTGAGCAAACGGCAATCTTGCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGggtgtggtgt                                                         <  2:44761/245‑1 (MQ=255)
                                                                                                                                                                               aTCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGgtgtggtgtg                                                        >  1:550642/1‑250 (MQ=255)
                                                                                                                                                                                tCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGAc                                                      >  1:738058/1‑251 (MQ=255)
                                                                                                                                                                                ccAGAACCGGCATTAACCCCGTCGCGCCAAACACCGCCACCGCGTGTTGAACGCCCATTACGGTTGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACGGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGAc                                                      <  2:342507/250‑1 (MQ=255)
                                                                                                                                                                                  aGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGCGCAAACGGCAATCGTTCATCCGGCGCGCCCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTCGGGTGTGGTGTGACg                                                     >  1:68781/1‑250 (MQ=255)
                                                                                                                                                                                     aTCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATg                                                  >  2:14943/1‑250 (MQ=255)
                                                                                                                                                                                     aTCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGGGTGGTGTGACGATg                                                  >  2:329686/1‑250 (MQ=255)
                                                                                                                                                                                     aTCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGCTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATg                                                  >  1:454977/1‑250 (MQ=255)
                                                                                                                                                                                     aTCGGCATTAACACAGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGCGTGACGATg                                                  >  1:187705/1‑250 (MQ=255)
                                                                                                                                                                                        ggCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGCAGGGTGTGGTGTGACGATGAATc                                              >  2:484418/1‑251 (MQ=255)
                                                                                                                                                                                         gCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCg                                             >  1:212105/1‑251 (MQ=255)
                                                                                                                                                                                          cATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGc                                            >  1:308330/1‑251 (MQ=255)
                                                                                                                                                                                          cATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGc                                            <  1:226838/251‑1 (MQ=255)
                                                                                                                                                                                            ttAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCtt                                          >  1:124141/1‑251 (MQ=255)
                                                                                                                                                                                                 acCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATg                                     >  1:572775/1‑251 (MQ=255)
                                                                                                                                                                                                    gTCGCGCAAAACATCGCCCCCGCGTGTTGAACGCCCATTACGGTCGTCTGAGCAAACGGCAAGCGTTCACCCGGCGCGCCCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGgcg                                  <  2:308330/251‑1 (MQ=255)
                                                                                                                                                                                                           cAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGCCCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACa                           <  2:42818/251‑1 (MQ=255)
                                                                                                                                                                                                            aaaCATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACCCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAACCCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGCCGCATCACCGCGTGATAACTGCGATCGAACCACACAAGCCCGTAGGGTTGGGGGTGACGATAACTCGCTTGGATGGCGCAAAAc                            >  2:567795/1‑249 (MQ=255)
                                                                                                                                                                                                                 ccGCCACCGCGTGTTGAACGCCCATTACTGGCGTCTGAGAAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGt                     <  1:252265/250‑1 (MQ=255)
                                                                                                                                                                                                                  cGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGCGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTc                    >  2:523689/1‑251 (MQ=255)
                                                                                                                                                                                                                  cGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTc                    >  1:104925/1‑251 (MQ=255)
                                                                                                                                                                                                                       ccGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTggg               <  1:14943/251‑1 (MQ=255)
                                                                                                                                                                                                                        cgcgTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAgggggt                                                              >  2:455360/1‑200 (MQ=255)
                                                                                                                                                                                                                        cgcgTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGggtgt                                                              <  1:455360/203‑1 (MQ=255)
                                                                                                                                                                                                                                tGAACGACCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCTTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGATTCGCTGGCGCAAAACAGAATGTCGTGGGTGCCGACg       >  2:349703/1‑251 (MQ=255)
                                                                                                                                                                                                                                     gCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAGTGTCGTGGGTGCCGACGCTcac  >  1:605651/1‑251 (MQ=255)
                                                                                                                                                                                                                                     ccccATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTcac  <  1:64442/250‑1 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                               
GCAATGCTGATATTCGGGTTAATGCCCTGACCGTTAAAGCCCGTCGCGGCGATCACCACGCCGACAAAAGCCGCGCTGGAGCCGAGATAACTGGGGACGCGCCCGCCGGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTCAC  >  NC_000913/1069505‑1069982

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