Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,430,348 A→G 78.5% T434A (ACG→GCG)  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,3480AG78.5% 376.0 / 90.2 227T434A (ACG→GCG) stfRRac prophage; putative tail fiber protein
Reads supporting (aligned to +/- strand):  ref base A (25/24);  new base G (85/92);  total (111/116)
Fisher's exact test for biased strand distribution p-value = 7.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.62e-01

Only 100 of 227 total aligned reads displayed.

CGGCGTCATCGGCATCATCGGCGTCTGCTTCAAAAGATGAGGCGACCAGACAAGCGTCAGCAGCGAAGAGCAGCGCCACGACGGCATCCACGAAGGCGACAGAGGCTGCTGGCAGTGCGACGGCGGCAGCTCAGAGCAAAAGTACGGCGGAATCCGCGGCAACGCGCGCCGAGACAGCAGCTAAACGGGCAGAGGATATTGCATCCGCCGTGGCGCTTGAGGATGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCT  >  NC_000913/1430205‑1430493
                                                                                                                                               |                                                                                                                                                 
              tcatcGGCGTCTGCTTCAAAAGATGAGGCGACCAGACAAGCGTCAGCAGCGAAGAGCAGCGCCACGACGGCATCCACGAAGGCGACAGAGGCTGCTGGCAGTGCGACGGCGGCAGCTCAGAGCAAAAGTACGGCGGAATCCGCGGCAAc                                                                                                                                <  2:259306/149‑1 (MQ=255)
                       tCTGCTTCAAAAGATGAGGCGACCAGACAAGCGTCAGCAGCGAAGAGCAGCGCCACGACGGCATCCACGAAGGCGACAGAGGCTGCTGGCAGTGCGACGGCGGCAGCTCAGAGCAAAAGTACGGCGGAATCCGCGGCAACGCGCGCCga                                                                                                                       <  5:72719/149‑1 (MQ=255)
                           cTTCAAAAGATGAGGCGACCAGACAAGCGTCAGCAGCGACGAGCAGCGCCACGACGGCATCCACGAAGGCGACAGAGGCTGCTGGCAGTGCGACGGCGGCAGCTCAGAGCAAAAGTACGGCGGAATCCGCGGCAACGCGCGCCgaga                                                                                                                     <  7:362412/147‑1 (MQ=255)
                            ttCAAAAGATGAGGCGACCAGACAAGCGTCAGCAGCGAAGAGCAGCGCCACGACGGCATCCACGAAGGCGACAGAGGCTGCTGGCAGTGCGACGGCGGCAGCTCAGAGCAAAAGTACGGCGGAATCCGCGGCAACGCGCGCCGAGAcag                                                                                                                  >  2:66812/1‑149 (MQ=255)
                                    aTGAGGCGACCAGACAAGCGTCAGCAGCGAAGAGCAGCGCCACGCCGGCATCCACGAAGGCGACAGAGGCTGCTGGCAGTGCGACGGCGGCAGCTCAGAGCAAAAGTACGGCGGAATCCGCGGCAACGCGCGCCGAGACAGCAGCTaaa                                                                                                          >  3:281975/1‑149 (MQ=255)
                                            aCCAGACAAGCGTCAGCAGCGAAGAGCAGCGCCACGACGGCATCCACGAAGGCGACAGAGGCTGCTGGCAGTGCGACGGCGGCAGCTCAGAGCAAAAGTACGGCGGAATCCGCGGCAACGCGCGCCGAGACAGCAGCTAAACGGGCaga                                                                                                  <  6:200430/149‑1 (MQ=255)
                                             ccAGACAAGCGTCAGCAGCGAAGAGCAGCGCCACGACGGCATCCACGAAGGCGACAGAGGCTGCTGGCAGTGCGACGGCGGCAGCTCAGAGCAAAAGTACGGCGGAATCCGCGGCAACGCGCGCCGAGACAGCAGCTAAACGGGCagag                                                                                                 <  1:122124/149‑1 (MQ=255)
                                              cAGACAAGCGTCAGCAGCGAAGAGCAGCGCCACGACGGCATCCACGAAGGCGACAGAGGCTGCTGGCAGTGCGACGGCGGCAGCTCAGAGCAAAAGTACGGCGGAATCCGCGGCAACGCTCGCCGAGACAGCAGCTAAACGGGCaga                                                                                                  <  5:189697/147‑1 (MQ=255)
                                              cAGACAAGCGTCAGCAGCGAACAGCAGCGCCACGACGGCATCCACGAAGGCGACAGAGGCTGCTGGCAGTGCGACGGCGGCAGCTCACAGCAAAACTACGGCGGAATCCGCGGCAACGCGCGCCGAGACAGAAGCGAAACGGGCaga                                                                                                  >  6:189697/1‑147 (MQ=255)
                                                          agcagcGAAGAGCAGCGCCACGACGGCATCCACGAAGGCGACAGAGGCTGCTGGCAGTGCGACGGCGGCAGCTCAGAGCAAAAGTACGGCGGAATCCGCGGCAACGCGCGCCGAGACAGCAGCTAAACGGGCAGAGGATATTGCATcc                                                                                     <  8:200914/148‑1 (MQ=255)
                                                              gcGAAGAGCAGCGCCACGACGGCATCCACGAAGGCGACAGAGGCTGCTGGCAGTGCGACGGCGGCAGCTCAGAGCAAAAGTACGGCGGAATCCGCGGCAACGCGCGCCGAGACAGCAGCTAAACGGGCAGAGGATATTGCATCCGCCGt                                                                                >  2:199060/1‑149 (MQ=255)
                                                                 aaGAGCAGCGCCACGACGGCATCCACGAAGGCGACAGAGGCTGCTGGCAGTGCGACGGCGGCAGCTCAGAGCAAAAGTACGGCGGAATCCGCGGCAACGCGCGCCGAGACAGCAGCTAAACGGGCAGAGGATATTGCATCCGCCGTGgc                                                                             <  2:101571/149‑1 (MQ=255)
                                                                         tgcgtcaacgGCATCCACGAAGGCGACAGAGGCTGCTGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGTCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTAGCGCTTGAgg                                                                     <  3:26832/142‑1 (MQ=255)
                                                                         tgcgtcaacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAgg                                                                     >  3:106782/8‑149 (MQ=255)
                                                                         tgcgtcaacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAgg                                                                     <  5:126857/142‑1 (MQ=255)
                                                                         tgcgtcaacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAgg                                                                     <  1:239943/142‑1 (MQ=255)
                                                                         tgcgtcaacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAgg                                                                     <  8:65603/142‑1 (MQ=255)
                                                                           cgtcaacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGa                                                                    <  8:131061/143‑1 (MQ=255)
                                                                           cgtcaacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGa                                                                    <  7:336253/143‑1 (MQ=255)
                                                                           cgtcaacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGAt                                                                   <  2:21513/144‑1 (MQ=255)
                                                                           cgtcaacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGAt                                                                   >  5:333784/6‑149 (MQ=255)
                                                                           cgtcaacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGAt                                                                   >  5:43800/6‑149 (MQ=255)
                                                                           cgtcaacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGAt                                                                   >  7:104337/6‑149 (MQ=255)
                                                                           cgtcaacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGAt                                                                   <  7:289320/144‑1 (MQ=255)
                                                                           cgtcaacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGAt                                                                   <  2:296004/144‑1 (MQ=255)
                                                                            gtcaacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGCAGCGGCGGCAATACGTGCAGAAAATTCGGCAACACGTGCAGAAGATATAGCTTCAGCTGTCGCGCatgatgatg                                                                  >  3:275674/5‑149 (MQ=255)
                                                                            gtcaacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTcag                                                                                    <  2:394331/127‑3 (MQ=255)
                                                                            gtcaacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATg                                                                  <  1:71614/145‑1 (MQ=255)
                                                                            gtcaacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATg                                                                  >  5:199934/5‑149 (MQ=255)
                                                                            gtcaacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATg                                                                  >  1:187272/5‑149 (MQ=255)
                                                                            gtcaacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATg                                                                  >  6:232575/5‑149 (MQ=255)
                                                                            gtcaacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATg                                                                  >  4:262158/5‑149 (MQ=255)
                                                                            gtcaacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATg                                                                  >  4:158701/5‑149 (MQ=255)
                                                                            gtcaacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATg                                                                  >  7:176931/5‑149 (MQ=255)
                                                                            gtcaacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATg                                                                  >  8:168605/5‑149 (MQ=255)
                                                                            gtcaacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATg                                                                  >  7:299347/5‑149 (MQ=255)
                                                                            gtcaacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAAATCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATg                                                                  >  2:338521/5‑149 (MQ=255)
                                                                            gtcaacgGCAGCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATg                                                                  <  8:210648/145‑1 (MQ=255)
                                                                             tcaacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCTATACGTGCAGAAAATTCGGGAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATGc                                                                 <  4:309687/146‑1 (MQ=255)
                                                                             tcaacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATGc                                                                 <  8:137396/146‑1 (MQ=255)
                                                                             tcaacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATGc                                                                 <  3:309446/146‑1 (MQ=255)
                                                                             tcaacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATGc                                                                 >  8:71757/4‑149 (MQ=255)
                                                                             tcaacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATGc                                                                 <  4:51507/146‑1 (MQ=255)
                                                                             tcaacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATGc                                                                 >  6:322011/4‑149 (MQ=255)
                                                                             tcaacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATGc                                                                 >  7:276/4‑149 (MQ=255)
                                                                             tcaacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATGc                                                                 >  6:382118/4‑149 (MQ=255)
                                                                             tcaacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATGc                                                                 <  6:7000/146‑1 (MQ=255)
                                                                             tcaacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATGc                                                                 >  3:62769/4‑149 (MQ=255)
                                                                             tcaacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATGc                                                                 >  7:383252/4‑149 (MQ=255)
                                                                             tcaacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGCCAAAACGTGCAGAACATATAGCTTCAGCTGTCGCGCTTGAGGATGc                                                                 >  2:379694/4‑149 (MQ=255)
                                                                             tcaacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGAGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATGc                                                                 >  1:288016/4‑149 (MQ=255)
                                                                              caacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTca                                                                                     <  2:373445/126‑2 (MQ=255)
                                                                              caacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTca                                                                                     >  1:373445/3‑127 (MQ=255)
                                                                              caacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATGcg                                                                >  8:92663/3‑148 (MQ=255)
                                                                               aacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCTGTTTCAGAATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATGcg                                                                <  7:147754/147‑2 (MQ=255)
                                                                               aacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATGcgg                                                               >  6:236227/2‑147 (MQ=255)
                                                                               aacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATGcg                                                                <  4:121810/147‑2 (MQ=255)
                                                                               aacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATGcg                                                                <  2:139288/147‑2 (MQ=255)
                                                                               aacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTGGAGGATGcgg                                                               >  5:237279/2‑147 (MQ=255)
                                                                               aacgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGAGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATGcg                                                                >  1:344557/2‑147 (MQ=255)
                                                                                acgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATGcgga                                                              <  3:102257/149‑4 (MQ=255)
                                                                                acgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATGcgga                                                              <  4:181570/149‑4 (MQ=255)
                                                                                acgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATGcgga                                                              >  8:277905/1‑146 (MQ=255)
                                                                                acgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATGcgga                                                              <  3:329714/149‑4 (MQ=255)
                                                                                acgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATGcgga                                                              <  4:335152/149‑4 (MQ=255)
                                                                                acgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATGcgga                                                              >  8:344573/1‑146 (MQ=255)
                                                                                acgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATGcgga                                                              <  4:76675/149‑4 (MQ=255)
                                                                                acgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATGcgga                                                              >  3:220667/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)
                                                                                 cgGCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATGcggac                                                             <  3:337818/149‑5 (MQ=255)
                                                                                  ggCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGCAGCGGCGGCAATACGTGCAGAAAATTCGGCAACACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATGcggaca                                                            >  1:47683/1‑144 (MQ=255)
                                                                                  ggCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATg                                                                  >  8:34204/1‑143 (MQ=255)
                                                                                  ggCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATGcggaca                                                            <  4:303024/149‑6 (MQ=255)
                                                                                  ggCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATGcggaca                                                            <  8:41763/149‑6 (MQ=255)
                                                                                  ggCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATGcggaca                                                            <  1:141370/149‑6 (MQ=255)
                                                                                  ggCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGACAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATGcggaca                                                            <  5:109496/149‑6 (MQ=255)
                                                                                   gCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCTATACGTGCAGAAGATATAGCTTGAGCTGTCGCGCTTGAGGATGcggacac                                                           <  8:363440/149‑7 (MQ=255)
                                                                                   gCATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATGcggacac                                                           <  2:375243/149‑7 (MQ=255)
                                                                                    cATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGa                                                                       >  2:383052/1‑136 (MQ=255)
                                                                                    cATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGa                                                                       <  1:383052/136‑1 (MQ=255)
                                                                                    cATCCACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATGcggacaca                                                          <  3:151534/149‑8 (MQ=255)
                                                                                        cACGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGtc                                                                               >  7:235218/1‑123 (MQ=255)
                                                                                         aCGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGctt                                                                                       <  7:342148/115‑3 (MQ=255)
                                                                                         aCGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATg                                                                  <  5:253363/136‑1 (MQ=255)
                                                                                          cGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGatat                                                                                            >  6:269326/1‑109 (MQ=255)
                                                                                          cGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGatat                                                                                            <  5:269326/109‑1 (MQ=255)
                                                                                          cGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAGAAGCGGCGGCAATACGTGCAGAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTGTCGCGCTTGAGGATGcg                                                                >  6:309728/1‑136 (MQ=255)
                                                                                                   cAGAGGCTGCTGGCAGTGCGACGGCGGCAGCTCAGAGCAAAAGTACGGCGGAATCCGCGGCAACGCGCGCCGAGACAGCAGCTAAACGGGCAGAGGATATTGCATCCGCCGTGGCGCTTGAGGATGCAAGTa                                                            >  6:2917/1‑132 (MQ=255)
                                                                                                   cAGAGGCTGCTGGCAGTGCGACGGCGGCAGCTCAGAGCAAAAGTACGGCGGAATCCGCGGCAACGCGCGCCGAGACAGCAGCTAAACGGGCAGAGGATATTGCATCCGCCGTGGCGCTTGAGGATGCAAGTa                                                            <  5:2917/132‑1 (MQ=255)
                                                                                                            ctgGCAGTGCGACGGCGGCAGCTCAGAGCAAAAGTACGGCGGAATCCGCGGCACCGCGCGCCGAGCCAGCAGCTAACCGGGCAGAGGATATTGCCTCCGCCGTGGCGCTTGAGGATGCAAGTACGACGAAAAAGGGGATAGTACAGCTc                                  >  3:241934/1‑149 (MQ=255)
                                                                                                              ggCAGTGCGACGGCGGCAGCTCAGAGCAAAAGTACGGCGGAATCCGCGGCAACGCGCGCCGAGACAGCAGCTAAACGGGCAGAGGata                                                                                             <  8:267913/88‑1 (MQ=255)
                                                                                                               gCAGTGCGACGGCGGCAGCTCAGAGCAAAAGTACGGCGGAATCCGCGGCAACGCGCGCCGAGACAGCAGCTAAACGGGCAGAGGATATTGCATCCGCCGTGGCGCTTGAGGATGCAAGTACGACGAAAAAGGGGATAGTACAGCTcag                                >  4:353147/1‑148 (MQ=255)
                                                                                                                    gCGACGGCGGCAGCTCAGAGCAAAAGTACGGCGGAATCCGCGGCAACGCGCGCCGAGACAGCAGCTAAACGGGCAGAGGATATTGCATCCGCCGTGGCGCTTGAGGATGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGc                          >  4:277787/1‑149 (MQ=255)
                                                                                                                         ggcggcAGCTCAGAGCAAAAGTACGGCGGAATCCGCGGCAACGCGCGCCGAGACAGCAGCTAAACGGGCAGAGGATATTGCATCCGCCGTGGCGCTTGAGGATGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGAcc                      >  7:58966/1‑148 (MQ=255)
                                                                                                                                  tCAGAGCAAAAGTACGGCGGAATCCGCGTCAACGCGCGCCGAGACAGAAGCTAAACGGGAAGAGGATATTGCATCCGCCGTGGCGCTTGAGGATGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGt            >  5:308218/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)
                                                                                                                                               |                                                                                                                                                 
CGGCGTCATCGGCATCATCGGCGTCTGCTTCAAAAGATGAGGCGACCAGACAAGCGTCAGCAGCGAAGAGCAGCGCCACGACGGCATCCACGAAGGCGACAGAGGCTGCTGGCAGTGCGACGGCGGCAGCTCAGAGCAAAAGTACGGCGGAATCCGCGGCAACGCGCGCCGAGACAGCAGCTAAACGGGCAGAGGATATTGCATCCGCCGTGGCGCTTGAGGATGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCT  >  NC_000913/1430205‑1430493

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