Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 265,582 T→A 7.6% intergenic (‑39/+38) 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,5820TA 7.6% 337.5 / 12.8 123intergenic (‑39/+38)ykfF/ykfBCP4‑6 prophage; uncharacterized protein/CP4‑6 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (57/56);  new base A (7/3);  total (64/59)
Fisher's exact test for biased strand distribution p-value = 3.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.47e-01

Only 100 of 124 total aligned reads displayed.

GAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGACGTAGACGTCGAACTGCCCGG  >  NC_000913/265433‑265729
                                                                                                                                                     |                                                                                                                                                   
     gACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCGGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTg                                                                                                                                                >  1:735452/1‑150 (MQ=255)
     gACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTAGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTg                                                                                                                                                >  1:1074164/1‑150 (MQ=255)
          cccTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCggg                                                                                                                                           <  1:1384139/150‑1 (MQ=255)
          cccTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGGGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCggg                                                                                                                                           <  1:1713509/150‑1 (MQ=255)
           ccTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGACAGCTAGCAGGCTGCCgggg                                                                                                                                          <  1:1304482/150‑1 (MQ=255)
           ccTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTAGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGACgggg                                                                                                                                          >  2:370103/1‑150 (MQ=255)
             tGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTg                                                                                                                                        <  2:1145924/150‑1 (MQ=255)
             tGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTg                                                                                                                                        <  2:1078076/150‑1 (MQ=255)
             tGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGACAGCTAGCAGGCTGCCGGGGtt                                                                                                                                        <  1:819342/150‑2 (MQ=255)
             cggTCGTCTTCGAGTGTGATATTGCTGTACGTGGAAGTGACCGCTTGCGCTTGGCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGAAAAAAGAAAACCCCGACAGCTAGCAGGCTGCCGGGGtt                                                                                                                                        <  2:771729/149‑2 (MQ=255)
                 cgtcTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAAtt                                                                                                                                    <  1:541090/150‑1 (MQ=255)
                 cggcTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAAtt                                                                                                                                    <  1:975464/147‑1 (MQ=255)
                  gtcTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAACTGTCGGGGTGAAttt                                                                                                                                   >  1:1743929/1‑150 (MQ=255)
                   tcTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTg                                                                                                                                  >  1:1068752/1‑150 (MQ=255)
                    cTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGc                                                                                                                                 >  2:1524383/1‑150 (MQ=255)
                          gtgtGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGa                                                                                                                           >  1:32443/1‑150 (MQ=255)
                          gtgtGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTAGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCAAGCTGTCGGGGTGAATTTGCTGGGGa                                                                                                                           >  2:914510/1‑150 (MQ=255)
                            gtgATATTGCTGTACGTGGTAGTGACCGCTTGCGCATGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAg                                                                                                                         <  1:818499/150‑1 (MQ=255)
                             tgATATTGCTGTACGTGGTAGTGCCCGCTTGCGCTTGTCTCCGTGTGAAATGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCATGATAAAAGAAAACCCCGGCAGCCTGCTAGATGTCGGGGTGAAGTTTCTGGGTAAGa                                                                                                                        >  1:403218/1‑150 (MQ=255)
                             tgATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGTCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGa                                                                                                                        <  1:1675134/150‑1 (MQ=255)
                              gATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGAt                                                                                                                       >  2:1749004/1‑150 (MQ=255)
                               atatTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGAta                                                                                                                      >  1:1522735/1‑150 (MQ=255)
                                  ttGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGAtacta                                                                                                                   >  2:1750199/1‑150 (MQ=255)
                                   tGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGAtactac                                                                                                                  >  2:1482863/1‑150 (MQ=255)
                                    gCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGAtact                                                                                                                    <  2:1140233/147‑1 (MQ=255)
                                    gCTGCACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGAtact                                                                                                                    <  2:909403/147‑1 (MQ=255)
                                    gCTGCACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGAtact                                                                                                                    >  1:909403/1‑147 (MQ=255)
                                     cTGTACGTGGTAGTGACCGCTTGCGCTGGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATCGGGCCAGTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATActacta                                                                                                                <  1:138829/150‑1 (MQ=255)
                                      tGTACTTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGTCCTCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTAt                                                                                                               <  2:174664/150‑1 (MQ=255)
                                       gTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATc                                                                                                              <  1:785353/150‑1 (MQ=255)
                                        tACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCa                                                                                                             <  1:1540458/150‑1 (MQ=255)
                                          cGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGt                                                                                                           >  1:1073260/1‑150 (MQ=255)
                                           gTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTc                                                                                                          >  1:882125/1‑150 (MQ=255)
                                            tGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGATGGGGAAGATACTACTATCAGTCg                                                                                                         >  2:170432/1‑150 (MQ=255)
                                              gTAGTTACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGtt                                                                                                       <  1:394940/150‑1 (MQ=255)
                                              gTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGtt                                                                                                       <  2:1169749/150‑1 (MQ=255)
                                                aGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTtgc                                                                                                     <  1:209643/150‑1 (MQ=255)
                                                aGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTtgc                                                                                                     <  1:1608442/150‑1 (MQ=255)
                                                 gTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGGTAGCTGTCGGGGTGAATTTGCTGGGGAAGAAACTACTATCAGTCGTtgct                                                                                                    >  1:1518160/1‑150 (MQ=255)
                                                 gTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTtgct                                                                                                    >  1:1155970/1‑150 (MQ=255)
                                                   gACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTtgctgc                                                                                                  >  1:507113/1‑150 (MQ=255)
                                                   gACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTATCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTtgctgc                                                                                                  <  2:44073/150‑1 (MQ=255)
                                                     ccGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAg                                                                                                <  2:410653/150‑1 (MQ=255)
                                                          tgcgcttgTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTcc                                                                                           >  1:626913/1‑150 (MQ=255)
                                                             gcttgTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCgag                                                                                        >  2:536819/1‑150 (MQ=255)
                                                                  tctcCGGGCGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTgg                                                                                   <  1:1207861/150‑1 (MQ=255)
                                                                   ctcCGGGTGAAAGGCCCCGTTAGCAGCAACACGGATTGGTTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGa                                                                                  <  1:50883/150‑1 (MQ=255)
                                                                     ccGGGGGAACGGCCCCGGTGGCCGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCGGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACa                                                                                <  1:1381134/150‑1 (MQ=255)
                                                                       gggTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGa                                                                              >  2:1714973/1‑150 (MQ=255)
                                                                        ggTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCTGTCTCCGAGAGTGGACAGaa                                                                             >  2:1797601/1‑150 (MQ=255)
                                                                          tGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACt                                                                           >  1:1224054/1‑150 (MQ=255)
                                                                                gCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCa                                                                     <  2:230260/150‑1 (MQ=255)
                                                                                 ccccGGTGGCAGCAACACGGATTGGGTCATTTTCTGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGCGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAg                                                                    >  2:1398305/1‑150 (MQ=255)
                                                                                 ccccGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAg                                                                    <  1:1595195/150‑1 (MQ=255)
                                                                                   ccGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCg                                                                  <  2:134235/150‑1 (MQ=255)
                                                                                      gTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGttc                                                               <  2:612438/150‑1 (MQ=255)
                                                                                         gcagcaACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGttcttc                                                            <  2:624123/150‑1 (MQ=255)
                                                                                          cagcaACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGGTGCAGTCTCCGAGAGTGGCCAGAACTTTCTCAGCGTTCTTcc                                                           >  2:1246424/1‑150 (MQ=255)
                                                                                               acacGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTAAGCGTTCTTCCGGTcc                                                      >  2:1383692/1‑150 (MQ=255)
                                                                                                      ttGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTTGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGa                                               <  1:291218/150‑1 (MQ=255)
                                                                                                         ggTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAgg                                            >  2:990861/1‑150 (MQ=255)
                                                                                                         ggTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAgg                                            >  1:925525/1‑150 (MQ=255)
                                                                                                         ggTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAgg                                            <  2:1053610/150‑1 (MQ=255)
                                                                                                              ttttCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTcccg                                       >  2:1625102/1‑149 (MQ=255)
                                                                                                               tttCGGGCTCCTGATAACAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGGCCCCGAGAGTGGACAGAACTTTCTCAGCGTGCTTCCGGTCCGCAGTGAAGGTCCCTg                                      <  2:245998/150‑1 (MQ=255)
                                                                                                                 tCGGGCTCCTGATACAAGACAACACCGGCAGCCTGCTAGCTGGCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCt                                    <  1:661091/150‑1 (MQ=255)
                                                                                                                       tCCTGATAAAAGAAAAGCGCGGCATCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACATAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTgg                              <  2:1225240/150‑1 (MQ=255)
                                                                                                                       tCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTgg                              <  2:1219043/150‑1 (MQ=255)
                                                                                                                        ccTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGt                             <  1:1000496/150‑1 (MQ=255)
                                                                                                                           gATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCAt                          >  1:1897687/1‑150 (MQ=255)
                                                                                                                           gATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCAt                          >  1:481825/1‑150 (MQ=255)
                                                                                                                           gATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCTGCTTTGTGGTCAtt                         <  2:1122471/150‑1 (MQ=255)
                                                                                                                             tAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGGCCCTGCTTTGTGGTCATTg                        >  2:216402/1‑150 (MQ=255)
                                                                                                                              aaaaGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGa                       <  2:1698646/150‑1 (MQ=255)
                                                                                                                               aaaGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGAc                      >  2:556591/1‑150 (MQ=255)
                                                                                                                                aaGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGACg                     <  1:602306/150‑1 (MQ=255)
                                                                                                                                 aGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGACGt                    <  2:620747/150‑1 (MQ=255)
                                                                                                                                  caaaaCCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGACGTa                   <  1:804187/149‑1 (MQ=255)
                                                                                                                                   aaaaCCCGGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGTGAGTGGACAGAACTTTCTCAGCGTTCTTCGGGTCCGCTGTGAAGGTCCCTGCTTTGTGGTCATTGACGTAg                  <  2:1518160/150‑1 (MQ=255)
                                                                                                                                   aaaaCCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGGCAGAACTTTCTCAGCGGTCTTCCGGTCCGCTGTGAAGGTCCCTGCTTTGTGGTCATTGACGTat                  >  2:942138/1‑149 (MQ=255)
                                                                                                                                   aaaaCCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGACGTAg                  <  2:1315441/150‑1 (MQ=255)
                                                                                                                                   aaaaCCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGACGTAg                  <  1:536819/150‑1 (MQ=255)
                                                                                                                                    tcaCCCCGACAGCTAGCAGGCTGCCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGACGTAGa                 >  1:695634/3‑150 (MQ=255)
                                                                                                                                    tcaCCCCGACAGCTAGCAGGCTGCCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGACGTAGa                 >  1:702291/3‑150 (MQ=255)
                                                                                                                                    aaaCCCCGGCAGCCTGCTAGCGGGCGGGGTGAATGTGCTGGGGAAGAAACTACTATCAGGCGTTGCGGCAGTCGCCGAGAGTGGACAGAACTTTCTCAGCGTTCATCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGACGTAGa                 <  2:225301/150‑1 (MQ=255)
                                                                                                                                     caCCCCGACAGCTAGCAGGCTGCCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGACGTAGAc                >  2:58910/2‑150 (MQ=255)
                                                                                                                                      aCCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGACGTAGACg               >  1:435373/1‑150 (MQ=255)
                                                                                                                                       ccccGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATAATACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGACGTAGACGt              >  2:1257139/1‑150 (MQ=255)
                                                                                                                                       ccccGACAGCTAGCAGGCTGCCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGACGTAGACGt              >  1:488240/1‑150 (MQ=255)
                                                                                                                                       ccccGACAGCTAGCAGGCTGCCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGACGTAGACGt              >  2:1756872/1‑150 (MQ=255)
                                                                                                                                        cccGTCAGCCTGCTAGCTGGCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGGCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGACGTAGACGTc             <  1:914510/150‑1 (MQ=255)
                                                                                                                                         ccGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGACGTAGACGTCg            <  2:1019647/150‑1 (MQ=255)
                                                                                                                                         ccGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGACGTAGACGTCg            >  1:156399/1‑150 (MQ=255)
                                                                                                                                            gCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGACGTAGACGTCGAAc         <  2:735452/150‑1 (MQ=255)
                                                                                                                                             cAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGACGTAGACGTCGAACt        >  2:1096846/1‑150 (MQ=255)
                                                                                                                                              aGCCTGCTAGCTGTGGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGAGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGACGTCTACGTcaatgac       <  2:323291/150‑7 (MQ=255)
                                                                                                                                               gCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGACGTAGACGTCGAACTGc      >  1:1138/1‑150 (MQ=255)
                                                                                                                                                 cTGCTAGCTGGCGGGGTGAATTTGCTGGGGACGACACTACTATCAGTCGATGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGACGTAGACGTCGAACTGccc    <  1:193819/150‑1 (MQ=255)
                                                                                                                                                   gCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGACGTAGACGTCGAACTGCCCgg  >  1:1420530/1‑150 (MQ=255)
                                                                                                                                                   gCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGGCGTTGCTGCAGGCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGACGTAGACGTCGAACTGCCCgg  <  2:1708673/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGACGTAGACGTCGAACTGCCCGG  >  NC_000913/265433‑265729

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: