Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 265,588 T→C 5.7% intergenic (‑45/+32) ykfF ← / ← ykfB CP4‑6 prophage; uncharacterized protein/CP4‑6 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913265,5880TC 5.7% 345.9 / 9.0 127intergenic (‑45/+32)ykfF/ykfBCP4‑6 prophage; uncharacterized protein/CP4‑6 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (61/54);  new base C (3/4);  total (66/61)
Fisher's exact test for biased strand distribution p-value = 7.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.74e-01

Only 100 of 128 total aligned reads displayed.

ACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGACGTAGACGTCGAACTGCCCGGCCTCAGA  >  NC_000913/265439‑265736
                                                                                                                                                     |                                                                                                                                                    
aCTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTATCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGt                                                                                                                                                      >  1:1885992/1‑150 (MQ=255)
    cccTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAACTGTCggg                                                                                                                                                  >  1:1493454/1‑150 (MQ=255)
    cccTGGTCGTCTTCGAGTGGGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGACAGCTAGCAGGCTGCCggg                                                                                                                                                  <  2:1284599/150‑1 (MQ=255)
    cccTGGTCGTCATCGAGTGTGATATTGCTGTACGTGGTAGTGAAAGCTTGCGCTTGTCTCCGGGTGAAAGGACCAGGTGGCAGCAACACGCATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCAGCTAGCTGTCggg                                                                                                                                                  >  2:758634/1‑150 (MQ=255)
    cccTGGTCGGCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCggg                                                                                                                                                  <  2:1163510/150‑1 (MQ=255)
    cccTGGGAGTCTTCGAGTGTTATATTGCTGTAAGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGCAAGGCCCCGGTGTAAGCAACACGGATTGGGTCATTTTCGGGAGCCTGATAAAAGAAAACCCCGACAGCTAGCATGCTGCCggg                                                                                                                                                  <  2:1426520/150‑1 (MQ=255)
     ccTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGTAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGACAGCTAGCAGGCTGCCgggg                                                                                                                                                 <  2:1007505/150‑1 (MQ=255)
     ccTGGTCGTCTTCGAGTGTGATATTGCTGTACGGGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCgggg                                                                                                                                                 <  1:1019878/150‑1 (MQ=255)
      cTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGt                                                                                                                                                <  1:1143091/150‑1 (MQ=255)
      cTGGTCGTCATCGAGTGGGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGGCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGt                                                                                                                                                <  1:78040/150‑1 (MQ=255)
       tGGTCGTCTTCGTGTGTGATATTGCTGTACGGGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCTTTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTg                                                                                                                                               <  1:2024893/150‑1 (MQ=255)
         gtcgtcTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGaa                                                                                                                                             >  1:358554/1‑150 (MQ=255)
           cgtcTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGGCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAAtt                                                                                                                                           <  2:780663/150‑1 (MQ=255)
                tCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTg                                                                                                                                      >  2:841242/1‑150 (MQ=255)
                 cGAGTGTGATATTGCTGTACGTTGTAGTGACCGCTTGCGCTTGTCTCAGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTCTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTgg                                                                                                                                     <  2:406434/150‑1 (MQ=255)
                        gATAGTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGGGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGAt                                                                                                                              <  2:1791680/150‑1 (MQ=255)
                           atTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGAtact                                                                                                                           >  1:867630/1‑150 (MQ=255)
                           atTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGAtact                                                                                                                           >  2:781638/1‑150 (MQ=255)
                            ttGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGAtacta                                                                                                                          >  1:504860/1‑150 (MQ=255)
                             tGCTGTCCGTGGTAGTGACCGCTTGCGATTGTCTACGGGTGAAAGGCCCCTGTGGCAGCAACACGGATTGGGTCATTTTCGGGATCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGGCGGGGTGAATTTGCTGGGGAAGAtactac                                                                                                                         <  2:384902/150‑1 (MQ=255)
                             tGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGAtactac                                                                                                                         <  2:1532509/150‑1 (MQ=255)
                              gCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTATCGGGCTCCTGATAAAAGCAAACCCCGGCAGCCTGCTAGCTGGCGGGGGGAATTTGCTGGGGAAGATactact                                                                                                                        >  2:122806/1‑150 (MQ=255)
                                     gTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTc                                                                                                                 <  2:504860/150‑1 (MQ=255)
                                      tGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCg                                                                                                                <  2:820508/150‑1 (MQ=255)
                                         tAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGGGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTTATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTtg                                                                                                             <  2:107183/150‑1 (MQ=255)
                                            tGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTtgctg                                                                                                          <  1:2020037/150‑1 (MQ=255)
                                             gACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTtgctgc                                                                                                         <  1:2201315/150‑1 (MQ=255)
                                             gACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTtgctgc                                                                                                         >  2:1522897/1‑150 (MQ=255)
                                             gACCGCGTGCGCTTGTCGCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCGCCTGATAAAAGAAAACCCCGGCAGGCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGGtgcggc                                                                                                         <  1:2037448/150‑1 (MQ=255)
                                              cccGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCa                                                                                                        <  2:622274/149‑1 (MQ=255)
                                               ccGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAg                                                                                                       <  1:1032979/150‑1 (MQ=255)
                                                 gcttgcgcttgTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGtc                                                                                                     >  2:900894/1‑150 (MQ=255)
                                                 gcttgcgcttgTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGCCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGtc                                                                                                     <  1:1443923/150‑1 (MQ=255)
                                                    tgcgcttgTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTcc                                                                                                  <  1:830875/150‑1 (MQ=255)
                                                      cgcttgTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCga                                                                                                >  1:1596463/1‑150 (MQ=255)
                                                      cgcttgTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTAGTTGCTGCAGTCTCCga                                                                                                >  2:1047918/1‑150 (MQ=255)
                                                       gcttgTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCgag                                                                                               <  1:773776/150‑1 (MQ=255)
                                                          tttctcCGGGTGAAAGGCCCCGGGGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGt                                                                                            <  1:293515/148‑1 (MQ=255)
                                                            tctcCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTgg                                                                                          >  1:1318944/1‑150 (MQ=255)
                                                             ctcCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAAATTCCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGa                                                                                         >  2:226616/1‑150 (MQ=255)
                                                             ctcCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGGCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGa                                                                                         <  1:1195388/150‑1 (MQ=255)
                                                             ctcCGGGTGAAACGCCCCGGTGGCAACAACACGGATTGGGTAATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCAGCAAGCTGACGGGCTGAATTTGCTCGGGAAGATACTACTATCAGTCCTTGCTGCAGTCTCCGAGAGTGGa                                                                                         >  1:1655300/1‑150 (MQ=255)
                                                              tcCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCAGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGAc                                                                                        <  1:162823/150‑1 (MQ=255)
                                                              tcCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCACGGCAGCCTGCTAGCTGTCGGGGAGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGCAc                                                                                        >  2:1759386/1‑150 (MQ=255)
                                                               ccGGGTGAAAGGCCCCGGTGGCAGCAACCCGGATTGGGTCATTTTCGGGCTCCTGAGAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACa                                                                                       >  1:146884/1‑150 (MQ=255)
                                                                     gAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAATGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACtt                                                                                 >  2:118515/1‑150 (MQ=255)
                                                                     gAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACtt                                                                                 >  2:1585929/1‑150 (MQ=255)
                                                                      aaaGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACttt                                                                                >  2:224205/1‑150 (MQ=255)
                                                                       aaGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTtc                                                                               >  2:889116/1‑150 (MQ=255)
                                                                        aGGCCCCGGTGGCAGCAACACGGTTTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCTGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTtct                                                                              <  1:550991/150‑1 (MQ=255)
                                                                        aGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTtct                                                                              >  1:1040665/1‑150 (MQ=255)
                                                                           ccccGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAg                                                                           <  2:568328/150‑1 (MQ=255)
                                                                               ggTGGCAGCGACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGtt                                                                       >  2:900392/1‑150 (MQ=255)
                                                                               ggTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGtt                                                                       >  1:1971113/1‑150 (MQ=255)
                                                                                gTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGttc                                                                      <  2:544442/150‑1 (MQ=255)
                                                                                 tGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGttct                                                                     >  2:1131478/1‑150 (MQ=255)
                                                                                  ggCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGttctt                                                                    <  2:2015389/150‑1 (MQ=255)
                                                                                    cagcaACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTcc                                                                  <  2:2184072/150‑1 (MQ=255)
                                                                                    cagcaACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTcc                                                                  <  2:2175659/150‑1 (MQ=255)
                                                                                    cagcaACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTcc                                                                  >  1:1372478/1‑150 (MQ=255)
                                                                                    cagcaACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTcc                                                                  <  1:781638/150‑1 (MQ=255)
                                                                                    cagcaACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTcc                                                                  >  2:1418763/1‑150 (MQ=255)
                                                                                    cagcaACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCCGGCTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTcc                                                                  >  1:1196752/1‑150 (MQ=255)
                                                                                             ggATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAg                                                         <  1:1138535/150‑1 (MQ=255)
                                                                                                ttGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGa                                                      <  1:1322073/150‑1 (MQ=255)
                                                                                                  gggTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAg                                                    >  2:1134873/1‑150 (MQ=255)
                                                                                                   ggTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGGCGGGGTGAATTTGCTGGGGATGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAgg                                                   <  2:1684881/150‑1 (MQ=255)
                                                                                                    gTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGt                                                  >  1:1633872/1‑150 (MQ=255)
                                                                                                      cATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTcc                                                >  1:1378049/1‑150 (MQ=255)
                                                                                                                cTCCTGATAAAAGAAACCCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTtgtg                                      <  2:1093908/150‑1 (MQ=255)
                                                                                                                cTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTtgtg                                      <  2:2161374/150‑1 (MQ=255)
                                                                                                                 tCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTgg                                     <  2:74939/150‑1 (MQ=255)
                                                                                                                 tCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTgg                                     >  2:461966/1‑150 (MQ=255)
                                                                                                                 tCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTgg                                     <  2:358554/150‑1 (MQ=255)
                                                                                                                   cTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTc                                   >  1:218693/1‑150 (MQ=255)
                                                                                                                   cTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTc                                   >  2:1429222/1‑150 (MQ=255)
                                                                                                                    tGATAAAAGAAAACCCCGGCAGCCTCCTAGATGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCa                                  >  2:475364/1‑150 (MQ=255)
                                                                                                                     gATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCAt                                 >  1:2140554/1‑150 (MQ=255)
                                                                                                                       tAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTg                               >  1:1762293/1‑150 (MQ=255)
                                                                                                                         aaaGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGAc                             >  2:1052388/1‑150 (MQ=255)
                                                                                                                         aaaGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGAc                             <  1:483042/150‑1 (MQ=255)
                                                                                                                         aaaGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGAc                             >  2:582280/1‑150 (MQ=255)
                                                                                                                          aaGAAAACCCCGGCTGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTGTATCAGCGTTCTTCCGGTCCGCAGTGATGGTCGCTGCTTTGTGGTCATTGACg                            >  1:1058601/1‑150 (MQ=255)
                                                                                                                            gAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGACGTa                          <  2:2076535/150‑1 (MQ=255)
                                                                                                                            gAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGACGTa                          >  1:1185893/1‑150 (MQ=255)
                                                                                                                             aaacCCCCGGCCGCCTGCTTGCTGTGGGGGTGAATTTGCTGGTGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGACGTAg                         <  2:1089695/150‑1 (MQ=255)
                                                                                                                              aaaCCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGACGTAgc                        >  1:1863228/1‑149 (MQ=255)
                                                                                                                               aaCCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGACGTAGAc                       <  1:977911/150‑1 (MQ=255)
                                                                                                                                aCCCCGACAGCTAGCAGGCTGCCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGATGGTCCCTGCTTTGTGGTCATTGACGTAGACg                      >  1:1107450/1‑150 (MQ=255)
                                                                                                                                aCCCCGACAGCTAGCAGGCTGCCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGACGTAGACg                      >  1:810770/1‑150 (MQ=255)
                                                                                                                                 ccccGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGACGTAGACGt                     >  1:907117/1‑150 (MQ=255)
                                                                                                                                  cccGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACAATCAGTCGTTCCTGCAGTCTCCGAGAGTGGTCAGCACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGACGTAGACGTc                    >  1:2047394/1‑150 (MQ=255)
                                                                                                                                   ccGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGACGTAGACGTCg                   >  1:523729/1‑150 (MQ=255)
                                                                                                                                    cGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAATATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGACGTAGACGTCGa                  <  1:841242/150‑1 (MQ=255)
                                                                                                                                         gCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAACGTCCCTGCTTTGTGGTCATTGACGTAGACGTCGAACTGc             >  1:2123684/1‑150 (MQ=255)
                                                                                                                                          ccTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGACGTAGACGTCGAACTGcc            >  1:1125735/1‑150 (MQ=255)
                                                                                                                                            tGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTTGTCATTGACGTAGACGTCGAACTGCACg          <  2:734728/150‑1 (MQ=255)
                                                                                                                                            tGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGACGTAGACGTCGAACTGCCCg          >  2:468807/1‑150 (MQ=255)
                                                                                                                                                aGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGACGTAGACGTCGAACTGCCCGGCCt      >  2:1129452/1‑150 (MQ=255)
                                                                                                                                                    gTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGACGTAGACGTCGAACTGCCCGGCCTCaga  >  2:1735846/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGACGTAGACGTCGAACTGCCCGGCCTCAGA  >  NC_000913/265439‑265736

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