Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 265,578 C→T 6.1% intergenic (‑35/+42) 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,5780CT 6.1% 411.9 / 9.7 131intergenic (‑35/+42)ykfF/ykfBCP4‑6 prophage; uncharacterized protein/CP4‑6 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (57/66);  new base T (3/5);  total (60/71)
Fisher's exact test for biased strand distribution p-value = 7.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.46e-01

Only 100 of 131 total aligned reads displayed.

ACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGACGTAGACGTCGAACTGC  >  NC_000913/265430‑265725
                                                                                                                                                    |                                                                                                                                                   
aCGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCt                                                                                                                                                    <  2:1017740/150‑1 (MQ=255)
aCGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTGGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGGCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCt                                                                                                                                                    <  1:1485293/150‑1 (MQ=255)
        gACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTg                                                                                                                                            <  1:1620196/150‑1 (MQ=255)
         aCTGCCCTGGTCGTCTTCGAGTGTGATATAGCTGTACGTGGTAGTGACCGCATGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCCGCCTGCTAGCTGt                                                                                                                                           >  1:297654/1‑150 (MQ=255)
           tGCCCTGGTCGTCTTCGAGTGGGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCg                                                                                                                                         <  2:1177509/150‑1 (MQ=255)
            gCCCTGGTCGTCTTCGAGTGTGATATTTCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCCGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCgg                                                                                                                                        >  1:1413235/1‑150 (MQ=255)
            gCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCgg                                                                                                                                        >  1:1074137/1‑150 (MQ=255)
             cccTGGTCGTCTTCTATTGTGATATTGCTGTACGTTGTAGTGACCGCTTGCGCTTGTCTACGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGACAGCTAGCATGCTGCCggg                                                                                                                                       <  1:1209733/150‑1 (MQ=255)
             cccTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCGACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCggg                                                                                                                                       >  1:34356/1‑150 (MQ=255)
             cccTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGACAGCTAGCAGGCTGCCggg                                                                                                                                       <  1:767786/150‑1 (MQ=255)
              ccTGGTCATGTTCTAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGAGCTTGTCTCCTGGTGAAATGCCCCGGTGGCAGAAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGACAGCTAGCAGGCTGCCgggg                                                                                                                                      <  2:1543171/150‑1 (MQ=255)
              catGGTCGTCTTCGCGTGTGATATTGCGGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGTCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCgggg                                                                                                                                      <  1:381145/148‑1 (MQ=255)
               cTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTGGCGCTTGTCTCCGGGTGAAAGGCCCCGGGGGCAGCAACACGGATTGGGGCATTTTCGGGCTCCTGATAAAAGAATACCCCGACAGCTAGCAGGCTGCCGGGGt                                                                                                                                     <  1:1371089/150‑1 (MQ=255)
                tGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTg                                                                                                                                    <  1:771859/150‑1 (MQ=255)
                 ggTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGAAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGa                                                                                                                                   >  2:1719938/1‑150 (MQ=255)
                     gtcTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGACAGCTAGCAGGCTGCCGGGGttttctt                                                                                                                               <  1:1624047/150‑7 (MQ=255)
                      gcTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCAGGGTGAAAGTCCCCGGTGGCATCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTg                                                                                                                              <  1:74120/149‑1 (MQ=255)
                       cTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGc                                                                                                                             >  2:166756/1‑150 (MQ=255)
                        ttCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGTGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCt                                                                                                                            <  2:1153484/150‑1 (MQ=255)
                          cGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTTAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTgg                                                                                                                          <  2:951091/150‑1 (MQ=255)
                          cGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTgg                                                                                                                          <  1:1339013/150‑1 (MQ=255)
                           gAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTggg                                                                                                                         <  1:479891/150‑1 (MQ=255)
                            aGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCGACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTgggg                                                                                                                        <  2:34356/150‑1 (MQ=255)
                             gtgtGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGa                                                                                                                       >  2:306309/1‑150 (MQ=255)
                               gtgATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAg                                                                                                                     <  1:166756/150‑1 (MQ=255)
                                tgATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATATAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAgc                                                                                                                    >  1:267502/1‑149 (MQ=255)
                                 gATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGa                                                                                                                                   >  2:163912/1‑134 (MQ=255)
                                 gATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGa                                                                                                                                   <  1:163912/134‑1 (MQ=255)
                                  atatTGCTGTCCGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGAta                                                                                                                  <  2:1166634/150‑1 (MQ=255)
                                  atatTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGAta                                                                                                                  <  2:1576049/150‑1 (MQ=255)
                                      tGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGAtactac                                                                                                              <  2:1074137/150‑1 (MQ=255)
                                       gCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATAGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATactact                                                                                                             >  2:1418031/1‑150 (MQ=255)
                                         tGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATc                                                                                                          <  2:1661075/150‑1 (MQ=255)
                                             cGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGt                                                                                                       >  2:143842/1‑150 (MQ=255)
                                                ggTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGt                                                                                                    <  2:540253/150‑1 (MQ=255)
                                                 gTAGTTACCGCTTGCGCTTGTCGCTGGGCGAAAGGCCCCGGTGGCAGCAACACGGATCGGGGCAATTTCGGGCTCCTGCTAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGtt                                                                                                   <  2:1008510/150‑1 (MQ=255)
                                                  tAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTtg                                                                                                  >  2:953461/1‑150 (MQ=255)
                                                    gTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTtgct                                                                                                <  1:188455/150‑1 (MQ=255)
                                                     tGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTtgctg                                                                                               >  2:740101/1‑150 (MQ=255)
                                                      gACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTtgctgc                                                                                              >  1:1238512/1‑150 (MQ=255)
                                                      gACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTtgctgc                                                                                              >  1:699203/1‑150 (MQ=255)
                                                      gACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTtgctgc                                                                                              <  1:416174/150‑1 (MQ=255)
                                                      gACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTtgctgc                                                                                              <  1:14002/150‑1 (MQ=255)
                                                       aCCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCa                                                                                             >  2:803501/1‑150 (MQ=255)
                                                        ccGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAg                                                                                            >  1:172804/1‑150 (MQ=255)
                                                          gcttgcgcttgTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGtc                                                                                          >  1:1664072/1‑150 (MQ=255)
                                                          gcttgcgcttgTCGCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGGCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGtc                                                                                          <  1:382080/150‑1 (MQ=255)
                                                          gcttgcgcttgGCTCCGGGTGAAAGGCCCCGGTGGCCGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAATACCCCGGCAGCCTGCTAGCTGTCGGGGGGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGtc                                                                                          >  1:895054/1‑150 (MQ=255)
                                                              gcgcttgTCTCCGGGTGAAATGCCCCGGTTTCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTATCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCg                                                                                      <  2:826065/150‑1 (MQ=255)
                                                              gcgcttgTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCg                                                                                      <  1:1500854/150‑1 (MQ=255)
                                                               cgcttgTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCga                                                                                     <  2:1330117/150‑1 (MQ=255)
                                                                gcttgTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCgag                                                                                    <  1:472693/150‑1 (MQ=255)
                                                                   tgTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGt                                                                                 >  2:1033433/1‑150 (MQ=255)
                                                                      ctcCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGCCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCGCCGACAGTggg                                                                              >  2:177241/1‑149 (MQ=255)
                                                                            gTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAAc                                                                        <  2:567278/150‑1 (MQ=255)
                                                                            gTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAAc                                                                        >  2:754471/1‑150 (MQ=255)
                                                                               aaaGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACttt                                                                     <  2:1649330/150‑1 (MQ=255)
                                                                                   gCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCa                                                                 >  2:1113233/1‑150 (MQ=255)
                                                                                    ccccGGTGGCAGCAACACGGATTGTGGCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGTAAGATTCTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAg                                                                <  1:1054607/150‑1 (MQ=255)
                                                                                     cccGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGc                                                               >  1:587228/1‑150 (MQ=255)
                                                                                      ccGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGACTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCg                                                              <  2:1271609/150‑1 (MQ=255)
                                                                                        ggTGGCAGCAACACGGATTGGTTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCATCCTTCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTATCTCAGCGtt                                                            <  2:1032673/150‑1 (MQ=255)
                                                                                            gcagcaACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGTGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGttctta                                                        >  2:658981/1‑149 (MQ=255)
                                                                                            gcagcaACACGGATTGGGTCATATTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGttcttc                                                        <  2:195823/150‑1 (MQ=255)
                                                                                             cagcaACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGATTTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTcc                                                       >  2:657160/1‑150 (MQ=255)
                                                                                             cagcaACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTcc                                                       <  2:1191641/150‑1 (MQ=255)
                                                                                                 aaCACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTc                                                   >  1:687153/1‑150 (MQ=255)
                                                                                                     cGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCTTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCa                                               >  2:88151/1‑150 (MQ=255)
                                                                                                     cGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCa                                               >  1:927522/1‑150 (MQ=255)
                                                                                                       gATTGGGGCATTTTCTGGCTCATGATAAAATAAAACCCCGGCAGCCTGCTAGAGGTCGGGGTGAATTTGCTGGGGAATATACTACTATCAGTCGTTGCGGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGt                                             <  2:297654/150‑1 (MQ=255)
                                                                                                            ggTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAgg                                        >  2:691484/1‑150 (MQ=255)
                                                                                                                aTTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTccc                                    <  2:927522/150‑1 (MQ=255)
                                                                                                                  tttCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGGTCTTCCGGTCCGCAGTGAAGGTCCCTg                                  <  2:802091/150‑1 (MQ=255)
                                                                                                                      gggCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCttt                              >  1:998238/1‑150 (MQ=255)
                                                                                                                       ggCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTtg                             <  1:740101/150‑1 (MQ=255)
                                                                                                                       ggCTCCTGATAAAAGAAAACCCCGGCAGCCAGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTtg                             <  2:73253/150‑1 (MQ=255)
                                                                                                                        gCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGATGGTCCCTGCTTtgt                            <  1:143842/150‑1 (MQ=255)
                                                                                                                        gCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTtgt                            >  2:35032/1‑150 (MQ=255)
                                                                                                                          tCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGGCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGGCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTgg                          <  2:851671/150‑1 (MQ=255)
                                                                                                                           ccTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCTGTCCGCAGTGAAGGTCCCTGCTTTGTGGt                         <  1:691484/150‑1 (MQ=255)
                                                                                                                            cTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTc                        <  2:1238512/150‑1 (MQ=255)
                                                                                                                             tGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATCTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCAt                      >  2:1531281/1‑150 (MQ=255)
                                                                                                                                tAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCAGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTg                    >  1:1487976/1‑150 (MQ=255)
                                                                                                                                 aaacGAACACCCCGGCAGCCTGCTAGCGGGCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGa                   <  2:267502/150‑1 (MQ=255)
                                                                                                                                 aaaaGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGa                   <  1:803501/150‑1 (MQ=255)
                                                                                                                                 aaaaGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGa                   <  1:953461/150‑1 (MQ=255)
                                                                                                                                  aaaGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGAc                  >  1:488391/1‑150 (MQ=255)
                                                                                                                                   aaGAAAACCCCGGCAGCCTGCTAGCTCTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGCACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGACg                 >  1:105872/1‑150 (MQ=255)
                                                                                                                                    aGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTAGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGACGt                >  2:1349801/1‑150 (MQ=255)
                                                                                                                                     gAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGACGTa               >  1:320344/1‑150 (MQ=255)
                                                                                                                                     gAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTATCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGACGTa               >  1:1679485/1‑150 (MQ=255)
                                                                                                                                        caCCCCGACAGCTAGCAGGCTGCCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCGTTGTGGTCATTGACGTAGAc            >  1:199219/2‑150 (MQ=255)
                                                                                                                                          ccccGACAGCTAGCAGGCTGCCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGACGTAGACGt          >  1:136834/1‑150 (MQ=255)
                                                                                                                                          ccccGACAGCTAGCAGGCTGCCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCAGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGACGTAGACGt          >  1:1209979/1‑150 (MQ=255)
                                                                                                                                           cccGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGACGTAGACGTc         <  2:883055/150‑1 (MQ=255)
                                                                                                                                           cccGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGACGTAGACGTc         >  1:1416364/1‑150 (MQ=255)
                                                                                                                                            ccGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGACGTAGACGTCg        >  2:1313260/1‑150 (MQ=255)
                                                                                                                                            ccGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGACGTAGACGTCg        >  1:112402/1‑150 (MQ=255)
                                                                                                                                            ccGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCCCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGACGTAGACGTCg        >  1:789653/1‑150 (MQ=255)
                                                                                                                                                 aGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACAATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGACGTAGACGTCGAACTg   <  2:14508/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGACGTAGACGTCGAACTGC  >  NC_000913/265430‑265725

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: