Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 265,583 A→G 5.4% intergenic (‑40/+37) 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,5830AG 5.4% 277.3 / 3.8 93intergenic (‑40/+37)ykfF/ykfBCP4‑6 prophage; uncharacterized protein/CP4‑6 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (49/38);  new base G (3/2);  total (53/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGAT‑AAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGACGTAGACGTCGAACTGCCCGG  >  NC_000913/265435‑265729
                                                                                                                                                     |                                                                                                                                                  
aGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGAT‑AAAAGAAAACCCCGGCAGCCTGCTAg                                                                                                                                                   >  2:1094358/1‑150 (MQ=255)
aGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGAT‑AAAAGAAAACCCCGGCAGCCTGCTAg                                                                                                                                                   >  1:199763/1‑150 (MQ=255)
aGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGAT‑AAAAGAAAACCCCGGCAGCCTGCTAg                                                                                                                                                   >  2:855459/1‑150 (MQ=255)
 gTGACCGCGCTGGTCGTCTTCGCGGGGGATATGGCTGGACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGACAGGCCCCGGTGGCAGCAACACGGATTGGGGCAGGTTCGGGCTCCTGAT‑AAAAGAAAACCCCGTCAGCCTGCTAGc                                                                                                                                                  <  1:156866/150‑1 (MQ=255)
     cTGCCCTGGTCGTCTTCGAGTGTGCTATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGAT‑AAAAGAAAACCCCGGCAGCCTGCTAGCTGtg                                                                                                                                              >  1:828320/1‑149 (MQ=255)
     cTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGAT‑AAAAGAAAACCCCGGCAGCCTGCTAGCTGTc                                                                                                                                              >  2:554233/1‑150 (MQ=255)
       gCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGAT‑AAAAGAAAACCCCGGCAGCCTGCTAGCTGTCgg                                                                                                                                            >  2:386853/1‑150 (MQ=255)
        cccTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGAT‑AAAAGAAAACCCCGACAGCTAGCAGGCTGCCggg                                                                                                                                           <  2:602789/150‑1 (MQ=255)
        cccTGGTCGTCTACGAGGGTGATATTGCTGTACGTGGTAGTGACCCCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGTCAGCCAGACGGATTGGGTCATTTTCGGCCTCCTGAT‑AAAAGCAAACCCCGGCAGCCTGCTAGATGTCggg                                                                                                                                           >  2:72572/1‑150 (MQ=255)
        cccGGGTCGTCTTCGAGTGTGATATTGCGGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGAT‑AAAAGAAAACCCCGACAGCTAGCAGGCTGCCggg                                                                                                                                           <  1:483203/150‑1 (MQ=255)
           tGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGAT‑AAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTg                                                                                                                                        <  1:501611/150‑1 (MQ=255)
             gtcgtcTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGAT‑AAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGaa                                                                                                                                      >  1:752610/1‑150 (MQ=255)
             gtcgtcTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCAGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGAT‑AAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGaa                                                                                                                                      >  1:175079/1‑150 (MQ=255)
                 tcTTCGAGTGTGATATTGCTGTACGTCGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGAT‑AAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTg                                                                                                                                  <  2:776965/150‑1 (MQ=255)
                  cTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGAT‑AAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTgg                                                                                                                                 >  2:15417/1‑149 (MQ=255)
                      gAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGAT‑AAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTggg                                                                                                                             <  1:17471/150‑1 (MQ=255)
                      gAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGAT‑AAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTggg                                                                                                                             >  2:1053289/1‑150 (MQ=255)
                      gAGTGTGATATTGCTGTACGTGGTAGGGACCGCTTGCGCTTGTCTCCGGGGGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGAT‑AAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTggg                                                                                                                             <  2:1100772/150‑1 (MQ=255)
                        gtgtGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTGATTTTCGGGCTCCTGAT‑AAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGa                                                                                                                           >  1:174730/1‑150 (MQ=255)
                        gtgtGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGAT‑AAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGa                                                                                                                           >  2:1105007/1‑150 (MQ=255)
                        gtgtGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGGCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGAT‑AAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGa                                                                                                                           <  1:1121061/150‑1 (MQ=255)
                          gtgATATTGCTGTCCGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGAT‑AAAAGAAATCCCCGGCAGCCTGCTGGCTGTCGGGGAGATTTTGCTGGGGacc                                                                                                                         >  1:57113/1‑148 (MQ=255)
                          gtgATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGAT‑AAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAg                                                                                                                         >  2:1030600/1‑150 (MQ=255)
                          gtgATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGCGTCAATTTCGGGCTCCTGAT‑AAAAGAAAACCCCGGCAGACTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAg                                                                                                                         >  2:859087/1‑150 (MQ=255)
                                ttGCAGTACGTGGTAGTTACCGCTTGCGCTTGTCTCATGGTGAAAGGCCCCGGTGGCAGCAACAAGGATTTGTGCATTTTCTTGATCCTGAT‑AAAAGAAAACCCCTGCAGCATGATAGCTGGCGGGGTGAATTTGCTGGGGAAGAtacta                                                                                                                   <  1:877908/150‑1 (MQ=255)
                                  ccagTACGTGGTAGTGACCGCTTGCGCTTGTCTCAGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGAT‑AAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATATGCTGGGGAAGAtactaca                                                                                                                 >  2:613498/4‑149 (MQ=255)
                                    tGTACGTGGTAGTGACCACTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGAT‑AAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTAt                                                                                                               >  1:873328/1‑150 (MQ=255)
                                     gTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGAT‑AAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATc                                                                                                              >  1:152231/1‑150 (MQ=255)
                                     gTACGTGGGAGAGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGGGGCAGCAACACGGATTGCGGCATTTGCGGGCTCCTGAT‑AAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATc                                                                                                              <  2:118408/150‑1 (MQ=255)
                                       aCGTGGTCCTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGAT‑AAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAg                                                                                                            <  2:54596/150‑1 (MQ=255)
                                       aCGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGAT‑AAAAGAAAACCCCGGCAGCCTGCTAGCTGTAGGGGTGAATTTGCTGGGGAAGATACAACTATCAg                                                                                                            >  1:649321/1‑150 (MQ=255)
                                          tGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGAT‑AAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCg                                                                                                         <  1:1064711/150‑1 (MQ=255)
                                              aGTGACCGCTTGCGCTTGGCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGAT‑AAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTtgc                                                                                                     <  2:646042/150‑1 (MQ=255)
                                               gTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACCCGGATTGGGTCATTTTCGGGCTCCTGAT‑AAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTtgct                                                                                                    >  1:417681/1‑150 (MQ=255)
                                                   ccGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCACCACGGATTGGGTCATTTTCGGGCTCCTGAT‑AAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAAGTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAg                                                                                                >  1:1126671/1‑150 (MQ=255)
                                                        tgcgcttgTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGAT‑AAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTcc                                                                                           >  2:580227/1‑150 (MQ=255)
                                                         gcgcttgTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGAT‑AAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCg                                                                                          >  1:993070/1‑150 (MQ=255)
                                                         gcgcttgTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGAT‑AAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCg                                                                                          <  2:561145/150‑1 (MQ=255)
                                                          cgcttgTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGAT‑AAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCga                                                                                         >  1:1313725/1‑150 (MQ=255)
                                                           gctttTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCGCCTGAT‑AAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCgag                                                                                        <  2:1179661/150‑1 (MQ=255)
                                                            cttgTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCgag                                                                                        <  1:757461/150‑1 (MQ=255)
                                                              tgTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTGCGGGCTCCTGAT‑AAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGt                                                                                     <  1:858730/150‑1 (MQ=255)
                                                               gTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGAT‑AAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTg                                                                                    >  1:857589/1‑150 (MQ=255)
                                                                tctcCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGAT‑AAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTgg                                                                                   <  1:399564/150‑1 (MQ=255)
                                                                 ctcCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGAT‑AAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCGCCGAGAGTggc                                                                                  >  1:510385/1‑149 (MQ=255)
                                                                     ggggGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTGCGGGCTCCTGAT‑AAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGa                                                                              <  1:155898/150‑1 (MQ=255)
                                                                       gTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGAT‑AAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAAc                                                                            >  1:1038719/1‑150 (MQ=255)
                                                                        tGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGAT‑AAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACt                                                                           >  2:1314783/1‑150 (MQ=255)
                                                                         gAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGAT‑AAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACtt                                                                          >  1:368004/1‑150 (MQ=255)
                                                                          aaaGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGAT‑AAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACttt                                                                         >  1:606846/1‑150 (MQ=255)
                                                                             ggCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGAT‑AAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTtctc                                                                      >  2:240965/1‑150 (MQ=255)
                                                                             ggCCCCGGTGGCAGCAACACGGATTGGGGCATTTTCGGGCTCCTGAT‑AAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTtctc                                                                      <  2:609141/150‑1 (MQ=255)
                                                                               ccccGGTGGCAGCAACACGGTTTGGGGCATTTTGGGGCTCCTGAT‑AAAAGAAAACCCCGGCAGCGTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAg                                                                    <  2:117549/150‑1 (MQ=255)
                                                                                cccGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGAT‑AAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAgg                                                                   >  1:878723/1‑149 (MQ=255)
                                                                                cccGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGAT‑AAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGc                                                                   >  2:1237961/1‑150 (MQ=255)
                                                                                 ccGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGAT‑AAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCg                                                                  <  1:709935/150‑1 (MQ=255)
                                                                                       gcagcaACACGGATTGGGTCATTTTCGGGCTCCTGAT‑AAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAAATTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGttcttc                                                            >  2:1022473/1‑150 (MQ=255)
                                                                                        cagcaACACGGATTGGGTCATTTTCGGGCTCCTGAT‑AAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTTTCAGCGTTCTTcc                                                           <  1:479806/150‑1 (MQ=255)
                                                                                            aaCCCGGATTGGGTCATTTTCGGGCTCCTGAT‑AAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTc                                                       <  1:1089389/150‑1 (MQ=255)
                                                                                            aaCACGTCTTGGTTCATTTTCGGGATCCTGAT‑AAAAGAAAACCCCGGCATCCTGCTAGCTGTCGGGGTGAATTTGCTTGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTc                                                       <  2:1154259/150‑1 (MQ=255)
                                                                                            aaCACGGATTGGGTCATTTTCGGGCTCCTGAT‑AAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGTCAGAACTTTCTCAGCGTTCTTCCGGTc                                                       >  1:550296/1‑150 (MQ=255)
                                                                                            aaCACGGATTGGGTCATTTTCGGGCTCCTGAT‑AAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACTGAACTTTCTCAGCGTTCTTCCGGTc                                                       >  1:545557/1‑150 (MQ=255)
                                                                                            aaCACGGATTGGGTCATTTTCGGGCTCCTGAT‑AAAAGAAAACCCCGGCAGCCGGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTc                                                       <  2:1149108/150‑1 (MQ=255)
                                                                                               acGGATTGGGTCATTTTCGGGCTCCTGAT‑AAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGc                                                    <  2:606846/150‑1 (MQ=255)
                                                                                                cGGATTGGGTCATTTTCGGGCTCCTGAT‑AAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCa                                                   >  1:300512/1‑150 (MQ=255)
                                                                                                 ggATTGGGTCATTTTCGGGCTCCTGAT‑AAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAg                                                  <  2:718748/150‑1 (MQ=255)
                                                                                                   aTTGGGTCATTTTCGGGCTCCTGAT‑AAAAGAAAACCCCTGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTg                                                <  1:519607/150‑1 (MQ=255)
                                                                                                       gttCATTTTATGGCTCCTGAT‑AAAAGAAAACCCCGGCAGCATGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGGCGTTGCTGCAGTCTCCGAGAGTGGCCAGTACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAgg                                            <  1:162337/148‑1 (MQ=255)
                                                                                                       ggTCATTTTCGGGCTCCTGAT‑AAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCTGTCGCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAgg                                            >  2:162700/1‑150 (MQ=255)
                                                                                                       ggTCATTTTCGGGCTCCTGAT‑AAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAgg                                            >  1:734128/1‑150 (MQ=255)
                                                                                                             tttCGGGCTCCTGAT‑AAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGCCAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTg                                      <  2:927879/150‑1 (MQ=255)
                                                                                                              ttCGGGCTCCTGAT‑AAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGc                                     <  1:706195/150‑1 (MQ=255)
                                                                                                                cGGGCTCCTGAT‑AAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCGCTGCtt                                   >  2:498878/1‑150 (MQ=255)
                                                                                                                  ggCTCCTGAT‑AAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTtg                                 >  2:231379/1‑150 (MQ=255)
                                                                                                                     tCCTGAT‑AAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTgg                              <  2:969759/150‑1 (MQ=255)
                                                                                                                      ccTGAT‑AAAAGAAAACCCCGGCAGCCTGCTTGCTGTCGGGGTGAATTTGCTGGGGAAGATTCTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGAAAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGGGTTGTGGt                             >  2:65965/1‑150 (MQ=255)
                                                                                                                      ccTGAT‑AAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGt                             <  1:1314783/150‑1 (MQ=255)
                                                                                                                      ccTGAT‑AAAAGAAAACCCCGGCAGCCTGCCAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGt                             >  1:1113963/1‑150 (MQ=255)
                                                                                                                        tGAT‑AAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCa                           >  1:366826/1‑150 (MQ=255)
                                                                                                                         gAT‑AAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCAt                          >  2:1102081/1‑150 (MQ=255)
                                                                                                                          aT‑AAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCAtt                         <  2:175079/150‑1 (MQ=255)
                                                                                                                           t‑AATAGAAAACCCCGGCACCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTCCTATCAGTCGGGGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTGCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTg                        <  2:649321/150‑1 (MQ=255)
                                                                                                                             aaaaGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGa                       <  1:377886/150‑1 (MQ=255)
                                                                                                                             aaaaGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATAAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTACGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGa                       >  1:602197/1‑150 (MQ=255)
                                                                                                                               aaGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGACg                     >  1:265617/1‑150 (MQ=255)
                                                                                                                                   aaaCCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGACGTAGa                 >  2:664675/1‑150 (MQ=255)
                                                                                                                                     aCCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGACGTAGACg               <  2:734128/150‑1 (MQ=255)
                                                                                                                                      ccccGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGACGTAGACGt              >  2:409549/1‑150 (MQ=255)
                                                                                                                                      ccccGACAGCTAGCAGGCTGCCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACCGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGACGTAGACGt              >  1:1019533/1‑150 (MQ=255)
                                                                                                                                      ccccGACAGCTAGCAGGCTGCCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGACGTAGACGt              >  2:778623/1‑150 (MQ=255)
                                                                                                                                        ccGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGACGTAGCCGTCg            >  1:1038255/1‑150 (MQ=255)
                                                                                                                                         cGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGACGTAGACGTCGa           <  1:302120/150‑1 (MQ=255)
                                                                                                                                                  gCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGACGTAGACGTCGAACTGCCCgg  <  1:1216937/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
AGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGAT‑AAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGACGTAGACGTCGAACTGCCCGG  >  NC_000913/265435‑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: ‑