Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,535,073 A→G intergenic (‑1020/‑1524) topAI ← / → ythB toxin of the TopAI‑YjhQ toxin‑antitoxin system, TopA inhibitor/protein YthB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,535,0730AG91.4% 141.6 / 4.9 59intergenic (‑1020/‑1524)topAI/ythBtoxin of the TopAI‑YjhQ toxin‑antitoxin system, TopA inhibitor/protein YthB
Reads supporting (aligned to +/- strand):  ref base A (1/4);  new base G (27/26);  total (28/31)
Fisher's exact test for biased strand distribution p-value = 3.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.91e-01

CTGGCAGCGTAATGCAGGTAGCCCAGTTCAACTCTCACTACCAGTATGATCCCAAGTTCGAGCGCGGAATGTATCTGTATGAGCATCGTCGCTGTTTCAATAACATCATCGACTACTGTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAACCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGCAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCT  >  NC_000913/4534828‑4535314
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                 
cTGGCAGCGTAATGCAGGTAGCCCAGTTCAACTCTCACTACCAGTATGATCCCAAGTTCGAGCGCGGAATGTATCTGTATGAGCATCGTCGCTGTTTCAATAACATCATCGACTACTGTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCctc                                                                                                                                                                                                                                              >  1:471365/1‑251 (MQ=255)
  ggCAGCGTAATGCAGGTAGCCCAGTTCAACTCTCACTACCAGTATGATCCCAAGTTCGAGCGCGGAATGTATCTGTATGAGCATCGTCGCTGTTTCAATAACATCATCGACTACTGTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAACCCTCtt                                                                                                                                                                                                                                            <  2:275381/251‑1 (MQ=255)
    cAGCGTAATGCAGGTAGCCCAGTTCAACTCTCACTACCAGTATGATCCCAAGTTCGAGCGCGGAATGTATCTGTATGAGCATCGTCGCTGTTTCAATAACATCATCGACTACTGTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTg                                                                                                                                                                                                                                           >  1:675856/1‑250 (MQ=255)
     aGCGTAATGCAGTTATCCCAGCTCATCTCTCACTTCCAGTATGATCCCAAGTTCGAGCGCGGAATGTATCTGTATGTGCATCGTCGCTGTTTCAATAACATCATCGACTACTGTAACTCACTGTGCGATCACGGGAAGGTGCAGACTAAAAGAGTGATTGAAAAAGGACGATTTTCCCCGTAATGGGATATATACATAGTGATGGTAGAGGCATGAAGCCAAATGGTGGGAGCCGGCAGAGCCCTCTTGaa                                                                                                                                                                                                                                         <  2:552057/251‑1 (MQ=255)
              cAGGTAGCCCAGTTCAACTCTCACTACCAGTATGATCCCAAGTTCGAGCGCGGAATGTATCTGTATGAGCATTGTCGCTGTTTCAATAACATCATCGACTACGGTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGATCCGTCATAGCCCTCTTGAAGCTGAAACa                                                                                                                                                                                                                                <  2:115322/251‑1 (MQ=255)
                            cAACTCTCACTACCAGTATGATCCCAAGTTCGAGCGCGGAATGTATCTGTATGAGCATCGTCGCTGTTTCAATAACATCATCGACTACTGTAACTCACTGTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATTGCAGCAGGGCt                                                                                                                                                                                                                  <  2:443871/251‑1 (MQ=255)
                                   cACTACCAGTATGATCCCAAGTTCGAGCGCGGCATGTATCTGTATGAGCATCGTCGCTGTTTCAATAACATCATCGACTACTGTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGc                                                                                                                                                                                                            >  2:105824/1‑250 (MQ=255)
                                        ccAGTATGATCCCAAGTTCGAGCGCGGAATGTATCTGTATGAGCATCGTCGCTGTTTCAATAACATCATCGACTACTGTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAACCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATaa                                                                                                                                                                                                      <  1:212183/251‑1 (MQ=255)
                                                  cccAAGTTCGAGCGCGGAATGTATCTGTATGAGCATCGTCGCTGTTTCAATAACATCATCGACTACTGTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGatat                                                                                                                                                                                             <  1:187092/250‑1 (MQ=255)
                                                   ccAAGTTCGAGCGCGGAATGTATCTGTATGAGCATCGTCGCTGTTTCAATAACATCATCGACTACTGTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGAtata                                                                                                                                                                                            >  2:341241/1‑250 (MQ=255)
                                                      aGTTCGAGCGCGGAAGGTATCTGTACGAGCATCGTCGCTGTTTCAATAACATCATCGACTACTGGAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGTGATGGAAAAAGGACGATTCTCCCCGCAATGGGATATTTACATATTGATGGGAGAGGCAGGAAGCCAAATGTAGGGAGCCGTCATAGCCCTCTTGAAGCGGAAACAATAGGAGCATGGCTGGTCGCACATAAGGATGATATAGAgc                                                                                                                                                                                        <  2:293240/251‑1 (MQ=255)
                                                               cccgAATCTATCTGTATGAGCATCGTTGCTGTTTCAATAACATCAACGACTATTGTAACTCACTTTGCTTTTCCGTTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGAGTTTCCCCGCAATGGGATATTTTCATATTGATGGGAGAGGCATGAAGCCAAATTTAGTGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATATCAGCATGGCTGTTGGCACATAAGGATGATATAGAACGCTAttccg                                                                                                                                                                               <  2:48886/248‑4 (MQ=255)
                                                                cgGAATGTATCTGTATTAGCATCGGCGCTGTTTCAATAACATCATCGACTACTTGAACTCACTTTGCTATCTCGGTAAGTTGCAGCTTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACgg                                                                                                                                                                              <  2:424981/251‑1 (MQ=255)
                                                                cgGAATGTATCTGTATGAGCATCGTCGCTGTTTCAATAACATCATCGACTACTGTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAACAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGTTGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAACAGCATGGCTGGTGGCACATAAGGATGATATAGAGGGCCagttcgg                                                                                                                                                                              >  2:14180/1‑245 (MQ=255)
                                                                  gAATGTATCTGTATGAGCATCGTCGCTGTTTCAATAACATCATCGACTACTGTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCAAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACCTAAGGATGATATAGAGCGCCATAACgggg                                                                                                                                                                            >  1:649448/1‑251 (MQ=255)
                                                                   aaTGTATCTGTATGAGCATCGTCGCTGTTTCAATAACATCATCGACTACTGTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACgggg                                                                                                                                                                            >  1:595268/1‑250 (MQ=255)
                                                                        aTCTGTATGAGCATCGTCGCTGTTTCAATAACATCATCGACTACTGTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGGGGAACATAAGGATGATATAGAGCGCAATTTCGGGGAACCAc                                                                                                                                                                      >  2:240990/1‑251 (MQ=255)
                                                                             tATGAGCATCGTCGCTGTTTCAATAACATCATCGACTACTGTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATag                                                                                                                                                                                           >  1:78552/1‑225 (MQ=255)
                                                                             tATGAGCATCGTCGCTGTGTCAATAACATCATCGACTACTGTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATag                                                                                                                                                                                           <  2:78552/225‑1 (MQ=255)
                                                                               tGAGCATCGTCGCTGTTTCAATAACATCATCGACTACTGTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCCCATAAGGATGATATAGAGCGCCATTACGGGGAACCCATGTACa                                                                                                                                                                >  2:587016/1‑250 (MQ=255)
                                                                                   cATCGTCGCTGTTTCAATAACATCATCGACTACTGTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGGCCATAAGGATGATATAGAGCGCCATTACGGGGAACCCCTGTACAAAgtt                                                                                                                                                           >  2:462902/1‑251 (MQ=255)
                                                                                    aTCGTCGCTGTTTCAATAACATCATCGACTACTGTAACTCACTTTGCTATCACGGAAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTAATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCAATTACGGGGAACCACTGTACAAAgttg                                                                                                                                                          >  1:581130/1‑251 (MQ=255)
                                                                                               ttCAATAACATCATCGACTACGGTAACTCACATTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATATCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCAGTACGGGGAACCACTGTACAAAGTTGTTGGtgttgt                                                                                                                                                <  2:661245/250‑1 (MQ=255)
                                                                                                 cAATAACATCATCGACTACTGTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGAc                                                                                                                                             <  1:364788/251‑1 (MQ=255)
                                                                                                       catcatCGACTACTGTAACTCACGTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCAtt                                                                                                                                       <  1:621560/251‑1 (MQ=255)
                                                                                                          catcGACTACTGTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCCAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGCTATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTTTGAAGCCATTTTc                                                                                                                                    >  2:589298/1‑251 (MQ=255)
                                                                                                          catcGACTACTGTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACAt                                                                                                                                                                                                        <  1:362354/183‑1 (MQ=255)
                                                                                                          catcGACTACTGTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGGGGCACAt                                                                                                                                                                                                        >  2:362354/1‑183 (MQ=255)
                                                                                                             cgaatactgttacctacgttGCTAGCACGGTAAGTTGCAGCCTAAAAGATGGATGGAAAAAGGACGATTTTCCTCGTAATGGGATAATTACATATGGATGGTAGAGGCATGAAGACAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCAGTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCgg                                                                                                                                  <  2:13859/232‑1 (MQ=255)
                                                                                                                 tactGTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGgcgt                                                                                                                               >  1:107013/1‑248 (MQ=255)
                                                                                                                     gTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTa                                                                                                                         <  1:587016/251‑1 (MQ=255)
                                                                                                                       aaCTCACTTTGCTATCACGGTAAGTTGCAGCCCAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAACCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGca                                                                                                                                                                                                           >  1:317905/1‑167 (MQ=255)
                                                                                                                       aaCTCACTTTGCTATCACGGTAAGTTGCAGCCCAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAACCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGca                                                                                                                                                                                                           <  2:317905/167‑1 (MQ=255)
                                                                                                                        aCTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATg                                                                                                                      <  2:209563/251‑1 (MQ=255)
                                                                                                                                   tATCCCGGTAAGTTTCAGCCTAAAATAGGTATAGAAAAAGGACGATTTTCCCCGCAAGGGGATAATTACATATTGATTGTAGAGGCATGAAGCCAAATGGAGTGAGCCGTCATAGCCCGCTTGAAGCTTAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAATTTATTGGTGTTGGGACGTCATTTTCTGCGTAGGTTAATGCTATCAAAAt                                                                                                            <  2:442059/250‑1 (MQ=255)
                                                                                                                                   tATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGCGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATg                                                                                                           >  1:405218/1‑251 (MQ=255)
                                                                                                                                          gTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAACCCGCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTTGTGTTGTGACGCCATTTTCGGCGCAGGTTAATGCTATCAAAATGTCAtt                                                                                                      <  2:177069/249‑1 (MQ=255)
                                                                                                                                                   aggcTAAAAGAGGGAATGAAAAAGGACGATCGTTCCCGCAATGGTATTTTTACACATGGATGGTAGAGGCATGAAGACAAATGGAGGGAGTAGTCATAGCCCTGTTGAATCGGAATCAATATTAGCATGGCTGGCGGCACAGCAGGAGGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTGCGGCGTAGGTTAATGCTAGCAAAAGGGCATTGCGTAAGCTgg                                                                                           <  2:476792/248‑1 (MQ=255)
                                                                                                                                                           agagGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCACTCTTGAAGCTGAAACAATAGCAGCATGGCTGGGGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACAACTGTACAAAGTTGTTGGTGTTGTGAAGCAATTTTCGGCGTAGGTGAATGCTATCAAAATGTAATTGAGTAGGCTGGAAAGTAAt                                                                                   >  2:328495/1‑251 (MQ=255)
                                                                                                                                                             agGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATgg                                                                                 >  1:332012/1‑251 (MQ=255)
                                                                                                                                                                   ggAAAAAGGACGATTTGCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTTTACAAAGTTGTTGGTGGTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAg                                                                            <  2:220114/250‑1 (MQ=255)
                                                                                                                                                                               aTTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGttt                                                               <  1:105824/251‑1 (MQ=255)
                                                                                                                                                                                ttttCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGga                                                                                                                                                                                                   >  1:466898/1‑118 (MQ=255)
                                                                                                                                                                                gtttCCCCGCAATGGGATATTTACATATTGATGTTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCGTGAAGCTGAAACAATAGCAGCAGGGTTGGTGGCACATAAGga                                                                                                                                                                                                   <  2:466898/117‑1 (MQ=255)
                                                                                                                                                                                       cGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGAc                                                         >  2:537858/1‑249 (MQ=255)
                                                                                                                                                                                         cAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGACCCACTGTACACAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAGGCTGAAATTTAATGGCAAAGACGAGCAGGGTTTACTGACAGTg                                                     >  2:400656/1‑251 (MQ=255)
                                                                                                                                                                                                 tatTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACtgt                                             >  1:178916/1‑251 (MQ=255)
                                                                                                                                                                                                      aCATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACt                                        <  1:537858/251‑1 (MQ=255)
                                                                                                                                                                                                          atTGATGGTAGAGGCAAGAAGCCAAATGGAGGGTGCCATCATAGCCCACGTGAAGCTGAAACAATAGCATCATGGCTGGTTGCAAATAAGGATTATATAGAGCGCCATTACGGGGAACCCCTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTAGCAAAATGTCAGTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGTTTTACTGACAGTGGGGACTGTGCACTCCtt                                    <  2:177826/251‑1 (MQ=255)
                                                                                                                                                                                                             gATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGGTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCTTga                                  <  2:391914/250‑2 (MQ=255)
                                                                                                                                                                                                              agggTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCCTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCTTGCAg                                <  2:47708/249‑1 (MQ=255)
                                                                                                                                                                                                              aTGGTAGAGGCATGAAGCCAAATGGCGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGTAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCTTGc                                  >  2:156372/1‑249 (MQ=255)
                                                                                                                                                                                                               tGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCTTGCAgg                               >  1:644050/1‑251 (MQ=255)
                                                                                                                                                                                                                            aaGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGc                  <  1:452551/251‑1 (MQ=255)
                                                                                                                                                                                                                            aaGCCAAATGGAGGGAGCAGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGc                  >  1:113520/1‑251 (MQ=255)
                                                                                                                                                                                                                                 aaaTGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTg             >  1:228209/1‑251 (MQ=255)
                                                                                                                                                                                                                                   agggagtgaggcggcATAGCCCTCTTGAAGCTGAAACAATAGCAGTAGTGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTATAAAGTTGTTGGTGTTGTGACGCCATTTTCGGTGTAGGTTAATGCTATCAAAATTTCATTGCGTAATCTGTAAATTAATGGCAAAGACTAGAAGGGTTTACTGACAGTGTGGACTGTGCCCTCCGTGCAGGGGGCGGAAAGGGCGTTTGtt           <  2:652751/237‑1 (MQ=255)
                                                                                                                                                                                                                                       agggagCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAAGGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTctct       <  1:400656/251‑1 (MQ=255)
                                                                                                                                                                                                                                            gCCGACATAGCCATCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGTAATTTCATGGCAAAGAAGAGCAGGGTTTACTGCCAGTGGGGACTGTGCACTCCTAGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTtct  >  2:126331/1‑251 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                 
CTGGCAGCGTAATGCAGGTAGCCCAGTTCAACTCTCACTACCAGTATGATCCCAAGTTCGAGCGCGGAATGTATCTGTATGAGCATCGTCGCTGTTTCAATAACATCATCGACTACTGTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAACCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGCAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCT  >  NC_000913/4534828‑4535314

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