Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 292,587 T→C Y121C (TAC→TGC)  yagK ← CP4‑6 prophage; conserved protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913292,5870TC98.2% 422.8 / ‑1.6 111Y121C (TAC→TGC) yagKCP4‑6 prophage; conserved protein
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base C (61/48);  total (61/50)
Fisher's exact test for biased strand distribution p-value = 2.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.13e-01

Only 100 of 121 total aligned reads displayed.

AGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAGGGCAACGGTAAATGCGTTTATCCTCACGCACCTTACGGGTACGGTCAGCCTCCAGCTTA  >  NC_000913/292439‑292733
                                                                                                                                                    |                                                                                                                                                  
 gTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGgcagt                                                                                                                                                 <  2:1315500/151‑1 (MQ=255)
  tACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGgcagt                                                                                                                                                 >  2:317038/1‑150 (MQ=255)
   aCTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGgcagta                                                                                                                                                <  1:1397078/150‑1 (MQ=255)
    cTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGgcagta                                                                                                                                                >  1:569233/1‑149 (MQ=255)
     tATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGc                                                                                                                                             >  2:1244438/1‑151 (MQ=255)
     tATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGc                                                                                                                                             >  2:674100/1‑151 (MQ=255)
       ttGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGCGt                                                                                                                                           >  1:1186447/1‑151 (MQ=255)
        tGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGCGTc                                                                                                                                          >  1:921708/1‑151 (MQ=255)
         ggTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGCGTc                                                                                                                                          <  1:286513/150‑1 (MQ=255)
         ggTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGCGTCt                                                                                                                                         >  2:126406/1‑151 (MQ=255)
         ggTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGCGTCt                                                                                                                                         <  1:461040/151‑1 (MQ=255)
            gtCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCttt                                                                                                                                       <  2:1091280/150‑1 (MQ=255)
              cTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGCGTCTttgtt                                                                                                                                    >  2:1569575/1‑151 (MQ=255)
               tAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGCGTCTttgttg                                                                                                                                   >  1:521177/1‑151 (MQ=255)
                aGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGCGTCTttgtt                                                                                                                                    <  1:338498/149‑1 (MQ=255)
                  cacaTACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGCGTCTTTGTTGaac                                                                                                                                >  1:939605/1‑151 (MQ=255)
                  cacaTACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGCGTCTTTGTTGaac                                                                                                                                <  1:1640271/151‑1 (MQ=255)
                    caTACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGCGTCTTTGTTGaacaa                                                                                                                              <  1:1150952/151‑1 (MQ=255)
                        cTTTCAGTTTTCAGGAAAATGGGCCAAACCCGGGTGATCCTCTCTTCCCAGTCTCAGGGCGCTGTACCCGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGCGTCTTTGTTGAacaaca                                                                                                                            <  2:86540/149‑1 (MQ=255)
                          tACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGCGTCTTTGTTGAACAACAGGc                                                                                                                         >  2:1639028/1‑150 (MQ=255)
                           aCAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGCGTCTTTGTTGAACAACAGGCa                                                                                                                        >  2:298211/1‑150 (MQ=255)
                             aGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGCGTCTTTGTTGAACAACAGGCAAAt                                                                                                                     >  2:41663/1‑151 (MQ=255)
                               ttttCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGCGTCTTTGTTGAACAACAGGCAAAtg                                                                                                                    >  2:960144/1‑150 (MQ=255)
                                   cAGGAAAATGGACCAGCCCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGAt                                                                                                                <  1:795394/150‑1 (MQ=255)
                                   cAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGAt                                                                                                                >  2:1239366/1‑150 (MQ=255)
                                   cAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGAt                                                                                                                >  2:245557/1‑150 (MQ=255)
                                     ggAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGt                                                                                                             >  2:710967/1‑151 (MQ=255)
                                       aaaaTGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTgg                                                                                                           <  1:133386/151‑1 (MQ=255)
                                        aaaTGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGc                                                                                                          >  1:1335072/1‑151 (MQ=255)
                                        aaaTGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCCGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGc                                                                                                          >  1:1515934/1‑151 (MQ=255)
                                         aaTGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCa                                                                                                         <  1:68065/151‑1 (MQ=255)
                                              accagaccagGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTAc                                                                                                    >  1:294820/1‑151 (MQ=255)
                                              accagaccagGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTAc                                                                                                    <  2:819976/151‑1 (MQ=255)
                                                cagaccagGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCt                                                                                                  >  1:417271/1‑151 (MQ=255)
                                                 agaccagGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTg                                                                                                 <  2:1410931/151‑1 (MQ=255)
                                                  gaccagGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGCGGCAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCt                                                                                                  >  1:769784/1‑149 (MQ=255)
                                                   accagGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGAt                                                                                               <  1:204050/151‑1 (MQ=255)
                                                    ccagGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGAtt                                                                                              <  2:69088/151‑1 (MQ=255)
                                                     cagGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTc                                                                                             >  2:778586/1‑151 (MQ=255)
                                                      agGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTc                                                                                             >  2:668463/1‑150 (MQ=255)
                                                      agGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCa                                                                                            <  2:1272390/151‑1 (MQ=255)
                                                        ggTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGa                                                                                          <  1:1141868/151‑1 (MQ=255)
                                                        ggTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGa                                                                                          <  2:817293/151‑1 (MQ=255)
                                                         gTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGa                                                                                          >  1:999210/1‑150 (MQ=255)
                                                            aTCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATAc                                                                                      >  1:336833/1‑151 (MQ=255)
                                                            aTCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATAc                                                                                      <  2:1557526/151‑1 (MQ=255)
                                                               gTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATActct                                                                                   <  1:317012/151‑1 (MQ=255)
                                                               gTCTCTATCCAGTCACAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGAGTCAGAATAcg                                                                                     <  2:1137525/149‑2 (MQ=255)
                                                                tctctctcCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCtt                                                                                  >  1:1051034/1‑151 (MQ=255)
                                                                tctctctcCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCtt                                                                                  >  2:1633033/1‑151 (MQ=255)
                                                                   ctctcCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCttt                                                                                 <  1:299115/149‑1 (MQ=255)
                                                                   ctctcCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCg                                                                               >  2:317502/1‑151 (MQ=255)
                                                                     ctcCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTc                                                                                <  2:393898/148‑1 (MQ=255)
                                                                     ctcCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGcc                                                                             >  1:1512456/1‑151 (MQ=255)
                                                                      tcCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGccc                                                                            >  2:1025945/1‑151 (MQ=255)
                                                                       ccAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGccc                                                                            >  1:739266/1‑150 (MQ=255)
                                                                         aGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCat                                                                          >  2:1159128/1‑150 (MQ=255)
                                                                              cAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATataat                                                                     <  2:1373164/150‑1 (MQ=255)
                                                                                gggCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATaataaa                                                                  <  2:529036/151‑1 (MQ=255)
                                                                                  gcgcTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATaataaata                                                                <  1:1234169/151‑1 (MQ=255)
                                                                                    gcTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAAtaaataaa                                                              <  1:1078153/151‑1 (MQ=255)
                                                                                      tGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAgg                                                            >  2:1294070/1‑151 (MQ=255)
                                                                                       gTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAggg                                                           <  2:925704/151‑1 (MQ=255)
                                                                                       gTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAgg                                                            <  1:1069536/150‑1 (MQ=255)
                                                                                         aCCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAGGGc                                                          <  1:1330147/150‑1 (MQ=255)
                                                                                          ccAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAGGGCaa                                                        <  1:309527/151‑1 (MQ=255)
                                                                                           cAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAGGGCAAc                                                       >  2:1243226/1‑151 (MQ=255)
                                                                                            aGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAGGGCAACg                                                      <  1:1115464/151‑1 (MQ=255)
                                                                                                 ccTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAGGGCAACGGTa                                                   <  2:814773/149‑1 (MQ=255)
                                                                                                  cTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAGGGCAACGGTAAAt                                                >  1:490683/1‑151 (MQ=255)
                                                                                                   tgtgATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAGGGCAACGGTAAATg                                               >  1:913076/1‑151 (MQ=255)
                                                                                                      gATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAGGGCAACGGTAAATGCg                                             <  1:1581456/150‑1 (MQ=255)
                                                                                                      gATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAGGGCAACGGTAAATGCGt                                            <  2:1473386/151‑1 (MQ=255)
                                                                                                       atcatcCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAGGGCAACGGTAAATGCGtt                                           <  1:1333552/151‑1 (MQ=255)
                                                                                                        tcatcCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAGGGCAACGGTAAATGCGttt                                          >  2:1573422/1‑151 (MQ=255)
                                                                                                          atcCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAGGGCAACGGTAAATGCGTTTAt                                        <  2:939675/151‑1 (MQ=255)
                                                                                                            cccTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAGGGCAACGGTAAATGCGTTTATcc                                      <  1:1638923/151‑1 (MQ=255)
                                                                                                             ccTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGCGTCTTTGTTGAACAACAGGCAAAGGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAAAAATTAAAGGGCAACGGTAAATGCGTTTAt                                        >  2:1326965/1‑148 (MQ=255)
                                                                                                               tctAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAGGGCAACGGTAAATGCGTTTATCCTc                                    >  1:263656/1‑150 (MQ=255)
                                                                                                               tctAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAGGGCAACGGTAAATGCGTTTATCCTc                                    >  1:1099367/1‑150 (MQ=255)
                                                                                                                  aaGGTTATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAGGGCAACGGTAAATGCGTTTATCCTCACGc                                >  2:1256109/1‑151 (MQ=255)
                                                                                                                   aGGTTATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAGGGCAACGGTAAATGCGTTTATCCTCACGCa                               <  2:769847/151‑1 (MQ=255)
                                                                                                                   aGGTTATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAGGGCAACGGTAAATGCGTTTATCCTCACGCa                               <  1:1457299/151‑1 (MQ=255)
                                                                                                                     gTTATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAGGGCAACGGTAAATGCGTTTATCCTCACGCAcc                             >  1:368948/1‑151 (MQ=255)
                                                                                                                      ttATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAGGGCAACGGTAAATGCGTTTATCCTCACGCACCt                            >  2:1542857/1‑151 (MQ=255)
                                                                                                                       tATCGTCCTGATCGTAATCCCCGAGGTGGCAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAGGGCAACGGTAAATGCGTTTATCCTCACGCACCtt                           >  2:184660/1‑151 (MQ=255)
                                                                                                                             ccTGATCGTAATCCCCGAGGTGGCAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAGGGCAACGGTAAATGCGTTTATCCTCACGCACCTTACGGGt                     >  2:1707200/1‑151 (MQ=255)
                                                                                                                                gATCGTAATCCCCGAGGTGGCAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAGGGCAACGGTAAATGCGTTTATCCTCACGCACCTTACGGGTACg                  <  1:778522/151‑1 (MQ=255)
                                                                                                                                gATCGTAATCCCCGAGGTGGCAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAGGGCAACGGTAAATGCGTTTATCCTCACGCACCTTACGGGTACg                  <  2:173511/151‑1 (MQ=255)
                                                                                                                                   cGTAATCCCCGAGGTGGCAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAGGGCAACGGTAAATGCGTTTATCCTCACGCACCTTACGGGTACGGt                <  2:1077763/150‑1 (MQ=255)
                                                                                                                                      aaTCCCCGAGGTGGCAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAGGGCAACGGTAAATGCGTTTATCCTCACGCACCTTACGGGTACGGTCAGc            <  2:148578/151‑1 (MQ=255)
                                                                                                                                      aaTCCCCGAGGTGGCAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAGGGCAACGGTAAATGCGTTTATCCTCACGCACCTTACGGGTACGGTCAGc            >  1:349852/1‑151 (MQ=255)
                                                                                                                                         ccccGAGGTGGCAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAGGGCAACGGTAAATGCGTTTATCCTCACGCACCTTACGGGTACGGTCAGCCTc         >  2:1489786/1‑151 (MQ=255)
                                                                                                                                           ccGAGGTGGCAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAGGGCAACGGTAAATGCGTTTATCCTCACGCACCTTACGGGTACGGTCAGCCTCCa       <  2:569284/151‑1 (MQ=255)
                                                                                                                                            cGAGGTGGCAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAGGGCAACGGTAAATGCGTTTATCCTCACGCACCTTACGGGTACGGTCAGCCTcc        >  1:130924/1‑149 (MQ=255)
                                                                                                                                              aGGTGGCAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAGGGCAACGGTAAATGCGTTTATCCTCACGCACCTTACGGGTACGGTCAGCCTCCAGCt    >  1:1143500/1‑151 (MQ=255)
                                                                                                                                               ggtggcAGTAAGCGTTTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAGGGCAACGGTAAATGCGTTTATCCTCACGCACCTTACGGGTACGGTCAGCCTCCAGCtt   <  2:1186536/151‑1 (MQ=255)
                                                                                                                                               ggtggcAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAGGGCAACGGTAAATGCGTTTATCCTCACGCACCTTACGGGTACGGTCAGCCTCCAGCtt   <  2:336859/151‑1 (MQ=255)
                                                                                                                                               ggtggcAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAGGGCAACGGTAAATGCGTTTATCCTCACGCACCTTACGGGTACGGTCAGCCTCCAGCtt   <  2:324410/151‑1 (MQ=255)
                                                                                                                                                gtggtAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAGGGCAACGGTAAATGCGTTTATCCTCACGCACCTTACGGGTACGGTCAGCCTCCAGCTTa  <  2:1548579/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
AGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAGGGCAACGGTAAATGCGTTTATCCTCACGCACCTTACGGGTACGGTCAGCCTCCAGCTTA  >  NC_000913/292439‑292733

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: