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.

TGACGTGGTGGACTACGAGGAGAGCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCG‑TACGTCAACTATTCCGCCAGCAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCTCACCGATGCCCGCGTGGAGGCCATCAGCCGGCAGCTGAT  >  NC_000913/286423‑286716
                                                                                                                                                   |                                                                                                                                                   
  aCGTGGTGGACTACGAGGAGAGCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCG‑TACGTCAACTATTCCGCcaacagcag                                                                                                                                               <  2:1575835/151‑1 (MQ=255)
    gtggtgGACTACGAGGAGAGCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGGGGGGCGGTGGTGGGCTGGATCCTGGCG‑TACGTCAACTATTCCGCCAACAGCAGCg                                                                                                                                             >  2:1457660/1‑151 (MQ=255)
    gtggtgGACTACGAGGAGAGCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCG‑TACGTCAACTATTCCGCCACCAGCAGCg                                                                                                                                             >  6:1534435/1‑151 (MQ=255)
    gtggtgGACTACGAGGAGAGCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCG‑TACGTCAACTATTCCGCCAACAGCAGCg                                                                                                                                             <  3:363198/151‑1 (MQ=255)
      ggtggACTACGAGGAGAGCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCG‑TACGTCAACTATTCCGCCAACAGCAGCGTg                                                                                                                                           <  8:1384203/151‑1 (MQ=255)
          gACTACGAGGAGAGCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCG‑TACGTCAACTATTCCGCCAACAGCAGCGTGCAGc                                                                                                                                       >  7:532643/1‑151 (MQ=255)
                gaggagAGCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCG‑TACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTg                                                                                                                                 <  8:248851/151‑1 (MQ=255)
                 aggagAGCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCG‑TACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGa                                                                                                                                <  3:1128737/151‑1 (MQ=255)
                    cgagCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGTCCTGGCGATTGGCGGGGTGGTGGTGGGCTGGATCCTGGCG‑TACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGt                                                                                                                             <  5:1333071/150‑1 (MQ=255)
                     gagCCGCAGCGGTCTCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGTATCCTGGCG‑TACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTg                                                                                                                            <  6:989134/151‑1 (MQ=255)
                     gagCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCG‑TACGTCAACTATTCCGCCAACAGCAGCGTGCAGGCGGTTGAGGTg                                                                                                                            >  2:569273/1‑151 (MQ=255)
                     gagCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCG‑TACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTg                                                                                                                            <  3:641203/151‑1 (MQ=255)
                       gCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCG‑TACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCt                                                                                                                          >  7:117491/1‑151 (MQ=255)
                       gCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCG‑TACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCt                                                                                                                          >  6:754243/1‑151 (MQ=255)
                          gCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCG‑TACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTcac                                                                                                                       >  5:1004929/1‑151 (MQ=255)
                            aGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCG‑TACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTcacc                                                                                                                      >  3:1181759/1‑150 (MQ=255)
                            aGCGGTCGCCGCCTCGACGGGCTGGTGTGCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCG‑TACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTcacca                                                                                                                     <  1:636749/151‑1 (MQ=255)
                            aGCGGTCGCCGCCTCGAAGGGCGGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCG‑TACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTcacca                                                                                                                     <  6:1663897/151‑1 (MQ=255)
                              gggTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGGGTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCG‑TACGTCAAGTATTCCGCCTACAGCAGCGTGCAGCCGGTTGAGGTGCTCaccacc                                                                                                                   <  7:301236/150‑1 (MQ=255)
                                gTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCG‑TACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCAt                                                                                                                 <  7:1231908/151‑1 (MQ=255)
                                  cgccgcCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCG‑TACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGGTGAGGTGCTCACCACCATc                                                                                                                <  8:1062860/150‑1 (MQ=255)
                                    ccgccTCGACGTGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTTGGCTGGATCCTTGCG‑TACGTCAACTATTCCGCCAACATCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCaaa                                                                                                             <  1:481322/151‑1 (MQ=255)
                                    ccgccTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCG‑TACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCaaa                                                                                                             <  8:1245655/151‑1 (MQ=255)
                                    ccgccTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCG‑TACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCaa                                                                                                              >  6:228811/1‑150 (MQ=255)
                                         tCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCG‑TACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAAttct                                                                                                        >  8:387328/1‑151 (MQ=255)
                                            ccGGGCTGGCGGTCTCCACCTACCTGTTCAGCCTGAAGATTCGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCG‑TACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAAttctgtt                                                                                                     <  4:1075020/150‑1 (MQ=255)
                                            aCGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCG‑TACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAAttctgtt                                                                                                     >  6:694982/1‑151 (MQ=255)
                                            aCGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCG‑TACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAAttctgtt                                                                                                     >  4:51871/1‑151 (MQ=255)
                                             cGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCG‑TACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAAttctgttc                                                                                                    >  8:1149542/1‑151 (MQ=255)
                                                gCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCG‑TACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGc                                                                                                 >  8:1506648/1‑151 (MQ=255)
                                                   ggTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCG‑TACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCgt                                                                                               <  7:1359503/150‑1 (MQ=255)
                                                    gtgtTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCG‑TACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCgtgg                                                                                             >  5:657206/1‑151 (MQ=255)
                                                     tgtTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCG‑TACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCgtggt                                                                                            >  5:1610635/1‑151 (MQ=255)
                                                      gtTCTCCTCCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCG‑TACGTCAACTAATCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCgtggt                                                                                            <  8:545455/150‑1 (MQ=255)
                                                        tctcCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCG‑TACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGcc                                                                                         >  1:1471670/1‑151 (MQ=255)
                                                           ccACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCG‑TACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCg                                                                                        <  5:1173744/149‑1 (MQ=255)
                                                                 acctGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCG‑TACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGgtggtg                                                                                  <  1:114911/149‑1 (MQ=255)
                                                                  cctGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCG‑TACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGctc                                                                               <  3:171631/151‑1 (MQ=255)
                                                                    tGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCG‑TACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGctct                                                                              >  4:864109/1‑150 (MQ=255)
                                                                    tGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCG‑TACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGctct                                                                              >  7:88671/1‑150 (MQ=255)
                                                                    tGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCG‑TACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGctct                                                                              >  7:980368/1‑150 (MQ=255)
                                                                      ttCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCG‑TACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTAcg                                                                           >  2:1029406/1‑151 (MQ=255)
                                                                       tCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCG‑TACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTAcg                                                                           <  3:1156728/150‑1 (MQ=255)
                                                                           ccTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGGTGGATCCTGGCG‑TACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGc                                                                      <  4:62764/151‑1 (MQ=255)
                                                                             tGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCG‑TACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCAt                                                                    >  7:850398/1‑151 (MQ=255)
                                                                              gAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCG‑TACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCa                                                                     <  2:222871/149‑1 (MQ=255)
                                                                                aGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCG‑TACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGtt                                                                 <  5:1502500/151‑1 (MQ=255)
                                                                                   ttGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCG‑TACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTtcat                                                              >  8:770282/1‑151 (MQ=255)
                                                                                     ggCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCG‑TACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTtcatca                                                            <  5:739977/151‑1 (MQ=255)
                                                                                      gCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCG‑TACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTcatcat                                                           >  3:1424535/1‑151 (MQ=255)
                                                                                          ggCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCG‑TACGTCAACTATTGCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTg                                                       <  6:14025/151‑1 (MQ=255)
                                                                                          ggCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCG‑TACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGc                                                         <  2:250924/149‑1 (MQ=255)
                                                                                          ggCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCG‑TACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGc                                                         <  5:92822/149‑1 (MQ=255)
                                                                                          ggCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCG‑TACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGc                                                         <  5:601678/149‑1 (MQ=255)
                                                                                          ggCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCG‑TACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCtt                                                       >  3:101993/1‑150 (MQ=255)
                                                                                            cGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCG‑TACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTc                                                     <  2:1525124/151‑1 (MQ=255)
                                                                                             gagtGGCGGGGCGGGGGTGGGCGGGATCCTGGCG‑TACGACAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAAATCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCg                                                    <  2:81480/148‑1 (MQ=255)
                                                                                             gATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCG‑TACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCg                                                    <  1:1221447/151‑1 (MQ=255)
                                                                                                 ggcggggcggTGGTGGGCTGGATCCTGGCG‑TACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGctct                                                >  5:193386/1‑151 (MQ=255)
                                                                                                 ggcggggcggTGGTGGGCTGGATCCTGGCG‑TACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGctc                                                 <  7:1525094/150‑1 (MQ=255)
                                                                                                    ggggcggTGGTGGGCTGGATCCTGGCG‑TACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTAc                                              >  1:1093558/1‑150 (MQ=255)
                                                                                                    ggggcggTGGTGGGCTGGATCCTGGCG‑TACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACa                                             <  4:445431/151‑1 (MQ=255)
                                                                                                     gggcggTGGTGGGCTGGATCCTGGCG‑TACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACaa                                            <  6:1129179/151‑1 (MQ=255)
                                                                                                      ggcggTGGTGGGCTGGATCCTGGCG‑TACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGAGTGGTGCCGGTGGTCCTCTACGAGGGCATGTTCATCATGCTGTCGCTCTACa                                             >  6:639983/1‑149 (MQ=255)
                                                                                                       gcggTGGTGGGCTGGATCCTGGCG‑TACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGc                                          <  3:1316889/151‑1 (MQ=255)
                                                                                                        gggtggtGGGCTGGATCCTGGCG‑TACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCt                                         <  2:741033/150‑1 (MQ=255)
                                                                                                        cggTGGTGGGCTGGATCCTGGCG‑TACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCt                                         >  3:996157/1‑151 (MQ=255)
                                                                                                         ggtggtGGGCTGGATCCTGGCG‑TACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCTc                                        <  6:491956/151‑1 (MQ=255)
                                                                                                          gtggtgGGCTGGATCCTGGCG‑TACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCTc                                        >  2:70965/1‑150 (MQ=255)
                                                                                                          gtggtgGGCTGGATCCTGGCG‑TACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCTc                                        >  5:430545/1‑150 (MQ=255)
                                                                                                          gtggtgGGCTGGATCCTGGCG‑TACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCTCa                                       <  1:1198129/151‑1 (MQ=255)
                                                                                                          gtggtgGGCTGGATCCTGGCG‑TACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCTCa                                       <  4:1215633/151‑1 (MQ=255)
                                                                                                           tggtggGCTGGATCCTGGCG‑TACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCTCAc                                      >  2:1384275/1‑151 (MQ=255)
                                                                                                            ggtggGCTGGATCCTGGCG‑TACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCTCAcc                                     <  8:608089/151‑1 (MQ=255)
                                                                                                               gggCGGGATCC‑GGCGCTACGTCAACTCTTGCGCCAAGATCAGCTCGCAGCCGGTGGAGGGTCGCACCACAATGAATATTCTGTGCTGCGTTGTGCCGGTGGTGCTGTACGCGGGCATGTTCATCATGCTGTCGCTCTCCAAGCTCACCGAt                                  <  6:1299030/151‑1 (MQ=255)
                                                                                                                ggCTGGATCCTGGCG‑TACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCTCACCGATg                                 <  4:1501884/151‑1 (MQ=255)
                                                                                                                ggCTGGATCCTGGCG‑TACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCTCACCGATg                                 <  4:225075/151‑1 (MQ=255)
                                                                                                                       tCCTGGCG‑TACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCTCACCGATGCCCGCGt                          >  3:1379037/1‑151 (MQ=255)
                                                                                                                        ccTGGCT‑TACGTCAATTATTCCGGCAAGAGCAGCGTGCAGCGGGTTGAGGTGGTCACCACGATCATAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCTCACCGATGCCCGCGt                          <  2:296230/150‑1 (MQ=255)
                                                                                                                        ccTGGCG‑TACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCTCACCGATGCCCGCGTg                         <  3:9314/151‑1 (MQ=255)
                                                                                                                         cTGGCG‑TACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCTCACCGATGCCCGCGTgg                        >  6:352147/1‑151 (MQ=255)
                                                                                                                          tGGCG‑TACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCTCACCGATGCCCGCGTgg                        >  4:1369248/1‑150 (MQ=255)
                                                                                                                           ggCG‑TACGTCAACTATTGCGCCAACACGAGCGTGCAGCCGGTTGAGGTGCTCACCTCCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCTCACCGATGCCCGCGTgg                        <  1:270449/149‑1 (MQ=255)
                                                                                                                           ggCG‑TACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCTCACCGATGCCCGCGTgg                        >  4:796558/1‑149 (MQ=255)
                                                                                                                           ggCG‑TACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCTCACCGATGCCCGCGTGGAg                      <  6:1234682/151‑1 (MQ=255)
                                                                                                                            gCG‑TACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCTCACCGATGCCCGCGTGGAgg                     >  1:1381940/1‑151 (MQ=255)
                                                                                                                            gCG‑TACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCTCACCGATGCCCGCGTGGAgg                     <  3:158318/151‑1 (MQ=255)
                                                                                                                              g‑TACGTCAACTATTCCGCCAACAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCTCACCGATGCCCGCGTGGAGGcc                   <  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)
                                                                                                                                                   |                                                                                                                                                   
TGACGTGGTGGACTACGAGGAGAGCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCG‑TACGTCAACTATTCCGCCAGCAGCAGCGTGCAGCCGGTTGAGGTGCTCACCACCATCAAAATTCTGTTCTGCGTGGTGCCGGTGGTGCTCTACGCGGGCATGTTCATCATGCTGTCGCTCTACAAGCTCACCGATGCCCGCGTGGAGGCCATCAGCCGGCAGCTGAT  >  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: ‑