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

AGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGACGTAGACGTCGAACTGCCCGG  >  NC_000913/265435‑265729
                                                                                                                                                     |                                                                                                                                                  
aGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAg                                                                                                                                                   >  2:1094358/1‑150 (MQ=255)
aGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAg                                                                                                                                                   >  1:199763/1‑150 (MQ=255)
aGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAg                                                                                                                                                   >  2:855459/1‑150 (MQ=255)
 gTGACCGCGCTGGTCGTCTTCGCGGGGGATATGGCTGGACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGACAGGCCCCGGTGGCAGCAACACGGATTGGGGCAGGTTCGGGCTCCTGATAAAAGAAAACCCCGTCAGCCTGCTAGc                                                                                                                                                  <  1:156866/150‑1 (MQ=255)
     cTGCCCTGGTCGTCTTCGAGTGTGCTATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGtg                                                                                                                                              >  1:828320/1‑149 (MQ=255)
     cTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTc                                                                                                                                              >  2:554233/1‑150 (MQ=255)
       gCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCgg                                                                                                                                            >  2:386853/1‑150 (MQ=255)
        cccTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGACAGCTAGCAGGCTGCCggg                                                                                                                                           <  2:602789/150‑1 (MQ=255)
        cccTGGTCGTCTACGAGGGTGATATTGCTGTACGTGGTAGTGACCCCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGTCAGCCAGACGGATTGGGTCATTTTCGGCCTCCTGATAAAAGCAAACCCCGGCAGCCTGCTAGATGTCggg                                                                                                                                           >  2:72572/1‑150 (MQ=255)
        cccGGGTCGTCTTCGAGTGTGATATTGCGGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGACAGCTAGCAGGCTGCCggg                                                                                                                                           <  1:483203/150‑1 (MQ=255)
           tGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTg                                                                                                                                        <  1:501611/150‑1 (MQ=255)
             gtcgtcTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGaa                                                                                                                                      >  1:752610/1‑150 (MQ=255)
             gtcgtcTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCAGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGaa                                                                                                                                      >  1:175079/1‑150 (MQ=255)
                 tcTTCGAGTGTGATATTGCTGTACGTCGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTg                                                                                                                                  <  2:776965/150‑1 (MQ=255)
                  cTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTgg                                                                                                                                 >  2:15417/1‑149 (MQ=255)
                      gAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTggg                                                                                                                             <  1:17471/150‑1 (MQ=255)
                      gAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTggg                                                                                                                             >  2:1053289/1‑150 (MQ=255)
                      gAGTGTGATATTGCTGTACGTGGTAGGGACCGCTTGCGCTTGTCTCCGGGGGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTggg                                                                                                                             <  2:1100772/150‑1 (MQ=255)
                        gtgtGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTGATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGa                                                                                                                           >  1:174730/1‑150 (MQ=255)
                        gtgtGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGa                                                                                                                           >  2:1105007/1‑150 (MQ=255)
                        gtgtGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGGCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGa                                                                                                                           <  1:1121061/150‑1 (MQ=255)
                          gtgATATTGCTGTCCGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAATCCCCGGCAGCCTGCTGGCTGTCGGGGAGATTTTGCTGGGGacc                                                                                                                         >  1:57113/1‑148 (MQ=255)
                          gtgATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAg                                                                                                                         >  2:1030600/1‑150 (MQ=255)
                          gtgATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGCGTCAATTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGACTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAg                                                                                                                         >  2:859087/1‑150 (MQ=255)
                                ttGCAGTACGTGGTAGTTACCGCTTGCGCTTGTCTCATGGTGAAAGGCCCCGGTGGCAGCAACAAGGATTTGTGCATTTTCTTGATCCTGATAAAAGAAAACCCCTGCAGCATGATAGCTGGCGGGGTGAATTTGCTGGGGAAGAtacta                                                                                                                   <  1:877908/150‑1 (MQ=255)
                                  ccagTACGTGGTAGTGACCGCTTGCGCTTGTCTCAGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATATGCTGGGGAAGAtactaca                                                                                                                 >  2:613498/4‑149 (MQ=255)
                                    tGTACGTGGTAGTGACCACTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTAt                                                                                                               >  1:873328/1‑150 (MQ=255)
                                     gTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATc                                                                                                              >  1:152231/1‑150 (MQ=255)
                                     gTACGTGGGAGAGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGGGGCAGCAACACGGATTGCGGCATTTGCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATc                                                                                                              <  2:118408/150‑1 (MQ=255)
                                       aCGTGGTCCTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAg                                                                                                            <  2:54596/150‑1 (MQ=255)
                                       aCGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTAGGGGTGAATTTGCTGGGGAAGATACAACTATCAg                                                                                                            >  1:649321/1‑150 (MQ=255)
                                          tGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCg                                                                                                         <  1:1064711/150‑1 (MQ=255)
                                              aGTGACCGCTTGCGCTTGGCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTtgc                                                                                                     <  2:646042/150‑1 (MQ=255)
                                               gTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACCCGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTtgct                                                                                                    >  1:417681/1‑150 (MQ=255)
                                                   ccGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCACCACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAAGTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAg                                                                                                >  1:1126671/1‑150 (MQ=255)
                                                        tgcgcttgTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTcc                                                                                           >  2:580227/1‑150 (MQ=255)
                                                         gcgcttgTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCg                                                                                          >  1:993070/1‑150 (MQ=255)
                                                         gcgcttgTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCg                                                                                          <  2:561145/150‑1 (MQ=255)
                                                          cgcttgTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCga                                                                                         >  1:1313725/1‑150 (MQ=255)
                                                           gctttTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCGCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCgag                                                                                        <  2:1179661/150‑1 (MQ=255)
                                                            cttgTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCgag                                                                                        <  1:757461/150‑1 (MQ=255)
                                                              tgTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTGCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGt                                                                                     <  1:858730/150‑1 (MQ=255)
                                                               gTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTg                                                                                    >  1:857589/1‑150 (MQ=255)
                                                                tctcCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTgg                                                                                   <  1:399564/150‑1 (MQ=255)
                                                                 ctcCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCGCCGAGAGTggc                                                                                  >  1:510385/1‑149 (MQ=255)
                                                                     ggggGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTGCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGa                                                                              <  1:155898/150‑1 (MQ=255)
                                                                       gTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAAc                                                                            >  1:1038719/1‑150 (MQ=255)
                                                                        tGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACt                                                                           >  2:1314783/1‑150 (MQ=255)
                                                                         gAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACtt                                                                          >  1:368004/1‑150 (MQ=255)
                                                                          aaaGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACttt                                                                         >  1:606846/1‑150 (MQ=255)
                                                                             ggCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTtctc                                                                      >  2:240965/1‑150 (MQ=255)
                                                                             ggCCCCGGTGGCAGCAACACGGATTGGGGCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTtctc                                                                      <  2:609141/150‑1 (MQ=255)
                                                                               ccccGGTGGCAGCAACACGGTTTGGGGCATTTTGGGGCTCCTGATAAAAGAAAACCCCGGCAGCGTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAg                                                                    <  2:117549/150‑1 (MQ=255)
                                                                                cccGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAgg                                                                   >  1:878723/1‑149 (MQ=255)
                                                                                cccGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGc                                                                   >  2:1237961/1‑150 (MQ=255)
                                                                                 ccGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCg                                                                  <  1:709935/150‑1 (MQ=255)
                                                                                       gcagcaACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAAATTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGttcttc                                                            >  2:1022473/1‑150 (MQ=255)
                                                                                        cagcaACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTTTCAGCGTTCTTcc                                                           <  1:479806/150‑1 (MQ=255)
                                                                                            aaCCCGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTc                                                       <  1:1089389/150‑1 (MQ=255)
                                                                                            aaCACGTCTTGGTTCATTTTCGGGATCCTGATAAAAGAAAACCCCGGCATCCTGCTAGCTGTCGGGGTGAATTTGCTTGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTc                                                       <  2:1154259/150‑1 (MQ=255)
                                                                                            aaCACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGTCAGAACTTTCTCAGCGTTCTTCCGGTc                                                       >  1:550296/1‑150 (MQ=255)
                                                                                            aaCACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACTGAACTTTCTCAGCGTTCTTCCGGTc                                                       >  1:545557/1‑150 (MQ=255)
                                                                                            aaCACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCGGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTc                                                       <  2:1149108/150‑1 (MQ=255)
                                                                                               acGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGc                                                    <  2:606846/150‑1 (MQ=255)
                                                                                                cGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCa                                                   >  1:300512/1‑150 (MQ=255)
                                                                                                 ggATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAg                                                  <  2:718748/150‑1 (MQ=255)
                                                                                                   aTTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCTGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTg                                                <  1:519607/150‑1 (MQ=255)
                                                                                                       gttCATTTTATGGCTCCTGATAAAAGAAAACCCCGGCAGCATGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGGCGTTGCTGCAGTCTCCGAGAGTGGCCAGTACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAgg                                            <  1:162337/148‑1 (MQ=255)
                                                                                                       ggTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCTGTCGCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAgg                                            >  2:162700/1‑150 (MQ=255)
                                                                                                       ggTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAgg                                            >  1:734128/1‑150 (MQ=255)
                                                                                                             tttCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGCCAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTg                                      <  2:927879/150‑1 (MQ=255)
                                                                                                              ttCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGc                                     <  1:706195/150‑1 (MQ=255)
                                                                                                                cGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCGCTGCtt                                   >  2:498878/1‑150 (MQ=255)
                                                                                                                  ggCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTtg                                 >  2:231379/1‑150 (MQ=255)
                                                                                                                     tCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTgg                              <  2:969759/150‑1 (MQ=255)
                                                                                                                      ccTGATAAAAGAAAACCCCGGCAGCCTGCTTGCTGTCGGGGTGAATTTGCTGGGGAAGATTCTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGAAAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGGGTTGTGGt                             >  2:65965/1‑150 (MQ=255)
                                                                                                                      ccTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGt                             <  1:1314783/150‑1 (MQ=255)
                                                                                                                      ccTGATAAAAGAAAACCCCGGCAGCCTGCCAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGt                             >  1:1113963/1‑150 (MQ=255)
                                                                                                                        tGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCa                           >  1:366826/1‑150 (MQ=255)
                                                                                                                         gATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCAt                          >  2:1102081/1‑150 (MQ=255)
                                                                                                                          aTAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCAtt                         <  2:175079/150‑1 (MQ=255)
                                                                                                                           tAATAGAAAACCCCGGCACCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTCCTATCAGTCGGGGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTGCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTg                        <  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)
                                                                                                                                                     |                                                                                                                                                  
AGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCCCGGTGGCAGCAACACGGATTGGGTCATTTTCGGGCTCCTGATAAAAGAAAACCCCGGCAGCCTGCTAGCTGTCGGGGTGAATTTGCTGGGGAAGATACTACTATCAGTCGTTGCTGCAGTCTCCGAGAGTGGACAGAACTTTCTCAGCGTTCTTCCGGTCCGCAGTGAAGGTCCCTGCTTTGTGGTCATTGACGTAGACGTCGAACTGCCCGG  >  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: