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

GCAATGCTGATATTCGGGTTAATGCCC‑TGACCGTTAAAGCCCGTCGCGGCGATCACCACGCCGACAAAAGCCGCGCTGGAGCCGAGATAACTGGGGACGCGCCCGCCGGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGT‑TGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTCAC  >  NC_000913/1069505‑1069982
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                               
gcAATGCGGATATTGGGGTTAATGCCC‑TGACCGTGAAAGCCCGTCGCGGCGACCACCCCCCCGAAAAAACCCGCGCGGGATCCGAGATAACTGGGGACGCGCCCGCCGGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCaaa                                                                                                                                                                                                                                     <  1:464838/251‑1 (MQ=255)
          atatCGGGGTAATGCCCCTGCCCTTAAAAGCCCGCCGCGGCGACCCCCCCGCCGAAAAAAGCCGCGCTGAAGCCGAGAAAACGGGGAACGCGCCCGCCGGTAATAAAGAAAAACAGTAACGGGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCg                                                                                                                                                                                                                            <  1:515082/248‑1 (MQ=255)
               ggtttAATGCCC‑GAACCTTTAAACCCCGTCGCGGCGACCCCCCCCCCGACAAAAGCCGCGGTGGACCCGAGATATCGGGGGACGCGCCCGCCGGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATc                                                                                                                                                                                                                      <  2:468638/248‑1 (MQ=255)
                      ggCCC‑TGACCGTAAAAGCCCGTCGCGGCAATCCCCACGCCGACAAAAGCCGCGCGGGAGCCGAGATAACTGGGGACGCGCCCGCCGGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGcgcg                                                                                                                                                                                                               <  1:620940/250‑1 (MQ=255)
                               ccGTTAAAGCCCGTCGCGGCGATCACCACGCCGACACAAGCCGCGCTGGAGCCGAGATAACTGGGGACGCGCCCGCCGGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCAGTAACACCGTTGCGCCAAACATCGCCACCGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAAAGGCAATCGTTCTTCCGGCGCGACCAcgcc                                                                                                                                                                                                       >  2:508687/1‑251 (MQ=255)
                                             cgcgGCGATCACCACGCCGACAAAAGCCGCGCTGGAGCCGAGATAACTGGGGACGCGCCCGCCGGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTCAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAgg                                                                                                                                                                                         >  1:383869/1‑251 (MQ=255)
                                                             ccGACAAAAGCCGCGCTGGAGCCGAGATAACTGGGGACGCGCCCGCCGGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGcc                                                                                                                                                                          >  1:643717/1‑250 (MQ=255)
                                                                    atccggcgtgtggcgaagaaatcggggggACGCCCACGACGGTAAAAAAGAAAACCAGTATCGTCCACATCCCCGACATTAAAATGGAAATATTGGGCTCCTGCCCAATTAAAATCGGAATTAACCCCGTCGCGCCAACCATCGCCACCGCGTGT‑TGAACGCCCTTTACGGTCGTTTGAGCAAGCGGCAATTTTTCATCCGGCGGGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGa                                                                                                                                                                  <  2:326017/226‑1 (MQ=255)
                                                                          cgcTGGAGCCGAGATAACTGGGGACGCGCCCGCCGGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGCAAACCg                                                                                                                                                            >  1:145243/1‑251 (MQ=255)
                                                                           ccggcgctgaaaaaactgggggagccgcccgccGGAAAAAAAGAAAAACAGTAACGCGCCGATCCCCGACATTAAAATGGAAAGATTGGGTTCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCCCCTGTTGT‑TGAACGCTCATTCCTGTCGTCTGAGCAAACGACAATCTTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGa                                                                                                                                                           <  1:291443/226‑1 (MQ=255)
                                                                              ggAGCCGAGATAACTGGGGACGCGCCCGCCGGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAAc                                                                                                                                                         >  1:42818/1‑250 (MQ=255)
                                                                                     ataaagtgggggccctgcccgccGGAAAAAAAGAAAAACAGTACCGGGCCGATCCCCGACATTAAAATGGAAAGATTGGTATCCAGCCCCATCAGAAGCGGCATTAACCCCGTCGCGCCAAACATCGCCACCGCGTTT‑TGAACGCCCATTTCTGTCGTCTGAGCAAACGGCAAGCGTTCATCCTGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACTGAACATTGCCAt                                                                                                                                                 <  2:42673/235‑1 (MQ=255)
                                                                                     agaTAACTGGGGACGCGCCCGCCGGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCAt                                                                                                                                                 >  1:607179/1‑251 (MQ=255)
                                                                                       aaaacggggcccgcccccgccGGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCATCCCCATCAGAATCGGCATTAACACCTTCGCGCCAAACATCGCCACCGCGTTT‑TGAACGCCCATTACGGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCa                                                                                                                                               <  1:710077/236‑1 (MQ=255)
                                                                                              gagagccacccgcaGGTAAAAAAGAAAAACAGAAACGTGCCGATCCCCGACATTAAAATGTAAATATGGGGATCCAGCCCCATCAGAATCGGCATTAACCCCTCCGCGCCAAACATCGCCTCCGCGTGT‑TGAACGCCCAGTACGGTCGGCTGAGCAAACGTCAATCGTTCATCCGGCGTAACCACGCCGCTCGCTGTAGAGGTCGATTTTAACTGCCAGTTAGGAAAACCGAACATTGCTATCAGCTGtct                                                                                                                                        <  2:466736/243‑1 (MQ=255)
                                                                                                      cccgccGGTAAAAAAAAAAAACAGGAACGTGCAGATCCCCGACATTATGATGGAAAGATGGGTATCGAGCCCCAGCAAAGTCGGCATTACCACCGGCCCGCTAAACATTGCCGCCGC‑TGTGTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAACCGTTCATCCGGCGCGTCCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAGACCGAACATTGCCATCAGCTGTCTCCTTAagt                                                                                                                                <  2:71861/251‑2 (MQ=255)
                                                                                                             gTAATAAAGAAAAACAGTACCGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTa                                                                                                                          >  2:662552/1‑250 (MQ=255)
                                                                                                              taataaAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGc                                                                                                                                                                           >  2:635795/1‑200 (MQ=255)
                                                                                                              taataaAGAAAAACAGTAACGTGCCGACCCCCGACATTAAAATGGAAAGATTGGGAGCCAGCCCCATCAGAAGCGTCATTAACACCGTCCCGCCAAACACCGCCACCGCTTGT‑TGAACGCCCATTACTGCCTTCTGAGCAAACGGCAATCGTTCACCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGc                                                                                                                                                                           <  1:635795/200‑1 (MQ=255)
                                                                                                                 taaAGAAAAACAGTAACGTGCCGATCCCCGCCATTAAACTGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGCGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACa                                                                                                                       >  2:82423/1‑249 (MQ=255)
                                                                                                                  aaaGAAAAACAGTAACGTGCCGAACCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCACCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGT‑TGAACGCCCATTCCTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGc                                                                                                                    <  1:110124/251‑1 (MQ=255)
                                                                                                                      aaaaaCAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAggg                                                                                                                >  1:282228/1‑251 (MQ=255)
                                                                                                                        aaaCAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGCACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGgcg                                                                                                              >  2:126192/1‑251 (MQ=255)
                                                                                                                            aGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATc                                                                                                          >  2:189264/1‑251 (MQ=255)
                                                                                                                            aGTAACGGGCCGTTCCCCGACATAAAAAGGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACCCCGTCGCGCCAAACATCGCCCCCGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATc                                                                                                          <  2:643717/251‑1 (MQ=255)
                                                                                                                               aaCGCGCCGATCCCCGACATTAAAAGGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACCCCTTCGCGCCAAACATCGCCGCCGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGGCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAgc                                                                                                       <  2:231691/251‑1 (MQ=255)
                                                                                                                                    gCCGATCCCCGCCCTCAAAATGGAAAGATTGGTATCCAGCCACCTCAGAATCGGCATTAACACCGTCGCGCCAAACATCCCCACCGCGTGT‑TGAACGCCCCTTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGCGGTCGATTTTAACTGCCAGTGCGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTACCAAGCAGGGCGCATCAGCGCGTg                                                                                                  >  2:717795/1‑251 (MQ=255)
                                                                                                                                     ccaaaccccGACATTAAAATGGAAAAATTGGGATCCAGCCCCATCAGAACCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCTTCCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGa                                                                                                 <  2:115863/246‑1 (MQ=255)
                                                                                                                                      cGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCACACATCGCAACCGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGACAACCGAACATTGCCATCAGCTGTCTCCTTCAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGAt                                                                                                >  1:44761/1‑251 (MQ=255)
                                                                                                                                            ccGACATTAAAAGGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGT‑TGAACGCCCATTACTGGCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCTACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGc                                                                                          <  1:666895/251‑1 (MQ=255)
                                                                                                                                             cGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTCAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCg                                                                                         >  2:541913/1‑251 (MQ=255)
                                                                                                                                              gACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCg                                                                                         >  1:342507/1‑250 (MQ=255)
                                                                                                                                                   tAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGa                                                                                    >  2:252265/1‑250 (MQ=255)
                                                                                                                                                    aaaaTGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAAc                                                                                  >  2:498336/1‑251 (MQ=255)
                                                                                                                                                    aaaaTGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTAGCGCCACACATCGCCACCGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGAGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGCGGTTCACAAGCAGGGCCCCTCAGCGCGTGATAACTGCGATCGAAc                                                                                  >  2:207444/1‑251 (MQ=255)
                                                                                                                                                      aaTGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAAcc                                                                                 >  2:262033/1‑250 (MQ=255)
                                                                                                                                                           aaaGATTGGGATCCAGCCCCATCAGAATCGGCATTCACACCGTCGCGCCACACATCGCCACCGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGCAAACCGAACATTGCCCTCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGTTAACTGCGATCGAACCACAcc                                                                            >  2:659346/1‑250 (MQ=255)
                                                                                                                                                             aGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGt                                                                                                   >  1:395791/1‑225 (MQ=255)
                                                                                                                                                             aGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGACGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAACCCGAACATGGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGACCGAACCACACCAGc                                                                         >  1:713546/1‑251 (MQ=255)
                                                                                                                                                              gatggggagccgcccccATCATAAACGGCATTAACCCCGTCGCGCCAAACACCGCCACCGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGcc                                                                        <  2:261959/238‑1 (MQ=255)
                                                                                                                                                              gATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGcc                                                                        <  1:710389/251‑1 (MQ=255)
                                                                                                                                                                ttGGGATCCAACCCCTTCCGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGT‑TGAACGCCCATTACTGTCGTCTCAGCACACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGCACATTGCCATCAGCTGTCTCCTTACGGAGGTTACCAAGCAGGGCGCATCAGCGCGTGCTAACTGCGCTCGAACCACACCAGCCCg                                                                      >  2:106004/1‑251 (MQ=255)
                                                                                                                                                                 tGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGt                                                                     >  1:577416/1‑251 (MQ=255)
                                                                                                                                                                    gATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACCTCGCCACCGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGGCACTCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAgg                                                                  >  2:734426/1‑251 (MQ=255)
                                                                                                                                                                     acccAGCCCCATCAGAATCGGCATTACCCCCGTCGCGCCAAACATCGCCACCGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCTTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAAcc                                                                                 <  1:704259/233‑1 (MQ=255)
                                                                                                                                                                     aTCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAAcc                                                                                 >  2:704259/1‑235 (MQ=255)
                                                                                                                                                                        cAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTATGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGggtgg                                                              >  2:722037/1‑250 (MQ=255)
                                                                                                                                                                             ccaacggAATCTGCATTAACACCGTCGCCCCAAACACCGCCCCCGCGTGT‑TGAACGCCCATGCCTGGCGTCTGAGCAAACGGCAATCTTGCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGggtgtggtgt                                                         <  2:44761/245‑1 (MQ=255)
                                                                                                                                                                               aTCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGgtgtggtgtg                                                        >  1:550642/1‑250 (MQ=255)
                                                                                                                                                                                tCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGAc                                                      >  1:738058/1‑251 (MQ=255)
                                                                                                                                                                                ccAGAACCGGCATTAACCCCGTCGCGCCAAACACCGCCACCGCGTGT‑TGAACGCCCATTACGGTTGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACGGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGAc                                                      <  2:342507/250‑1 (MQ=255)
                                                                                                                                                                                  aGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGT‑TGAACGCCCATTACTGTCGTCTGCGCAAACGGCAATCGTTCATCCGGCGCGCCCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTCGGGTGTGGTGTGACg                                                     >  1:68781/1‑250 (MQ=255)
                                                                                                                                                                                     aTCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATg                                                  >  2:14943/1‑250 (MQ=255)
                                                                                                                                                                                     aTCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGGGTGGTGTGACGATg                                                  >  2:329686/1‑250 (MQ=255)
                                                                                                                                                                                     aTCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGCTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATg                                                  >  1:454977/1‑250 (MQ=255)
                                                                                                                                                                                     aTCGGCATTAACACAGTCGCGCCAAACATCGCCACCGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGCGTGACGATg                                                  >  1:187705/1‑250 (MQ=255)
                                                                                                                                                                                        ggCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGCAGGGTGTGGTGTGACGATGAATc                                              >  2:484418/1‑251 (MQ=255)
                                                                                                                                                                                         gCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCg                                             >  1:212105/1‑251 (MQ=255)
                                                                                                                                                                                          cATTAACACCGTCGCGCCAAACATCGCCACCGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGc                                            >  1:308330/1‑251 (MQ=255)
                                                                                                                                                                                          cATTAACACCGTCGCGCCAAACATCGCCACCGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGc                                            <  1:226838/251‑1 (MQ=255)
                                                                                                                                                                                            ttAACACCGTCGCGCCAAACATCGCCACCGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCtt                                          >  1:124141/1‑251 (MQ=255)
                                                                                                                                                                                                 acCGTCGCGCCAAACATCGCCACCGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATg                                     >  1:572775/1‑251 (MQ=255)
                                                                                                                                                                                                    gTCGCGCAAAACATCGCCCCCGCGTGT‑TGAACGCCCATTACGGTCGTCTGAGCAAACGGCAAGCGTTCACCCGGCGCGCCCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGgcg                                  <  2:308330/251‑1 (MQ=255)
                                                                                                                                                                                                           cAAACATCGCCACCGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGCCCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACa                           <  2:42818/251‑1 (MQ=255)
                                                                                                                                                                                                            aaaCATCGCCACCGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACCCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAACCCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGCCGCATCACCGCGTGATAACTGCGATCGAACCACACAAGCCCGTAGGGTTGGGGGTGACGATAACTCGCTTGGATGGCGCAAAAc                            >  2:567795/1‑249 (MQ=255)
                                                                                                                                                                                                                 ccGCCACCGCGTGT‑TGAACGCCCATTACTGGCGTCTGAGAAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGt                     <  1:252265/250‑1 (MQ=255)
                                                                                                                                                                                                                  cGCCACCGCGTGT‑TGAACGCCCATTACTGTCGTCTGCGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTc                    >  2:523689/1‑251 (MQ=255)
                                                                                                                                                                                                                  cGCCACCGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTc                    >  1:104925/1‑251 (MQ=255)
                                                                                                                                                                                                                       ccGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTggg               <  1:14943/251‑1 (MQ=255)
                                                                                                                                                                                                                        cgcgTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAgggggt                                                              >  2:455360/1‑200 (MQ=255)
                                                                                                                                                                                                                        cgcgTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGggtgt                                                              <  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)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                               
GCAATGCTGATATTCGGGTTAATGCCC‑TGACCGTTAAAGCCCGTCGCGGCGATCACCACGCCGACAAAAGCCGCGCTGGAGCCGAGATAACTGGGGACGCGCCCGCCGGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGT‑TGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTCAC  >  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: ‑