Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 286,569 G→A S392N (AGC→AAC)  yagG → CP4‑6 prophage; putative sugar transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913286,5690GA96.9% 1172.0 / 17.7 355S392N (AGC→AAC) yagGCP4‑6 prophage; putative sugar transporter
Reads supporting (aligned to +/- strand):  ref base G (3/8);  new base A (165/178);  total (169/186)
Fisher's exact test for biased strand distribution p-value = 2.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.57e-01

Only 100 of 364 total aligned reads displayed.

TGACGTGGTGGACTACGAGGAGAGCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAGCAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCTCACCGATGCCCGCGTGGAGGCCATCAGCCGGCAGCTGAT  >  NC_000913/286423‑286716
                                                                                                                                                   |                                                                                                                                                   
  aCGTGGTGGACTACGAGGAGAGCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCcaacagcag                                                                                                                                               <  2:1575835/151‑1 (MQ=255)
    gtggtgGACTACGAGGAGAGCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGGGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCg                                                                                                                                             >  2:1457660/1‑151 (MQ=255)
    gtggtgGACTACGAGGAGAGCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCACCAGCAGCg                                                                                                                                             >  6:1534435/1‑151 (MQ=255)
    gtggtgGACTACGAGGAGAGCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCg                                                                                                                                             <  3:363198/151‑1 (MQ=255)
      ggtggACTACGAGGAGAGCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTg                                                                                                                                           <  8:1384203/151‑1 (MQ=255)
          gACTACGAGGAGAGCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGc                                                                                                                                       >  7:532643/1‑151 (MQ=255)
                gaggagAGCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTg                                                                                                                                 <  8:248851/151‑1 (MQ=255)
                 aggagAGCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGa                                                                                                                                <  3:1128737/151‑1 (MQ=255)
                    cgagCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGTCCTGGCGATTGGCGGGGTGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGt                                                                                                                             <  5:1333071/150‑1 (MQ=255)
                     gagCCGCAGCGGTCTCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGTATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTg                                                                                                                            <  6:989134/151‑1 (MQ=255)
                     gagCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGGCGGTTGAGGTg                                                                                                                            >  2:569273/1‑151 (MQ=255)
                     gagCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTg                                                                                                                            <  3:641203/151‑1 (MQ=255)
                       gCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCt                                                                                                                          >  7:117491/1‑151 (MQ=255)
                       gCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCt                                                                                                                          >  6:754243/1‑151 (MQ=255)
                          gCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTcac                                                                                                                       >  5:1004929/1‑151 (MQ=255)
                            aGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTcacc                                                                                                                      >  3:1181759/1‑150 (MQ=255)
                            aGCGGTCGCCGCCTCGACGGGCTGGTGTGCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTcacca                                                                                                                     <  1:636749/151‑1 (MQ=255)
                            aGCGGTCGCCGCCTCGAAGGGCGGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTcacca                                                                                                                     <  6:1663897/151‑1 (MQ=255)
                              gggTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGGGTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAAGTATTCCGCCTACAGCAGCGTGCAGCCGGTTGAGGTGCTCaccacc                                                                                                                   <  7:301236/150‑1 (MQ=255)
                                gTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCAt                                                                                                                 <  7:1231908/151‑1 (MQ=255)
                                  cgccgcCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGGTGAGGTGCTCACCACCATc                                                                                                                <  8:1062860/150‑1 (MQ=255)
                                    ccgccTCGACGTGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTTGGCTGGATCCTTGCGTACGTCAACTATTCCGCCAACATCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCaaa                                                                                                             <  1:481322/151‑1 (MQ=255)
                                    ccgccTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCaaa                                                                                                             <  8:1245655/151‑1 (MQ=255)
                                    ccgccTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCaa                                                                                                              >  6:228811/1‑150 (MQ=255)
                                         tCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAAttct                                                                                                        >  8:387328/1‑151 (MQ=255)
                                            ccGGGCTGGCGGTCTCCACCTACCTGTTCAGCCTGAAGATTCGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAAttctgtt                                                                                                     <  4:1075020/150‑1 (MQ=255)
                                            aCGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAAttctgtt                                                                                                     >  6:694982/1‑151 (MQ=255)
                                            aCGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAAttctgtt                                                                                                     >  4:51871/1‑151 (MQ=255)
                                             cGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAAttctgttc                                                                                                    >  8:1149542/1‑151 (MQ=255)
                                                gCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGc                                                                                                 >  8:1506648/1‑151 (MQ=255)
                                                   ggTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCgt                                                                                               <  7:1359503/150‑1 (MQ=255)
                                                    gtgtTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCgtgg                                                                                             >  5:657206/1‑151 (MQ=255)
                                                     tgtTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCgtggt                                                                                            >  5:1610635/1‑151 (MQ=255)
                                                      gtTCTCCTCCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTAATCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCgtggt                                                                                            <  8:545455/150‑1 (MQ=255)
                                                        tctcCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGcc                                                                                         >  1:1471670/1‑151 (MQ=255)
                                                           ccACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCg                                                                                        <  5:1173744/149‑1 (MQ=255)
                                                                 acctGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGgtggtg                                                                                  <  1:114911/149‑1 (MQ=255)
                                                                  cctGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGctc                                                                               <  3:171631/151‑1 (MQ=255)
                                                                    tGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGctct                                                                              >  4:864109/1‑150 (MQ=255)
                                                                    tGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGctct                                                                              >  7:88671/1‑150 (MQ=255)
                                                                    tGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGctct                                                                              >  7:980368/1‑150 (MQ=255)
                                                                      ttCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTAcg                                                                           >  2:1029406/1‑151 (MQ=255)
                                                                       tCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTAcg                                                                           <  3:1156728/150‑1 (MQ=255)
                                                                           ccTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGGTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGc                                                                      <  4:62764/151‑1 (MQ=255)
                                                                             tGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCAt                                                                    >  7:850398/1‑151 (MQ=255)
                                                                              gAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCa                                                                     <  2:222871/149‑1 (MQ=255)
                                                                                aGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGtt                                                                 <  5:1502500/151‑1 (MQ=255)
                                                                                   ttGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTtcat                                                              >  8:770282/1‑151 (MQ=255)
                                                                                     ggCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTtcatca                                                            <  5:739977/151‑1 (MQ=255)
                                                                                      gCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTcatcat                                                           >  3:1424535/1‑151 (MQ=255)
                                                                                          ggCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTGCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTg                                                       <  6:14025/151‑1 (MQ=255)
                                                                                          ggCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGc                                                         <  2:250924/149‑1 (MQ=255)
                                                                                          ggCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGc                                                         <  5:92822/149‑1 (MQ=255)
                                                                                          ggCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGc                                                         <  5:601678/149‑1 (MQ=255)
                                                                                          ggCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCtt                                                       >  3:101993/1‑150 (MQ=255)
                                                                                            cGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTc                                                     <  2:1525124/151‑1 (MQ=255)
                                                                                             gagtGGCGGGGCGGGGGTGGGCGGGATCCTGGCGTACGACAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAAATCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCg                                                    <  2:81480/148‑1 (MQ=255)
                                                                                             gATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCg                                                    <  1:1221447/151‑1 (MQ=255)
                                                                                                 ggcggggcggTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGctct                                                >  5:193386/1‑151 (MQ=255)
                                                                                                 ggcggggcggTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGctc                                                 <  7:1525094/150‑1 (MQ=255)
                                                                                                    ggggcggTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTAc                                              >  1:1093558/1‑150 (MQ=255)
                                                                                                    ggggcggTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACa                                             <  4:445431/151‑1 (MQ=255)
                                                                                                     gggcggTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACaa                                            <  6:1129179/151‑1 (MQ=255)
                                                                                                      ggcggTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGAGTGGTGCCGGTGGTCCTCTACGAGGGCATGTTCATCATGCTGTCGCTCTACa                                             >  6:639983/1‑149 (MQ=255)
                                                                                                       gcggTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGc                                          <  3:1316889/151‑1 (MQ=255)
                                                                                                        gggtggtGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCt                                         <  2:741033/150‑1 (MQ=255)
                                                                                                        cggTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCt                                         >  3:996157/1‑151 (MQ=255)
                                                                                                         ggtggtGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCTc                                        <  6:491956/151‑1 (MQ=255)
                                                                                                          gtggtgGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCTc                                        >  2:70965/1‑150 (MQ=255)
                                                                                                          gtggtgGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCTc                                        >  5:430545/1‑150 (MQ=255)
                                                                                                          gtggtgGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCTCa                                       <  1:1198129/151‑1 (MQ=255)
                                                                                                          gtggtgGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCTCa                                       <  4:1215633/151‑1 (MQ=255)
                                                                                                           tggtggGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCTCAc                                      >  2:1384275/1‑151 (MQ=255)
                                                                                                            ggtggGCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCTCAcc                                     <  8:608089/151‑1 (MQ=255)
                                                                                                               gggCGGGATCCGGCGCTACGTCAACTCTTGCGCCAAGATCAGCTCGCAGCCGGTGGAGGGTCGCACCACAATGAATATTCTGTGCTGCGTTGTGCCGGTGGTGCTGTACGCGGGCATGTTCATCATGCTGTCGCTCTCCAAGCTCACCGAt                                  <  6:1299030/151‑1 (MQ=255)
                                                                                                                ggCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCTCACCGATg                                 <  4:1501884/151‑1 (MQ=255)
                                                                                                                ggCTGGATCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCTCACCGATg                                 <  4:225075/151‑1 (MQ=255)
                                                                                                                       tCCTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCTCACCGATGCCCGCGt                          >  3:1379037/1‑151 (MQ=255)
                                                                                                                        ccTGGCTTACGTCAATTATTCCGGCAAGAGCAGCGTGCAGCGGGTTGAGGTGGTCACCACGATCATAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCTCACCGATGCCCGCGt                          <  2:296230/150‑1 (MQ=255)
                                                                                                                        ccTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCTCACCGATGCCCGCGTg                         <  3:9314/151‑1 (MQ=255)
                                                                                                                         cTGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCTCACCGATGCCCGCGTgg                        >  6:352147/1‑151 (MQ=255)
                                                                                                                          tGGCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCTCACCGATGCCCGCGTgg                        >  4:1369248/1‑150 (MQ=255)
                                                                                                                           ggCGTACGTCAACTATTGCGCCAACACGAGCGTGCAGCCGGTTGAGGTGCTCACCTCCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCTCACCGATGCCCGCGTgg                        <  1:270449/149‑1 (MQ=255)
                                                                                                                           ggCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCTCACCGATGCCCGCGTgg                        >  4:796558/1‑149 (MQ=255)
                                                                                                                           ggCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCTCACCGATGCCCGCGTGGAg                      <  6:1234682/151‑1 (MQ=255)
                                                                                                                            gCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCTCACCGATGCCCGCGTGGAgg                     >  1:1381940/1‑151 (MQ=255)
                                                                                                                            gCGTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCTCACCGATGCCCGCGTGGAgg                     <  3:158318/151‑1 (MQ=255)
                                                                                                                              gTACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCTCACCGATGCCCGCGTGGAGGcc                   <  5:1070555/151‑1 (MQ=255)
                                                                                                                                tACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCTCACCGATGCCCGCGTGGAGGcc                   >  4:702059/1‑150 (MQ=255)
                                                                                                                                  cttCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCTCACCGATGCCCGCGTGGAGGCCATc                >  8:203092/3‑151 (MQ=255)
                                                                                                                                     cAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCTCACCGATGCCCGCGTGGAGGCCATCAg              <  1:328373/150‑1 (MQ=255)
                                                                                                                                     cAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCTCACCGATGCCCGCGTGGAGGCCATCAGc             <  1:925531/151‑1 (MQ=255)
                                                                                                                                       aCTATTCCGCCAACAGCAGCGTGCAGCCTGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCTCACCGATGCCCGCGTGGAGGCCATCAGCCg           <  2:1204073/151‑1 (MQ=255)
                                                                                                                                       aCTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCTCACCGATGCCCGCGTGGAGGCCATCAGCCg           <  3:566564/151‑1 (MQ=255)
                                                                                                                                         tATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCCCCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCTCACCGATGCCCGCGTGGCGGCAATCAGCCGGc         >  3:84478/1‑151 (MQ=255)
                                                                                                                                           ttCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCTCACCGATGCCCGCGTGGAGGCCATCAGCCGGc         >  4:2614/1‑149 (MQ=255)
                                                                                                                                            tCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCTCACCGATGCCCGCGTGGAGGCCATCAGCCGGCAg       <  3:28766/150‑1 (MQ=255)
                                                                                                                                               gCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCTCACCGATGCCCGCGTGGAGGCCATCAGCCGGCAGCTGa   <  1:1613281/151‑1 (MQ=255)
                                                                                                                                               gCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCTCACCGATGCCCGCGTGGAGGCCATCAGCCGGCAGCTGa   <  3:71232/151‑1 (MQ=255)
                                                                                                                                                ccAACAGCAGCGGGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCTCACCGATGCCCGCGTGGAGGCCATCAGCCGGCAGCTGAt  <  7:764992/151‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TGACGTGGTGGACTACGAGGAGAGCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAGCAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCTCACCGATGCCCGCGTGGAGGCCATCAGCCGGCAGCTGAT  >  NC_000913/286423‑286716

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