Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,298,199 C→T A390T (GCA→ACA)  fdhF ← formate dehydrogenase‑H, selenopolypeptide subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,298,1990CT97.8% 1378.9 / 11.1 407A390T (GCA→ACA) fdhFformate dehydrogenase‑H, selenopolypeptide subunit
Reads supporting (aligned to +/- strand):  ref base C (5/4);  new base T (193/203);  total (199/208)
Fisher's exact test for biased strand distribution p-value = 7.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.85e-01

Only 100 of 429 total aligned reads displayed.

GCCGCCGACGCGGTTTTGGTCATAAAGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGCGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGC  >  NC_000913/4298051‑4298346
                                                                                                                                                    |                                                                                                                                                   
gccgccgACGCGGTTTTGGTCATAAAGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCAtgcgctgcgc                                                                                                                                                   >  8:167607/1‑151 (MQ=255)
 ccgccgACGCGGTTTTGGTCATAAAGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGtgcg                                                                                                                                                  >  2:1429846/1‑151 (MQ=255)
  cgccgACGCGGTTTTGGTCATAAAGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCgg                                                                                                                                                 >  1:1184561/1‑151 (MQ=255)
   gccgACGCGGTTTTGGTCATAAAGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGt                                                                                                                                                >  4:861580/1‑151 (MQ=255)
     cgACGCGGTTTTGGTCATAAAGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGc                                                                                                                                              >  7:1523330/1‑151 (MQ=255)
       aCGCGGTTTTGGTCATAAAGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGtgcggtgcgg                                                                                                                                            <  4:430642/151‑1 (MQ=255)
          cgGTTTTGGTCATAAAGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGc                                                                                                                                           >  6:350117/1‑149 (MQ=255)
           ggTTTTGGTCATAAAGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGc                                                                                                                                        <  2:430832/151‑1 (MQ=255)
              tttGGTCATAAAGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTc                                                                                                                                      <  8:632069/150‑1 (MQ=255)
                  gTCATAAAGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTCCTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGa                                                                                                                                 <  8:848942/151‑1 (MQ=255)
                     aTAAAGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCtgctgc                                                                                                                              >  4:338378/1‑151 (MQ=255)
                     aTAAAGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGc                                                                                                                              <  5:1193319/151‑1 (MQ=255)
                       aaaGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGa                                                                                                                            <  6:1000030/151‑1 (MQ=255)
                        aaGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGAt                                                                                                                           <  6:1667699/151‑1 (MQ=255)
                         aGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATa                                                                                                                          <  1:543067/151‑1 (MQ=255)
                             gTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCg                                                                                                                      <  8:1389955/151‑1 (MQ=255)
                               ccTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGt                                                                                                                    <  3:1548357/151‑1 (MQ=255)
                               ccTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGt                                                                                                                    <  4:1525808/151‑1 (MQ=255)
                                  gAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATg                                                                                                                 <  7:1048333/151‑1 (MQ=255)
                                  gAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATg                                                                                                                 <  7:883538/151‑1 (MQ=255)
                                       aTAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTg                                                                                                            >  5:1316931/1‑151 (MQ=255)
                                       aTAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTg                                                                                                            <  7:255916/151‑1 (MQ=255)
                                        tAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTgg                                                                                                           >  5:849655/1‑151 (MQ=255)
                                             aGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGc                                                                                                      >  7:728486/1‑151 (MQ=255)
                                             aGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGc                                                                                                      >  3:1423939/1‑151 (MQ=255)
                                              gTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCt                                                                                                     >  2:1396093/1‑151 (MQ=255)
                                              gTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCt                                                                                                     >  6:94704/1‑151 (MQ=255)
                                                 ccAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTc                                                                                                  >  8:752367/1‑151 (MQ=255)
                                                   aGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTcc                                                                                                 >  6:1144311/1‑150 (MQ=255)
                                                   aGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTcc                                                                                                 >  2:254606/1‑150 (MQ=255)
                                                         tCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCAc                                                                                               <  6:861452/146‑1 (MQ=255)
                                                          cAACGGCTTTACGTACTGCGGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCa                                                                                         <  6:484291/151‑1 (MQ=255)
                                                          cAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGcccc                                                                                          >  1:101016/1‑150 (MQ=255)
                                                          cAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGcccc                                                                                          <  3:334085/150‑1 (MQ=255)
                                                          cAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCa                                                                                         <  4:1291691/151‑1 (MQ=255)
                                                           aaaGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAg                                                                                        <  5:428350/151‑1 (MQ=255)
                                                              ggCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCt                                                                                     <  1:1427168/151‑1 (MQ=255)
                                                               gCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCtt                                                                                    <  4:957475/151‑1 (MQ=255)
                                                                 tttACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCGCAGGCTTTg                                                                                  >  5:1381839/1‑151 (MQ=255)
                                                                 tttACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTg                                                                                  >  4:1676158/1‑151 (MQ=255)
                                                                  ttccGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCttt                                                                                   <  2:853411/146‑1 (MQ=255)
                                                                    aCGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCg                                                                               <  7:605951/151‑1 (MQ=255)
                                                                    aCGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCg                                                                               <  2:904049/151‑1 (MQ=255)
                                                                       tACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAAt                                                                            >  3:926779/1‑151 (MQ=255)
                                                                       tACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAAt                                                                            >  8:687437/1‑151 (MQ=255)
                                                                       tACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAAt                                                                            >  1:312529/1‑151 (MQ=255)
                                                                        aCTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCGATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAAtt                                                                           <  4:1083784/151‑1 (MQ=255)
                                                                        aCTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAAtt                                                                           <  1:1501762/151‑1 (MQ=255)
                                                                        aCTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAAtt                                                                           >  2:102727/1‑151 (MQ=255)
                                                                         cTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAAttt                                                                          >  8:1523347/1‑151 (MQ=255)
                                                                           gCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAAttt                                                                          <  6:582971/149‑1 (MQ=255)
                                                                           gCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTtct                                                                        <  2:1539894/151‑1 (MQ=255)
                                                                           gCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTtct                                                                        >  3:1020732/1‑151 (MQ=255)
                                                                            ccGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTtctc                                                                       <  5:295233/151‑1 (MQ=255)
                                                                              gACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTcgc                                                                     <  4:657092/151‑1 (MQ=255)
                                                                               aCAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTcgcg                                                                    >  1:1575278/1‑151 (MQ=255)
                                                                                   gtcggtgTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGTCTTTGGCGAATTTCTCGCGGTTa                                                                <  1:1101524/145‑1 (MQ=255)
                                                                                     ccGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTa                                                                <  5:1479351/149‑1 (MQ=255)
                                                                                      cgcgTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGcc                                                             <  1:851743/151‑1 (MQ=255)
                                                                                           cAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATc                                                        <  2:1345457/151‑1 (MQ=255)
                                                                                           cAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATc                                                        >  4:1203691/1‑151 (MQ=255)
                                                                                           cAGTTTGTAGCGGATCTTCGCCCATAATGTACGCCGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATc                                                        >  3:834131/1‑151 (MQ=255)
                                                                                            aGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCt                                                       >  2:1645107/1‑151 (MQ=255)
                                                                                            aGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCt                                                       <  4:1037074/151‑1 (MQ=255)
                                                                                            aGTTTCTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCt                                                       >  6:1422156/1‑151 (MQ=255)
                                                                                             gTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCtt                                                      >  5:173540/1‑151 (MQ=255)
                                                                                             gTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCtt                                                      >  7:1136438/1‑151 (MQ=255)
                                                                                                tGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTc                                                    >  7:376128/1‑150 (MQ=255)
                                                                                                tGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTc                                                    <  7:1371157/150‑1 (MQ=255)
                                                                                                tGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTc                                                    <  2:581633/150‑1 (MQ=255)
                                                                                                 gTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCa                                                   <  1:707496/150‑1 (MQ=255)
                                                                                                   aGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGt                                                <  2:368824/151‑1 (MQ=255)
                                                                                                   aGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGt                                                <  8:1343451/151‑1 (MQ=255)
                                                                                                      ggATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACt                                             >  4:758178/1‑151 (MQ=255)
                                                                                                         tCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGt                                          >  8:1089322/1‑151 (MQ=255)
                                                                                                           ttCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTAt                                        >  8:1375328/1‑151 (MQ=255)
                                                                                                            tCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTAt                                        <  3:548494/150‑1 (MQ=255)
                                                                                                            tCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGCGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTg                                            <  7:1194978/146‑1 (MQ=255)
                                                                                                            tCGCCCATAATCTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGt                                          <  1:493654/148‑1 (MQ=255)
                                                                                                             cGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATcc                                      >  8:455079/1‑151 (MQ=255)
                                                                                                              gCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATccc                                     <  8:233901/151‑1 (MQ=255)
                                                                                                               cccATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCg                                    <  7:190579/151‑1 (MQ=255)
                                                                                                                ccATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGTGGCAGCTCGCTGATGCGATAGACGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCgg                                   >  3:179451/1‑151 (MQ=255)
                                                                                                                ccATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCgg                                   >  7:1601334/1‑151 (MQ=255)
                                                                                                                  aTAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGAt                                 >  8:397026/1‑151 (MQ=255)
                                                                                                                  aTAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGAt                                 >  7:445202/1‑151 (MQ=255)
                                                                                                                     aTGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACg                              >  2:1297569/1‑151 (MQ=255)
                                                                                                                         aCGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTAt                           <  7:989560/150‑1 (MQ=255)
                                                                                                                         aCGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTAt                           <  5:1065898/150‑1 (MQ=255)
                                                                                                                              gcacgcacTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCa                     <  8:1562922/151‑1 (MQ=255)
                                                                                                                              gcacgcacTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCa                     <  5:518881/151‑1 (MQ=255)
                                                                                                                               cacgcacTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAg                    <  8:1044047/151‑1 (MQ=255)
                                                                                                                                acgcacTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAgc                   <  1:1680263/151‑1 (MQ=255)
                                                                                                                                  gcacTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAgcgc                 >  4:1395609/1‑151 (MQ=255)
                                                                                                                                      ttCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGccc             >  4:238856/1‑151 (MQ=255)
                                                                                                                                       tCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGccc             <  1:135333/150‑1 (MQ=255)
                                                                                                                                         gCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCata          <  3:1164447/151‑1 (MQ=255)
                                                                                                                                         gCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCata          <  6:1549714/151‑1 (MQ=255)
                                                                                                                                          ccATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCatat         <  3:439904/151‑1 (MQ=255)
                                                                                                                                               cgctgtgcGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAg    <  2:701255/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GCCGCCGACGCGGTTTTGGTCATAAAGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGCGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGC  >  NC_000913/4298051‑4298346

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