Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,430,368 G→A 79.2% T440T (ACG→ACA) ‡ stfR → Rac prophage; putative tail fiber protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,430,3680GA79.2% 377.9 / 82.1 222T440T (ACG→ACA) ‡stfRRac prophage; putative tail fiber protein
Reads supporting (aligned to +/- strand):  ref base G (24/22);  new base A (85/90);  total (109/113)
Fisher's exact test for biased strand distribution p-value = 7.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01

Only 100 of 222 total aligned reads displayed.

TCTGCTTCAAAAGATGAGGCGACCAGACAAGCGTCAGCAGCGAAGAGCAGCGCCACGACGGCATCCACGAAGGCGACAGAGGCTGCTGGCAGTGCGACGGCGGCAGCTCAGAGCAAAAGTACGGCGGAATCCGCGGCAACGCGCGCCGAGACAGCAGCTAAACGGGCAGAGGATATTGCATCCGCCGTGGCGCTTGAGGATGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAG  >  NC_000913/1430228‑1430513
                                                                                                                                            |                                                                                                                                                 
tCTGCTTCAAAAGATGAGGCGACCAGACAAGCGTCAGCAGCGAAGAGCAGCGCCACGACGGCATCCACGAAGGCGACAGAGGCTGCTGGCAGTGCGACGGCGGCAGCTCAGAGCAAAAGTACGGCGGAATCCGCGGCAACGCGCGCCga                                                                                                                                           <  5:72719/149‑1 (MQ=255)
    cTTCAAAAGATGAGGCGACCAGACAAGCGTCAGCAGCGACGAGCAGCGCCACGACGGCATCCACGAAGGCGACAGAGGCTGCTGGCAGTGCGACGGCGGCAGCTCAGAGCAAAAGTACGGCGGAATCCGCGGCAACGCGCGCCgaga                                                                                                                                         <  7:362412/147‑1 (MQ=255)
     ttCAAAAGATGAGGCGACCAGACAAGCGTCAGCAGCGAAGAGCAGCGCCACGACGGCATCCACGAAGGCGACAGAGGCTGCTGGCAGTGCGACGGCGGCAGCTCAGAGCAAAAGTACGGCGGAATCCGCGGCAACGCGCGCCGAGAcag                                                                                                                                      >  2:66812/1‑149 (MQ=255)
        aaaaGATGAGGCGACCAGACAAGCGTCAGCAGCGAAGAGCAGCGCCACGACGGCATCCACGAAGGCGACAGAGGCTGCTGGCAGTGCGACGGCGGCAGCTCAGAGCAAAAGTACGGCGGAATCCGCGGCAACGCGCGCCGAGAcagcag                                                                                                                                   >  5:199783/1‑149 (MQ=255)
                      ccAGACAAGCGTCAGCAGCGAAGAGCAGCGCCACGACGGCATCCACGAAGGCGACAGAGGCTGCTGGCAGTGCGACGGCGGCAGCTCAGAGCAAAAGTACGGCGGAATCCGCGGCAACGCGCGCCGAGACAGCAGCTAAACGGGCagag                                                                                                                     <  1:122124/149‑1 (MQ=255)
                       cAGACAAGCGTCAGCAGCGAAGAGCAGCGCCACGACGGCATCCACGAAGGCGACAGAGGCTGCTGGCAGTGCGACGGCGGCAGCTCAGAGCAAAAGTACGGCGGAATCCGCGGCAACGCTCGCCGAGACAGCAGCTAAACGGGCaga                                                                                                                      <  5:189697/147‑1 (MQ=255)
                       cAGACAAGCGTCAGCAGCGAAGAGCAGCGCCACGACGGAATCCACGAAGGCGACAGAGGCTGCTGGCAGTGCGACGGCGGCAGCTCAGAGCAAAAGTACGGCGGAATCCGCGGCAACGCGCGCCGAGACAGCAGCTAAACGGGCAGAgg                                                                                                                    <  6:199783/149‑1 (MQ=255)
                                   agcagcGAAGAGCAGCGCCACGACGGCATCCACGAAGGCGACAGAGGCTGCTGGCAGTGCGACGGCGGCAGCTCAGAGCAAAAGTACGGCGGAATCCGCGGCAACGCGCGCCGAGACAGCAGCTAAACGGGCAGAGGATATTGCATcc                                                                                                         <  8:200914/148‑1 (MQ=255)
                                          aaGAGCAGCGCCACGACGGCATCCACGAAGGCGACAGAGGCTGCTGGCAGTGCGACGGCGGCAGCTCAGAGCAAAAGTACGGCGGAATCCGCGGCAACGCGCGCCGAGACAGCAGCTAAACGGGCAGAGGATATTGCATCCGCCGTGgc                                                                                                 <  2:101571/149‑1 (MQ=255)
                                                  tgcgtcaacgGCATCCACGAAGGCGACAGAGGCTGCTGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGTCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTAGCGCTTGAgg                                                                                         <  3:26832/142‑1 (MQ=255)
                                                  tgcgtcaacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAgg                                                                                         <  5:126857/142‑1 (MQ=255)
                                                  tgcgtcaacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAgg                                                                                         <  3:8874/142‑1 (MQ=255)
                                                  tgcgtcaacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAgg                                                                                         <  8:65603/142‑1 (MQ=255)
                                                  tgcgtcaacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGCGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAgg                                                                                         <  6:323602/142‑1 (MQ=255)
                                                  tgcgtcaacgGCATACACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAg                                                                                          >  3:327226/8‑148 (MQ=255)
                                                    cgtcaacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGa                                                                                        <  8:131061/143‑1 (MQ=255)
                                                    cgtcaacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGa                                                                                        <  7:336253/143‑1 (MQ=255)
                                                    cgtcaacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGAt                                                                                       >  8:164404/6‑149 (MQ=255)
                                                    cgtcaacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGAt                                                                                       <  2:296004/144‑1 (MQ=255)
                                                    cgtcaacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGAt                                                                                       <  2:335215/144‑1 (MQ=255)
                                                    cgtcaacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGAt                                                                                       >  5:333784/6‑149 (MQ=255)
                                                    cgtcaacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGAt                                                                                       <  7:289320/144‑1 (MQ=255)
                                                     gtcaacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGCAGCGGCGGCAATACGTGCAGAAAATTCGGCAACACGTGCAGAAGATATAGCTTCAGCTGTCGCGCatgatgatg                                                                                      >  3:275674/5‑149 (MQ=255)
                                                     gtcaacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTcag                                                                                                        >  1:394331/5‑129 (MQ=255)
                                                     gtcaacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTcag                                                                                                        <  2:394331/127‑3 (MQ=255)
                                                     gtcaacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGag                                                                                       <  5:262017/144‑2 (MQ=255)
                                                     gtcaacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGAt                                                                                       >  6:262017/5‑148 (MQ=255)
                                                     gtcaacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATg                                                                                      >  4:262158/5‑149 (MQ=255)
                                                     gtcaacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATg                                                                                      <  1:305636/145‑1 (MQ=255)
                                                     gtcaacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATg                                                                                      >  6:232575/5‑149 (MQ=255)
                                                     gtcaacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATg                                                                                      >  1:187272/5‑149 (MQ=255)
                                                     gtcaacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATg                                                                                      >  8:43592/5‑149 (MQ=255)
                                                     gtcaacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATg                                                                                      >  7:58040/5‑149 (MQ=255)
                                                     gtcaacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATg                                                                                      >  7:299347/5‑149 (MQ=255)
                                                     gtcaacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATg                                                                                      <  3:145777/145‑1 (MQ=255)
                                                     gtcaacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATg                                                                                      >  7:176931/5‑149 (MQ=255)
                                                     gtcaacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGCGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATg                                                                                      <  5:360804/145‑1 (MQ=255)
                                                     gtcaacgGCAGCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATg                                                                                      <  8:210648/145‑1 (MQ=255)
                                                     cACGACGGCATCCACGAAGGCGACAGAGGCTGCTGGCAGTGCGACGGCGGCAGCTCAGAGCAAAAGTACGGCGGAATCCGCGGCAACGCGCGCCGAGACAGCAGCTAAACGGGCAGAGGATATTGCATCCGCCGTGGCGCTTGAGGATg                                                                                      <  1:66812/149‑1 (MQ=255)
                                                      tcaacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATGc                                                                                     >  3:188594/4‑149 (MQ=255)
                                                      tcaacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATGc                                                                                     <  6:7000/146‑1 (MQ=255)
                                                      tcaacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATGc                                                                                     <  3:334423/146‑1 (MQ=255)
                                                      tcaacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATGc                                                                                     >  7:383252/4‑149 (MQ=255)
                                                      tcaacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATGc                                                                                     <  5:249025/146‑1 (MQ=255)
                                                      tcaacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATGc                                                                                     >  8:71757/4‑149 (MQ=255)
                                                      tcaacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATGc                                                                                     >  6:382118/4‑149 (MQ=255)
                                                       caacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTca                                                                                                         >  1:373445/3‑127 (MQ=255)
                                                       caacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTca                                                                                                         <  2:373445/126‑2 (MQ=255)
                                                       caacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATGcg                                                                                    >  8:92663/3‑148 (MQ=255)
                                                        gacgGCATCCACGAAGGCGACAGAGGCTGCTGGCAGTGCGACGGCGGCAGCTCAGAGCAAAAGTACGGCGGAATCCGCGGCAACGCGCGCCGAGACAGCAGCTAAACGGGCAGAGGATATTGCATCCGCCGTGGCGCTTGAGGATGCaa                                                                                   <  4:304380/149‑1 (MQ=255)
                                                        aacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCTGTTTCAGAATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATGcg                                                                                    <  7:147754/147‑2 (MQ=255)
                                                        aacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCTGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATGcgg                                                                                   <  4:158916/148‑3 (MQ=255)
                                                        aacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATGcgg                                                                                   >  3:227985/2‑147 (MQ=255)
                                                        aacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATGcgg                                                                                   >  8:278426/2‑147 (MQ=255)
                                                        aacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATGcgg                                                                                   >  6:236227/2‑147 (MQ=255)
                                                        aacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATGcgg                                                                                   >  8:19811/2‑147 (MQ=255)
                                                        aacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATGcgg                                                                                   >  2:210807/2‑147 (MQ=255)
                                                        aacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATGcg                                                                                    <  3:103559/147‑2 (MQ=255)
                                                        aacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATGcg                                                                                    >  4:310155/2‑147 (MQ=255)
                                                        aacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTAGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTAGCGCTTGAGGATGcg                                                                                    >  7:75955/2‑147 (MQ=255)
                                                         acgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATGcgga                                                                                  <  4:335152/149‑4 (MQ=255)
                                                         acgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATGcgga                                                                                  <  2:380539/149‑4 (MQ=255)
                                                         acgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATGcgga                                                                                  >  7:137663/1‑146 (MQ=255)
                                                         acgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATGcgga                                                                                  >  8:344573/1‑146 (MQ=255)
                                                          cgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCgc                                                                                                 <  5:76211/133‑1 (MQ=255)
                                                          cgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCgc                                                                                                 >  6:76211/1‑133 (MQ=255)
                                                          cgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATGcggac                                                                                 <  1:78147/149‑5 (MQ=255)
                                                          cgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATGcgga                                                                                  >  2:17408/1‑145 (MQ=255)
                                                          cgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATGcgga                                                                                  >  4:153528/1‑145 (MQ=255)
                                                          cgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATGcgga                                                                                  >  3:307353/1‑145 (MQ=255)
                                                           ggCATCCACGAAGGCGACAGAGGCTGCTGGCAGTGCGACGGCGGCAGCTCAGAGCAAAAGTACGGCGGAATCCGCGGCAACGCGCGCCGAGACAGCAGCTAAACGGGCAGAGGa                                                                                                                   >  8:269362/1‑114 (MQ=255)
                                                           ggCATCCACGAAGGCGACAGAGGCTGCTGGCAGTGCGACGGCGGCAGCTCAGAGCAAAAGTACGGCGGAATCCGCGGCAACGCGCGCCGAGACAGCAGCTAAACGGGCAGAGGa                                                                                                                   <  7:269362/114‑1 (MQ=255)
                                                           ggCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATg                                                                                      <  7:34204/143‑1 (MQ=255)
                                                           ggCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATg                                                                                      >  8:34204/1‑143 (MQ=255)
                                                           ggCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATGcggaca                                                                                <  1:141370/149‑6 (MQ=255)
                                                           ggCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATGcggaca                                                                                <  2:376673/149‑6 (MQ=255)
                                                           ggCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATGcggaca                                                                                >  5:215415/1‑144 (MQ=255)
                                                           ggCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGACAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATGcggaca                                                                                <  5:109496/149‑6 (MQ=255)
                                                            gCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCTATACGTGCAGAAGATATAGCTTGAGCTGTCGCGCTTGAGGATGcggacac                                                                               <  8:363440/149‑7 (MQ=255)
                                                            gCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATGcggacac                                                                               <  2:93215/149‑7 (MQ=255)
                                                            gCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATGcggacac                                                                               <  1:38084/149‑7 (MQ=255)
                                                            gCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATGcggaca                                                                                >  5:108253/1‑143 (MQ=255)
                                                             cATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGa                                                                                           >  2:383052/1‑136 (MQ=255)
                                                               tCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCt                                                                                              <  4:54563/131‑1 (MQ=255)
                                                               tCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCt                                                                                              >  3:54563/1‑131 (MQ=255)
                                                                ccACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGtc                                                                                                   >  2:215712/1‑124 (MQ=255)
                                                                ccACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGtc                                                                                                   <  1:215712/125‑2 (MQ=255)
                                                                 cACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGtc                                                                                                   >  7:235218/1‑123 (MQ=255)
                                                                   cGAAGGCGACAGAGGCTGCTGGCAGTGCGACGGCGGCAGCTCAGAGCAAAAGTACGGCGTAATCCGCGGCAACGCGCGCCGAGACAGCAGCTAAACGGGCAGAGGATATTGCATCCGCCGTGGCGCTTGAGGATGCAAGTACGACGaaa                                                                        <  4:198679/149‑1 (MQ=255)
                                                                   cGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATGcg                                                                                    >  6:309728/1‑136 (MQ=255)
                                                                     aaGGCGACAGAGGCTGCTGGCAGTGCGACGGCGGCAGCTCAGAGCAAAAGTACGGCGGAATCCGCGGCAACGCGCGCCGAGACAGCAGCTAAACGGGCAGAGGATATTGCATCCGCCGTGGCGCTTGAGGATGCAAGTACGACGaaaaa                                                                      <  4:281975/149‑1 (MQ=255)
                                                                            cAGAGGCTGCTGGCAGTGCGACGGCGGCAGCTCAGAGCAAAAGTACGGCGGAATCCGCGGCAACGCGCGCCGAGACAGCAGCTAAACGGGCAGAGGATATTGCATCCGCCGTGGCGCTTGAGGATGCAAGTa                                                                                >  6:2917/1‑132 (MQ=255)
                                                                            cAGAGGCTGCTGGCAGTGCGACGGCGGCAGCTCAGAGCAAAAGTACGGCGGAATCCGCGGCAACGCGCGCCGAGACAGCAGCTAAACGGGCAGAGGATATTGCATCCGCCGTGGCGCTTGAGGATGCAAGTa                                                                                <  5:2917/132‑1 (MQ=255)
                                                                                     ctgGCAGTGCGACGGCGGCAGCTCAGAGCAAAAGTACGGCGGAATCCGCGGCACCGCGCGCCGAGCCAGCAGCTAACCGGGCAGAGGATATTGCCTCCGCCGTGGCGCTTGAGGATGCAAGTACGACGAAAAAGGGGATAGTACAGCTc                                                      >  3:241934/1‑149 (MQ=255)
                                                                                        gCAGTGCGACGGCGGCAGCTCAGAGCAAAAGTACGGCGGAATCCGCGGCAACGCGCGCCGAGACAGCAGCTAAACGGGCAGAGGATATTGCATCCGCCGTGGCGCTTGAGGATGCAAGTACGACGAAAAAGGGGATAGTACAGCTcag                                                    >  4:353147/1‑148 (MQ=255)
                                                                                                           tCAGAGCAAAAGTACGGCGGAATCCGCGTCAACGCGCGCCGAGACAGAAGCTAAACGGGAAGAGGATATTGCATCCGCCGTGGCGCTTGAGGATGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGt                                >  5:308218/1‑149 (MQ=255)
                                                                                                           tCAGAGCAAAAGTACGGCGGAATCCGCGGCAACGCGCGCCGAGACAGCAGCTAAACGGGCAGAGGATATTGCATCCGCCGTGGCGCTTGAGGATGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGt                                >  6:182187/1‑149 (MQ=255)
                                                                                                             agagCAAAAGTACGGCGGAATCCGCGGCAACGCGCGCCGAGACCGCAGCTAAACGGGCAGAGGATATTGCATCCGCCGTGGCGCTTGAGGATGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAAAAGTACGTCt                              >  5:315493/1‑149 (MQ=255)
                                                                                                             agagCAAAAGTACGGCGGAATCCGCGGCAACGCGCGCCGAGACAGCAGCTAAACGGGCAGAGGATATTGCATCCGCCGTGGCGCTTGAGGATGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCt                              <  1:199060/149‑1 (MQ=255)
                                                                                                              gagCAAAAGTACGGCGGAATCCGCGGCAACGCGCGCCGAGACAGCAGCTAAACGGGCAGAGGATATTGCATCCGCCGTGGCGCTTGAGGATGCAAGTACGACGAAAACGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCt                              >  2:116082/1‑148 (MQ=255)
                                                                                                                                            |                                                                                                                                                 
TCTGCTTCAAAAGATGAGGCGACCAGACAAGCGTCAGCAGCGAAGAGCAGCGCCACGACGGCATCCACGAAGGCGACAGAGGCTGCTGGCAGTGCGACGGCGGCAGCTCAGAGCAAAAGTACGGCGGAATCCGCGGCAACGCGCGCCGAGACAGCAGCTAAACGGGCAGAGGATATTGCATCCGCCGTGGCGCTTGAGGATGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAG  >  NC_000913/1430228‑1430513

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