Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 695,330 G→A C77Y (TGC→TAC)  ubiF → 3‑demethoxyubiquinol 3‑hydroxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913695,3300GA86.8% 156.1 / 15.7 68C77Y (TGC→TAC) ubiF3‑demethoxyubiquinol 3‑hydroxylase
Reads supporting (aligned to +/- strand):  ref base G (5/4);  new base A (33/26);  total (38/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.97e-02

AGGATGATTGACCATGACAAATCAACCAACGGAAATTGCCATTGTCGGCGGAGGAATGGTCGGCGGCGCACTGGCGCTGGGGCTGGCACAGCACGGATTTGCGGTAACGGTGATCGAGCACGCAGAACCAGCGCCGTTTGTCGCTGATAGCCAACCGGACGTGCGGATCTCGGCGATCAGCGCGGCTTCGGTATCATTGCTTAAAGGGTTAGGGGTCTGGGATGCAGTACAGGCTATGCGTTGCCATCCTTACCGCAGACTGGAAACGTGGGAGTGGGAAACGGCGCATGTGGTGTTTGACGCCGCTGAACTTAAGCTACCGCTGCTTGGCTATATGGTGGAAAACACTGTCCTGCAACAGGCGTTGTGGCAGGCGCTGGAAGCGCATCCGAAAGTAACGTTACGTGTGCCAGGCTCGCTGATTGCGCTGCATCGCCATGATGATCTTCAGGAGCTGGAGCTGAAAGGCGGTGAAGTGATTCGC  >  NC_000913/695088‑695571
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                 
aGGATGATTGACCATGACAAATCAACCAACGGAAATTGCCATTGTCGGCGGAGGAATGGTCGGCGGCGCACTGGCGCTGGGGCTGGCACAGCACGGATTTGCGGTAACGGTGATCGAGCACGCAGAACCAGCGCCGTTTGTCGCTGATAGCCAACCGGACGTGCGGATCTCGGCGATCAGCGCGGCTTCGGTATCATTGCTTAAAGGGTTAGGGGTCTGGGATGCAGTACAGGCTATGCGTTACCATCCtt                                                                                                                                                                                                                                           >  2:51596/1‑251 (MQ=255)
  gatgatTGACCATGACAAATCAACCAACGGAAATTGCCATTGTCGGCGGAGGAATGGTCGGCGGCGCACTGGCGCTGGGGCTGGCACAGCACGGATTTGCGGTAACGGTGATCGAGCACGCAGAACCAGCGCCGTTTGTCGCTGATAGCCAACCGGACGTGCGGATCTCGGCGATCAGCGCGGCTTCGGTATCATTGCTTAAAGGGTTAGGGGTCTGGGATGCAGTACAGGCTATGCGTTGCCATCCTTAc                                                                                                                                                                                                                                         >  2:528516/1‑251 (MQ=255)
         gACCATGACAACTCAACCAACGGAACTTGCCATTGTCGGCGGAGGAATGGTCGGCGGCGCACTGGCGCTGGGGCTGGCACAGCACGGATTTGCGGTAACGGTGATCGAGCACGCAGAACCAGCGCCGTTTGTCGCTGATAGCCAACCGGACGTGCGGATCTAGGCGATCAGCGCGGCTTCGGTATCATTGCTTAAAGGGTTAGGGGTCTGGGATGCAGTACAGGCTATGCGTTACCATCCTTACCGCAGAc                                                                                                                                                                                                                                  >  2:33719/1‑251 (MQ=255)
           ccATGACAAATCAACCAACGGAAATTGCCATTGTCGGCGGAGGAATGGTCGGCGGCGCACTGGCGCTGGGGCTGGCACAGCACGGATTTGCGGTAACGGTGATCGAGCACGCAGAACCAGCGCCGTTTGTCGCTGATAGCCAACCGGACGTGCGGATCTCGGCGATCAGCGCGGCTTCGGTATCATTGCTTAAAGGGTTAGGGGTCTGGGATGCAGTACAGGCTATGCGTTACCATCATTACCGCAGACTg                                                                                                                                                                                                                                >  1:581155/1‑251 (MQ=255)
            cATGACAAATCAACCAACGGAAATTGCCATTGTCGGCGGAGGAATGGTCGGCGGCGCACTGGCGCTGGGGCTGGCACAGAACGGATTTGCGGTAACGGTGATCGAGCACGCAGAACCAGCGCCGTTTGTCGCTGATAGCCAACCGGACGTGCGGATCTCGGCGATCAGCGCGGCTTCGGTCTCATTGCATAAAGGGTTAGGGGTCTGGGATGCAGTACAGGCTATGCGTTACCATCCTTACCGCAGACTgg                                                                                                                                                                                                                               >  2:481925/1‑251 (MQ=255)
                 cAAATCAACCAACGGAAATTGCCATTGTCGGCGGAGGAATGGTCGGCGGCGCACTGGCGCTGGGGCTGGCACAGCACGGATTTGCGGTAACGGTGATCGAGCACGCAGAACCAGCGCCGTTTGTCGCTGATAGCCAACCGGACGAGCGGATCTCGGCGATCAGCGCGGCTTCGGTATCATTGCTTTAAGGGTTAGGGGTCTGGGATGCAGTACCGGCTATGCGTTACCATCCTTACCGCAGACTGGAAACg                                                                                                                                                                                                                          >  2:111792/1‑251 (MQ=255)
                  aaaTCAACCACCGGAAATTGCCATTGTCGGAGAAGGACTGCTCGGCGGCGCACTCGCGATGCAGCTGGCACAACACGGATTAGCAGTACCGGTGATCGAGAACGCAGACCCAGCGCCGTTTGTCGCTGATAGCCAACCGGACGTGCGGATCTCGGCGATCAGCGCGGCTTCGGTATCATTGCTTAACGGGTTAGGGGTCTGGGATGCAGTACAGGCTATGCGTTACCATCCTTACCGCAGACTGGaaaagg                                                                                                                                                                                                                         >  2:292145/1‑248 (MQ=255)
                  aaaTCAACCACCGGAAAATGACATTGTCTGCGGCGGAATGGGCGACGGCGCACTGGAGCAGAGGCAGGCACAGCACGGATTTGCGGGAACGGTGATCAAGTACGAAGAACCAGAGCCGTTTGTCGCTGATAGCCAACCGGACGTGCGGATCTCGGAGCTCAGCGCGGATGCGGGATCCTTGCTTAGAGGGTTAGGGGTCTGAGATGCAGTACAGGCCATGCGTTAACGTCCTTACCgacgacttgcaacgt                                                                                                                                                                                                                         >  2:371900/1‑237 (MQ=255)
                       aaccaacGGAAATTGCCATTGTCGGCGGAGGAATGGTCGGCGGCGCACTGGCGCTGGGGCTGGCACAGCACGGATTTGCGGTAACGGTGATCGAGCACGCAGAACCAGCGCCGTTTGTCGCTGATAGCCAACCGGACGTGCGGATCTCGGCGATCAGCGCGGCTTCGGTATCATTGCTTAAAGGGTTAGGGGTCTGGGATGCAGTACAGGCTATGCGTTACCATCCTTACCGCAGACTGGAAACgtggg                                                                                                                                                                                                                      >  1:13291/1‑249 (MQ=255)
                              ggAAATTGCCATTGTCGGCGGAGGAATGGTCGGCGGCGCACTGGCGCTGGGGCTGGCACAGCACGGATTTGCGGTAACGGTGATCGAGCACGCAGAACCAGCGCCGTTTGTCGCTGATAGCCAACCGGACGTGCGGATCTCGGCGTTCAGCGCGGCTTCGGTATCATTGCTTAAAGGGTTAGGGGTCTGGGATGCAGTACAGGCTATGCGTTACCATc                                                                                                                                                                                                                                              >  1:331010/1‑218 (MQ=255)
                              ggAAATTGCCATTGTCGGCGGAGGAATGGGCGGCGGCGCACTGGCGCTGGGGCTGGCACAGCACGGATTTGCGGTAACGGTGATCGAGCACGCAGAACCAGCGCCGTTTGTCGCTGATAGCCAACCGGACGTGCGGATCTCGGCGTTCAGCGCGGCTTCGGTATCATTGCTTAAAGGGTTAGGGGTCTGGGATGCAGTACAGGCTATGCGTTACCATc                                                                                                                                                                                                                                              <  2:331010/218‑1 (MQ=255)
                                  gttGCCAGTGTCGGCGGAGGGATGGGCGGTGGCGCACGGGCTCCGGGGCTGGCACAGCACGGATTTGCGGTAACGGGGATCGAGCACGCAGAACCAGGGCCGTTTGTCGCTGATAGCCAACCGGACGTGCTGATCTCGGCGATCAGCGCGGCTTCGGTATCATTGGTTAAAGGGTTAGGGGTCTGGGATGCAGTACAGGCTATGCGTTACCAACCTTGCCGCAGACTGGAAACGTGGGAGAGGGAAGCGgc                                                                                                                                                                                                         <  2:438027/250‑1 (MQ=255)
                                    tGCCATTGTGGACGGTGGAATGGTCGGCGGCGCACTGGCGCGTGGGCGGGCACAGCACTGATTTGCGGTAACGGGGATCGAGCCCGCAGAATCAGCGCCGTTTGTCGCTGATAGCCAACCGGACGTTCGGATCTCGGCGATCAGCGCGGCTTCGGTATCATTGCTTAAAGGGTTAGGGGTCTGGGATGCAGTACAGGCTATGCGTTACCATCCTTACCGCAGACTGGAAACGTGGGAGTGGGAAACGgcg                                                                                                                                                                                                        <  2:80145/250‑1 (MQ=255)
                                      ccATTGTCGGCGGAGGAATGGTCGGCGGCGCACTGGCGCTGGGGCTGGCACAGCACGGATTTGCGGTAACGGTGATCGAGCACGCAGAACCAGCGCCGTTTGTCGCTGATAGCCAACCGGACGTGCGGATCTCGGCGATCAGCGCGGCTTCGGTATCATTGCTTAAAGGGTTAGGGGTCTGGGATGCAGTACAGGCTATGCGTTACCATCCTTACCGCAGACTGGAAACGTGGGAGTGTGAAACGGCGCAt                                                                                                                                                                                                     <  2:69062/251‑1 (MQ=255)
                                           gTCGGCGGAGGAATGGTCGGCGGCGCACTGGCGCTGGGGCTGGCACAGCACGGATTTGCGGTAACGGTGATCGAGCACGCAGAACCAGCGCCGTTTGTCGCTGATAGCCAACCGGACGTGCGGATCTCGGCGATCAGCGCGGCTTCGGTATCATTGCTTAAAGGGTTAGGGGTCTGGGATGCAGTACAGGCTATGCGTTACCATCCTTACCGCAGACTGGAAACGTGGGAGTGGGAAACGGCGCATgtggt                                                                                                                                                                                                >  1:19040/1‑251 (MQ=255)
                                                cggAGGAATGGTCGGCGGCGCACTGGCGCTGGGGCTGGCACAGCACGGATTTGCGGTAACGGTGATCGAGCACGCAGAACCAGCGCCGTTTGTCGCTGATAGCCAACCGGACGTGCGGATCTCGGCGATCAGCGCGGCTTCGGTATCATTGCTTAAAGGGTTAGGGGTCTGGGATGCAGTACAGGCTATGCGTTGCCATCCTTACCGCAGACTGGAAACGTGGGAGTGGGAAACGGCGCATGTGGTGTTTg                                                                                                                                                                                           >  1:332113/1‑251 (MQ=255)
                                                 ggaggcATGGTCGGCGGCGCACTGGCGCTGGGGCTGGCACAGCACGGATTTGCGGTAACGGTGATCGAGCACGCAGAACCAGCGCCGTTTGTCGCTGATAGCCAACCGGACGTGCGGATCTCGGCGATCAGCGCGGCTTCGGTATCATTGCTTAAAGGGTTAGGGGTCTGGGATGCAGTACAGGCTATGCGTTACCATCCTTACCGCAGACTGGAAACGTGGGGGTGGGAAACGGCGCATGTGGTGGTTgg                                                                                                                                                                                          >  2:325795/1‑250 (MQ=255)
                                                 ggaggaATGGTCGGCGGCGCACTGGCGCTGGGGCTGGCACAGCACGGATTTGCGGTAACGGTGATCGAGCACGCAGAACCAGCGCCGTTTGTCGCTGATAGCCAACCGGACGTGCGGATCTCGGCGATCAGCGCGGCTTCGGTATCATTGCTTCAAGGGTTAGGGGTCTGGGATGCAGTACAGGCTATGCGTTACCATCCTTACCGCAGACTGGACACGTGGGAGTGGGAAACGGCGCCTGTGGTGTTTgt                                                                                                                                                                                          >  2:450620/1‑250 (MQ=255)
                                                 ggaggaATGGTCGGCGGCGCACTGGCGCTGGGGCTGGCACAGCACGGATTTGCGGTAACGGTGATCGAGCACGCAGAACCAGCGCCGTTTGTCGCTGATAGCCAACCGGACGTGCGGATCTCGGCGATCAGCGCGGCTTCGGTATCATTGCTTAAAGGGTTAGGGGTCTGGGATGCAGTACAGGCTATGCGTTGCCATCCTTACCGCAGACTGGAAACGTGGGAGTGGGAAACGGCGCATGTGGTGTTTg                                                                                                                                                                                           >  1:389668/1‑250 (MQ=255)
                                                 ggaggaATGGTCGGCGGCGCACTGGCGCTGGGGCTGGCACAGCACGGATTTGCGGTAACGGTGATCGAGCACGCAGAACCAGCGCCGTTTGTCGCTGATAGCCAACCGGACGTGCGGATCTCGGCGATCAGCGCGGCTTCGGTATCATTGCTTAAAGGGTTAGGGGTCTGGGATGCAGTACAGGCTATGCGTTACCATCCTTACCGCAGACTGGAAACGTGGGAGTGGGAAACGGCGCATGTGGGGTTTg                                                                                                                                                                                           >  2:369217/1‑250 (MQ=255)
                                                             ggctgcgcACTGGCGCTGGGGCTGGCACAGCACGGATTTGCGGTAACGGTGATCGAGCACGCAGAACCAGCGCCGTTTGTCGCTGATAGCCAACCGGACGTGCGGATCTCGGCGATCAGCGCGGCTTCGGTATCATTGCTTAAAGGGTTAGGGGTCTGGGATGCAGTACAGGCTATGCGTTACCATCCTTACCGCAGACTGGAAACGTGGGAGTGGGAAACGGCGCATGTGGTGTTTGACGCCGCTGAACt                                                                                                                                                                              >  2:156073/5‑251 (MQ=255)
                                                             ggcggcGCACTGGCGCTGGGGCTGGCACAGCACGGATTTGCGGTAACGGTGATCGAGCACGCAGAACCAGCGCCGTTTGTCGCTGATAGCCAACCGGACGTGCGGATCTCGGCGATCAGCTCGGCTTCGGTATCATTGCTTAAAGGGTTAGGGGTCTGGGATGCAGTACAGGCTATGCGTTACCATCCTTACCGCAGACTGGAAACGTGGGAGTGGGAAACGGCGCATGTGGTGTTTGACGCCGCTGAACt                                                                                                                                                                              >  1:152065/1‑251 (MQ=255)
                                                              gcggcgCACTGGCGCTGGGGCTGGCACAGCACGGATTTGCGGTAACGGTGATCGAGCACGCAGAACCAGCGCCGTTTGTCGCTGATAGCCAACCGGACGTGCGGATCTCGGCGATCAGCGCGGCTTCGGTATCATTGCTTAAAGGGTTAGGGGTCTGGGATGCAGTACAGGCTATGCGTTACCATCCTTACCGCAGACTGGAAACGTGGGAGTGGGAAACGGCGCATGTGGTGTTTGACGCCGCTGAACtt                                                                                                                                                                             <  1:617422/251‑1 (MQ=255)
                                                                        ggCGCTGGGGCTGGCACAGCACGGATTTGCGGTAACGGTGATCGAGCACGCAGAACCAGCGCCGTTTGTCGCTGATAGCCAACCGGACGTGCGGATCTCGGCGATCAGCGCGGCTTCGGTATCATTGCTTAAAGGGTTAGGGGTCTGGGATGCAGTACAGGCTATGCGTTACCATCCTTACCGCAGACTGGAAACGTGGGAGTGGGAAACGGCGCATGTGGTGTTTGACGCCGCTGAACTTAAGCTAcc                                                                                                                                                                     <  1:338531/249‑1 (MQ=255)
                                                                         gcgcTGGGGCTGGCACAGCACGGATTTGCGGTAACGGTGATCGAGCACGCAGAACCAGCGCCGTTTGTCGCTGATAGCCAACCGGACGTGCGGATCTCGGCGATCAGCGCGGCTTCGGTATCATTGCTTAAAGGGTTAGGGGTCTGGGATGCAGTACAGGCTATGCGTTGCCATCCTTACCGCAGACTGGAAACGTGGGAGTGGGAAACGGCGCATGTGGTGTTTGACGCCGCTGAACTTAAGCTACCgct                                                                                                                                                                  <  2:332113/251‑1 (MQ=255)
                                                                         gcgcTGGGGCTGGCACAGCACGGATTTGCGGTAACGGTGATCGAGCACGCAGAACCAGCGCCGTTTGTCGCTGATAGCCAACCGGACGTGCGGATCTCGGCGATCAGCGCGGCTTCGGTATCATTGCTTAAAGGGTTAGGGGTCTGGGATGCAGTACAGGCTATGCGTTGCCAGCCTTACCGCAGACTGGAAACGTGGGAGTGGGAAACGGCGCATGTGGTGTTTGACGCCGCTGAACTTAAGCTACCgct                                                                                                                                                                  <  1:528516/251‑1 (MQ=255)
                                                                               ggacTGGCACAGCACGGATTTGCGGTAACGGAGATCGAGCCCGCCGAACCAGCGCCGTTTGTCGCTGATAGCCAACCGGACGTGCGGATCTCGGCGATCAGCGCGGCTTAGGTATCATTGCTTAAAGGGTTAGGGGTCTGGGATGCAGTACAGGCTATGCGTTACCATCATTACCGCCGACTGGAAACGGGGGAGTGGGAAACGGCGCATGTGGTGTTTGACGCCGCTGAACTTAAGCTACCGCTGCTTgg                                                                                                                                                            >  2:363124/4‑251 (MQ=255)
                                                                                 gCTGGCACAGCACGGATTTGCGGTAACGGTGATCGAGCACGCAGAACCAGCGCCGTTTGTCGCTGATAGCCAACCGGACGTGCGGATCTCGGCGATCAGCGCGGCTTCGGTATCATTGCTTAAAGGGTTAGGGGTCTGGGATGCAGTACAGGCTATGCGTTACCATCCTTACCGCAGACTGGAAACGTGGGAGTGGGAAACGGCGCATGTGGTGTTTGACGCCGCTGAACTTAAGCTACCGCTGCTTGGCt                                                                                                                                                          <  1:494467/251‑1 (MQ=255)
                                                                                    ggCACAGCACGGATTTGCGGTAACGGTGATCGAGCACGCAGAACCAGCGCCGTTTGTCGCTGATAGCCAACCGGACGTGCGGATCTCGGCGATCAGCGCGGCTTCGGTATCATTGCTTAAAGGGTTAGGGGTCTGGGATGCAGTACAGGCTATGCGTTACCATCCTTACCGCAGACTGGAAACGTGGGAGTGGGAAACGGCGCATGTGGTGTTTGACGCCGCTGAACTTAAGCTACCGCTGCTTGGCtata                                                                                                                                                       <  1:369217/251‑1 (MQ=255)
                                                                                      cacaGCACGGATTTGCGGTAACGGTGATCGAGCACGCAGAACCAGCGCCGTTTGTCGCTGATAGCCAACCGGACGTGCGGATCTCGGCGATCAGCGCGGCTTCGGTATCATTGCTTAAAGGGTTAGGGGTCTGGGATGCAGTACAGGCTATGCGTTACCATCCTTACCGCAGACTGGAAACGTGGGAGTGGGAACCGGCGCATGTGGTGTTTGACGCCGCTGAACTTAAGCTACCGCTGCTTGGCTATAtg                                                                                                                                                     >  1:294050/1‑251 (MQ=255)
                                                                                           cACGGATTTGCGGTAACGGTGATCGAGCACGCAGAACCAGCGCCGTTTGTCGCTGATAGCCAACCGGACGTGCGGATCTCGGCGATCAGCGCGGCTTCGGTATCATTGCTTAAAGGGTTAGGGGTCTGGGATGCAGTACAGGCTATGCGTTACCATCCTTACCGCAGACTGGAAACGTGGGAGTGGGAAACGGCGCATGTGGTGTTTGACGCCGCTGAACTTAAGCTACCGCTGCTTGGCTATAtggtgg                                                                                                                                                 >  1:213974/1‑250 (MQ=255)
                                                                                           cACGGATTTGCGGTAACGGTGATCGAGCACGCAGAACCAGCGCCGTTTGTCGCTGATAGCCAACCGGACGTGCGGATCTCGGCGATCAGCGCGGCTTCGGTATCATTGCTTAAAGGGTTAGGGGTCTGGGATGCAGTACAGGCTATGCGTTACCATCCTTACCGCAGACTGGAAACGTGGGAGTGGGAAACGGCGCATGTGGTGTTTGACGCCGCTGAACGTAAGCTACCGCTGCTTGGCTATATGGTGGa                                                                                                                                                <  1:499310/251‑1 (MQ=255)
                                                                                               gATTTGCGGTAACGGTGATCGAGCACGCAGAACCAGCGCCGTTTGTCGCTGATAGCCAACCGGACGTGCGGATCTCGGCGATCAGCGCGGCTTCGGTATCATTGCTTAAAGGGTTAGGGGTCTGGGATGCAGTACAGGCTATGCGTTACCATCCTTACCGCAGACTGGAAACGTGGGAGTGGGAAACGGCGCATGTGGTGTTTGACGCCGCTGAACTTAAGCTACCGCTGCTTGGCTATATGGTGGAAAac                                                                                                                                            <  1:431078/251‑1 (MQ=255)
                                                                                                         aaCGGTGATCGAGCACGCAGAACCAGCGCCGTTTGTCGCTGATAGCCAACCGGACGTGCGGATCTCGGCGATCAGCGCGGCTTCGGTATCATTGCTTAAAGGGTTAGGGGTCTGGGATGCAGTACAGGCTATGCGTTACCATCCTTACCGCAGACTGGAAACGTGGGAGTGGGAAACGGCGCATGTGGTGTTTGACGCCGCTGAACTTAAGCTACCGCTGCTTGGCTATATGGTGGAAAACACTGTCCTGc                                                                                                                                  <  1:274748/251‑1 (MQ=255)
                                                                                                           cGGTGATCGCGCACGCAGAACCAGTGCCGTTTGGCGCTGATAGCCAACCGGACGTGCGGATCTCGGGGATCAGCGCGGCTTCGGTATCATTGCGTAAAGGGTTAGGGGTCTGGGAGGCAGTACAGGCTATGCGTTACCATCCTTACCGCAGACGTGAAACGTGGGAGGGGGAAACGGAGCATGTGGTGGTTGACGCCGCGGACCTTAAGATACAGCTGCTTGGCTATATGGTGGAAAACACTGTccggtaa                                                                                                                                <  2:184891/251‑6 (MQ=255)
                                                                                                            ggTGATCGAGTACGCAGACCCAGCGCCGTTTGGCGCTGATAGCCAACCGGACGTGCGGATCTCGTCGATAAGCGCGGCTTCGGTATCATTGCTTAAAGGGTTTGGGGTCTGGGATGCAGTACAGGCTATGCGTTACCATCCTTACCGCAGACTGGAAACGTGGGAGTGGGAAACGGCGCATGTGGTGTTTGACGCCGCTGAACTTAAGCTACCGCTGCTTGGCTATATGGTGGAAAACACTTTCCTGCAAc                                                                                                                               <  2:354299/251‑1 (MQ=255)
                                                                                                            ggTGATCGAGCACGCAGAACCAGCGCCGTTTGTCGCTGATAGCCAACCGGACGTGCGGATCTCGGCGATCAGCGCGGCTTCGGTATCATTGCTTAAAGGGTTAGGGGTCTGGGATGCAGTACAGGCTATGCGTTACCATCCTTACCGCAGACTGGAAACGTGGGAGTGGGAAACGGCGCATGTGGTGTTTGACGCCGCTGAACTTAAGCTACCGCTGCTTGGCTATATGGTGGAAAACACTGTCCTGCAAc                                                                                                                               >  1:466198/1‑251 (MQ=255)
                                                                                                              tGATCGAGCACGCAGAACCAGCGCCGTTTGTCGCTGATAGCCAACCGGACGTGCGGATCTCGGCGATCAGCGCGGCTTCGGTATCATTGCTTAAAGGGTTAGGGGTCTGGGATGCAGTACAGGCTATGCGTTACCATCCTTACCGCAGACTGGAAACGTGGGAGTGGGAAACGGCGCATGTGGTGTTTGACGCCGCTGAACTTAAGCTACCGCTGCTTGGCTATATGGTGGAAAACACTGTCCTGCAAc                                                                                                                               >  1:216418/1‑249 (MQ=255)
                                                                                                              tGATCGAGCACGCAGAACCAGCGCCGTTTGTCGCTGATAGCCAACCGGACGTGCGGATCTCGGCGATCAGCGCGGCTTCGGTATCATTGCTTAAAGGGTTAGGGGTCTGGGATGCAGTACAGGCTATGCGTTACCATCCTTACCGCAGACTGGAAACGTGGGAGTGGGAAACGGCGCATGTGGTGTTTGACGCCGCTGAACTTAAGCTACCGCTGCTTGGCTATATGGCGGAAAACACTGTCCTGCAAc                                                                                                                               >  2:375659/1‑249 (MQ=255)
                                                                                                                  cGAGCACGCAGAACCAGCGCCGTTTGTCGCTGATAGCCAACCGGACGTGCGGATCTCGGCGATCAGCGCGGCTTCGGTATCATTGCTTAAAGGGTTAGGGGTCTGGGATGCAGTACAGGCTATGCGTTGCCATCCTTACCGCAGACTGGAAACGTGGGAGTGGGAAACGGCGCATGTGGTGTTTGACGCCGCTGAACTTAAGCTACCGCTGCTTGGCTATATGGTGGAAAACACTGTCCTGCAACAGGCGt                                                                                                                         >  1:19758/1‑251 (MQ=255)
                                                                                                                       aCGCAGAACCAGCGCCGTTTGTCGCTGATAGCCAACCGGACGTGCGGATCTCGGCGATCAGCGCGGCTTCGGTATCATTGCTTAAAGGGTTAGGGGTCTGGGATGCAGTACAGGCTATGCGTTACCATCCTTACCGCAGACTGGAAACGTGGGAGTGGGAAACGGCGCATGTGGTGTTTGACGCCGCTGAACTTAAGCTACCGCTGCTTGGCTATATGGTGGAAAACACTGTCCTGCAACAGGCGTTGTgg                                                                                                                    >  1:334244/1‑251 (MQ=255)
                                                                                                                       aCGCAGAACCAGCGCCGTTTGTCGCTGATAGCCAACCGGACGTGCGGATCTCGGCGATCAGCGCGGCTTCGGTATCATTGCTTAAAGGGTGAGGGGTCTGGGATGCAGTACAGGCTATGCGTTACCATCCTTACCGCAGACTGGAAACGTGGGAGTGGGAAACGGCGCATGTGGTGTTTGACGCCGCTGAACTTAAGCTACCGCTGCTTGGCTATATGGTGGAAAACACTGTCCTGCAACAGGCGTTGTgg                                                                                                                    <  2:131492/251‑1 (MQ=255)
                                                                                                                             aaCCAGCGCCGTTTGTCGCTGATAGCCACCCGGACGTGCGGATCTCGGCGATCAGCGCGGATTCGGTATCATTGCTTAAAGGGTTAGGGGTCTGGGATGCCGTACAGGCTCTGCGTTACCATCCTTACCGCAGACTGGAAACGTGGGAGTGGGAAACGGCGCATGTGGTGTTTGACGCCGCTGAACTTAAGCTACCGCTGCTTGGCTATATGGTGGAAAACACTGTCCTGCAACAGGCGTTGTGGCAGgcg                                                                                                              >  2:396671/1‑251 (MQ=255)
                                                                                                                                  gcgcCGTTTGTCGCTGAGAGCCAACCGGACGTGCGGATCTCGGCGATCAGCGCGGCTTCGGTATCATTGCTTAAAGGGTTAGGGGTCTGGGATGCAGGACAGGCTATGCGTTACCATCCTTACCGCAGACTGGAAACGTGGGAGTGGGAAACGGCGCATGTGGTGTTGGTCGCCGCTGAACTTAAGCTACCGCTGCTTGGCTATATGGTGGAAAACACTGTCCTGCAACAGGCGTTGTGGCAGGCGCTGGa                                                                                                         <  2:486897/251‑1 (MQ=255)
                                                                                                                                       gTTTGTCGCTGATAGCCAACCGGACGTGCGGATCTCGGCGATCAGCGCGGCTTCGGTATCATTGCTTAAAGGGTTAGGGGTCTGGGATGCAGTACAGGCTATGCGTTACCATCCTTACCGCAGACTGGAAACGTGGGAGTGGGAAACGGCGCATGTGGTGTTTGACGCCGCTGAACTTAAGCTACCGCTGCTTGGCTATATGGTGGAAAACACTGTCCTGCAACAGGCGTTGTGGCAGGCGCTGGAAgcgc                                                                                                    <  1:246126/251‑1 (MQ=255)
                                                                                                                                         ttGTCGCTGATAGCCAACCGGACGTGCGGATCTCGGCGATCAGCGCGGCTTCGGTATCATTGCTTAAAGGGTTAGGGGTCTGGGATGCAGTACAGGCTATGCGTTACCATCCTTACCGCAGACTGGAAACGTGGGAGTGGGAAACGGCGCATGTGGTGTTTGACGCCGCTGAACTTAAGCTACCGCTGCTTGGCTATATGGTGGAAAACACTGTCCTGCAACAGGCGTTGTGGCAGGGGCGGGAAGCGCAt                                                                                                  >  2:398275/1‑251 (MQ=255)
                                                                                                                                             cGCTGATAGCCAACCGGACGTGCGGATCTCGGCGATCAGCGCGGCTTCGGTATCATTGCTTAAAGGGTTCGGGGTCTGGGATGCAGTACAGGCTATGCGTTACCATCCTTACCGCAGACTGGAAACGTGGGAGTGGGAAACGGCGCATGTGGTGTTTGACGCCGCTGAACTTAAGCTACCGCTGCTTGGCTATATGGTGGAAAACACTGTCCTGCAACAGGCGTTGTGGCAGGCGCTGGAAGGGCATCCg                                                                                               >  2:44454/1‑250 (MQ=255)
                                                                                                                                                    aGCCAACCGGACGTGCGGATCTCGGCGATCAGCGCGGCTTCGGTATCATTGCTTAAAGGGTTAGGGGTCTGGGATGCAGTACAGGCTATGCGTTACCATCCTTACCGCAGACTGGAAACGTGGGAGTGGGAAACGGCGCATGTGGTGTTTGACGCCGCTGAACTTAAGCTACCGCTGCTTGGCTATATGGTGGAAAACACTGTCCTGCAACAGGCGTTGTGGCAGGCGCTGGAAGCGCATCCGAAAGTAAc                                                                                       <  1:461758/251‑1 (MQ=255)
                                                                                                                                                    aGCCAACCGGACGTGCGGATCTCGGCGATCAGCGCGGCTTCGGTATCATTGCTTAAAGGGTTAGGGGTCTGGGATGCAGTACAGGCTATGCGTTACCATCCTTACCGCAGACTGGAAACGTGGGAGTGGGAAACGGCGCATGTGGTGTTTGACGCCGCTGAACTTAAGCTACCGCTGCTTGGCTATATGGCGGAAAACACTGTCCTGCAACAGGCGTTGTGGCAGGCGCTGGAAGCGCATCCGAAAGTaa                                                                                        <  1:548477/250‑1 (MQ=255)
                                                                                                                                                       cAACCGGACGTGCGGATCTCGGCGATCAGCGCGGCTTCGGTATCATTGCTTAAAGGGTTAGGGGTCTGGGATGCAGTACAGGCTATGCGTTACCATCCTTACCGCAGACTGGAAACGTGGGAGTGGGAAACGGCGCATGTGGTGTTTGACGCCGCTGAACTTAAGCTACCGCTGCTTGGCTATATGGTGGAAAACACTGTCCTGCAACAGGCGTTGTGGCAGGCGCTGGAAGCGCATCCGAAAGTAACGtt                                                                                    <  1:51596/251‑1 (MQ=255)
                                                                                                                                                                  gCGGATCTCGGCGATCAGCGCGGCTTCGGTATCATTGCTTAAAGGGTTAGGGGTCTGGGATGCAGTACAGGCTATGCGTTACCATCCTTACCGCAGACTGGAAACGTGGGAGTGGGAAACGGCGCATGTGGTGTTTGACGCCGCTGAACGTAAGCTACCGCTGCTTGGCTATATGGTGGAAAACACTGTCCTGCAACAGGCGTTGTGGCAGGCGCTGGAAGCGCTTCCGAAAGTAACGTTACGTGTGCCAg                                                                         <  2:119788/251‑1 (MQ=255)
                                                                                                                                                                          cGGCGATCAGCGCGGCTTCGGTATCATTGCTTAAAGGGTTAGGGGTCTTGGATGCAGTACAGGCTATGCGTTACCATCCTTACCGCAGACTGGAAACGTGGGAGTGGGAAACGGCGCATGTGGTGTTTGACGCCGCTGAACTTAAGCTACCGCTGCTTGGCTATATGGTGGAAAACACTGTCCTGCAACAGGCGTTGTGGCAGGCGCTGGAAGCGCATCCGAAAGTAACGTTACGTGTGCCAGGCTCg                                                                    <  1:612184/248‑1 (MQ=255)
                                                                                                                                                                          cGGCGATCAGCGCGGCTTCGGTATCATTGCTTAAAGGGTTAGGGGTCTGGGATGCAGTACAGGCTATGCGTTACCATCCTTACCGCAGACTGGAAACGTGGGAGTGGGAAACGGCGCATGTGGTGTTTGACGCCGCTGAACTTAAGCTACCGCTGCTTGGCTATATGGTGGAAAACACTGTCCTGCAACAGGCGTTGTGGCAGGCGCTGGAAGCGCATCCGAAAGTAACGTTACGTGTGCCAGGCTc                                                                     <  1:78365/247‑1 (MQ=255)
                                                                                                                                                                          cGGCGATCAGCGCGGCTTCGGTATCATTGCTTAAAGGGTTAGGGGTCTGGGATGCAGTACAGGCTATGCGTTACCATCCTTACCGCAGACTGGAAACGTGGGAGTGGGAAACGGCGCATGTGGTGTTTGACGCCGCTGAACTTAAGCTACCGCTGCTTGGCTATATGGTGGAAAACACTGTCCTGCAACAGGCGTTGTGGCAGGCGCTGGAAGCGCATCCGAAAGTAACGTTACGTGTGCCAGGCTCg                                                                    <  1:33719/248‑1 (MQ=255)
                                                                                                                                                                                      cgGCTTCGGTATCATTGCTTAAAGGGTTAGGGGTCTGGGATGCAGTACAGGCTATGCGTTACCATCCTTACCGCAGACTGGAAACGTGGGAGTGGGAAACGGCGCATGTGGTGTTTGACGCCGCTGAACTTAAGCTACCGCTGCTTGGCTATATGGTGGAAAACACTGTCCTGCAACAGGCGTTGTGGCAGGCGCTGGAAGCGCATCCGAAAGTAACGTTACGTGTGCCAGGCTCGCTGATTGCGCTGCAt                                                     <  2:157097/251‑1 (MQ=255)
                                                                                                                                                                                                aTCATTGCTTAAAGGGTTAGGGGTCTGGGATGCAGTACAGGCTATGCGTTGCCATCCTTACCGCAGACTGGAAACGTGGGAGTGGGAAACGGCGCATGTGGTGTTTGACGCCGCTGAACTTAAGCTACCGCTGCTTGGCTATATGGTGGAAAACACTGTCCTGCAACAGGCGTTGTGGCAGGCGCTGGAAGCGCATCCGAAAGTAACGTTACGTGTGCCAGGCTCGCTGATTGCGCTGCATCGCCatgatg                                           <  1:487550/251‑1 (MQ=255)
                                                                                                                                                                                                  cATTGCTTAAAGGGTTAGGGGTCTGGGATGCAGTACAGGCTATGCGTTACCATCCTTACCGCAGACTGGAAACGTGGGAGTGGGAAACGGCGCATGTGGTGTTTGACGCCGCTGAACTTAAGCTACCGCTGCTTGGCTATATGGTGGAAAACACTGTCCTGCACCAGGCGTTGTGGCAGGCGCTGGACGCGCATCCGAAAGTAACGTTACGTGAGCCAGGCTCGCTGATTGCGGTGCATCGACATaatgat                                         >  1:54524/1‑251 (MQ=255)
                                                                                                                                                                                                           aaGGGTTATGGGTCTGGGATGCAGTACAGGCTATGCGTTACCATCCTTACCGCAGACTGGAAACGTGGGAGTGGGAAACGGCGCATGTGGTGTTTGACGCCGCTGAACTTAAGCTACCGCTGCTTGGCTATATGGTGGAAAACACTGTCCTGCAACAGGCGTTGTGGCAGGCGCTGGAAGCGCATCCGAAAGTAACGTTACGTGTGCCAGGCTCGCTGATTGCGCTGCATCGCCATGATGATCTTCAggag                                <  1:422772/251‑1 (MQ=255)
                                                                                                                                                                                                             gggTTAGGGGTCTGGGATGCAGTACAGGCTATGCGTTACCATCCTTACCGCAGACTGGAAACGTGGGAGTGGGAAACGGCGCATGTGGTGTTTGACGCCGCTGAACTTAAGCTACCGCTGCTTGGCTATATGGTGGAAAACACTGTCCTGCAACAGGCGTTGTGGCAGGCGCTGGAAGCGCATCCGAAAGTAACGTTACGTGTGCCAGGCTCGCTGATTGCGCTGCATCGCCATGATGATCTTCAggagct                              <  1:176158/251‑1 (MQ=255)
                                                                                                                                                                                                                   ggggTCTGGGATGCAGTACAGGCTATGCGTTGCCATCCTTACCGCAGACTGGAAACGTGGGAGTGGGAAACGGCGCATGTGGTGTTTGACGCCGCTGAACTTAAGCTACCGCTGCTTGGCTATCTGGTGGAAAACACTGTCCTGCAACAGGCGTTGTGGCAGGCGCTGGAAGCGCATCCGAAAGTAACGTTACGTGTGCCAGGCTCGCTGATTGCGCTGCATCGCCATGATGATCTTCAggagctggagct                        >  1:392394/1‑251 (MQ=255)
                                                                                                                                                                                                                     ggTCTGGGATGCAGTACAGGCTATGCGTTGCCATCCTTACCGCAGACTGGAAACGTGGGAGTGGGAAACGGCGCATGTGGTGTTTGACGCCGCTGAACTTAAGCTACCGCTGCTTGGCTATATGGTGGAAAACACTGTCCTGCAACAGGCGTTGTGGCAGGCGCTGGAAGCGCATCCGAAAGTAACGTTACGTGTGCCAGGCTCGCTGATTGCGCTGCATCGCCATGATGATCTTCAGGAGCTGGAGCTGa                      <  1:135386/251‑1 (MQ=255)
                                                                                                                                                                                                                     ggTCTGAGATGCAGTACAGGCTATGCGTTACCATCCTTACCGCAGACTGGAAACGTGGGAGTGGGAAACGGCGCATGTGGTGTTTGACGCCGCTGAACTTAAGCTACCGCTGCTTGGCTATATGGTGGAAAACACTGTCCTGCAACAGGCGTTGTGGCAGGCGCTGGAAGCGCATCCGAAAGTAACGTTACGTGTGCCAGGCTCGCTGATTGCGCTGCATCGCCATGATGATCTTCAGgagctggagctg                       >  1:366041/1‑250 (MQ=255)
                                                                                                                                                                                                                                  gTACCGGCTATGCGTTACCATCCTTACCGCAGACTGGAAACGTGGGAGTGGGAAACGGCGCATGTGGTGTTTGACGCCGCTGAACTTAAGCTACCGCTGCTTGGCTATATGGTGGAAAACACTGTCCTGCAACAGGCGTTGTGGCAGGCGCTGGGACCGCATCCGAAAGGAACGTTACGTGTGCCAGGCGCGCTGATTGCGGTGCATCGCCATGGTTATCTTCAGGAGGTGGAGCTGAACGGCGGGGAAgg         >  2:345780/1‑250 (MQ=255)
                                                                                                                                                                                                                                   tACAGGCTATGCGTTACCATCCTTACCGCAGACTGGAAACGTGGGAGTGGGAAACGGCGCATGTGGTGTTTGACGCCGCTGAACTTAAGCTACCGCTGCTTGGCTATATGGTGGAAAACACTGTCCTGCAACAGGCGTTGTGGCAGGCGCTGGAAGCGCATCCGAAAGTAACGTTACGTGTGCCAGGCTCGCTGATTGCGCTGCATCGCCATGATGATCTTCAGGAGCTGGAGCTGAAAGGCGGTGAAGTg        >  1:400105/1‑251 (MQ=255)
                                                                                                                                                                                                                                    aCAGGCTATGCGTTACCATCCTTACCGCAGACTGGAAACGTGGGAGTGGGAAACGGCGCATGTGGCGTTTGACGCCGCTGACCTTAAGCTACCGCTGCTTGGCTATATGGTGGAAAACACTGTCCTGCAACAGGAGTTGTGGCAGGCGCTACAAGCGCATCCGAAAGTAACGCTACGTGGGGCAGGCACGCCGATTGCGCTGCATCGCCATGATGCTCTTCTggagcgggcgttgaaaagaggggacgtgg       >  2:3916/1‑227 (MQ=255)
                                                                                                                                                                                                                                      aGGCTATGCGTTACCATCCTTACCGCAGACTGGAAACGTGGGAGTGGGAAACGGCGCATGTGGTGTTTGACGCCGCTGAACTTAAGCTACCGCTGCTTGGCTATATGGTGGAAAACACTGTCCTGCAACAGGCGTTGTGGCAGGCGCTGGAAGCGCATCCGAAAGTAAAGTTACGTGTGCCAGGCTCGCTGATTGCGCTGAATCGCCATGATGATCTTCAGGAGCTGGAGCTGAAAGGGGGGGAAGGGAtt     >  1:512397/1‑251 (MQ=255)
                                                                                                                                                                                                                                       ggCTATGCGTTACCATCCTTACCGCAGACTGGAAACGTGGGAGTGGGAAACGGCGCATGTGGTGTTTGACGCCGCTGAACTTAAGCTACCGCTGCTTGGCTATATGGTGGAAAACACTGTCCTGCAACAGGCGTTGTGGCAGGCGCTGGAAGCGCATCCGAAAGTAACGTTACGTGTGCCAGGCTCGCTGATTGCGCTGCATCGCCATGATGATCTTCAGGAGCTGGAGCTGAAAGGCGGTGAAGTGATTc    >  1:112910/1‑251 (MQ=255)
                                                                                                                                                                                                                                         cTATGCGTTACCATCCTTACCGCAGACTGGAAACGTGGGAGTGGGAAACGGCGCATGTGGTGTTTGACGCCGCTGAACTTAAGCTACCGCTGCTTGGCTATATGGTGGAAAACACTGTCCTGCAACAGGCGTTGTGGCAGGCGCTGGAAGCGCATCCGAAAGTAACGTTACGTGTGCCAGGCTCGCTGATTGCGCTGCATCGCCATGATGATCTTCAGGAGCTGGAGCTGAAAGGCGGTGAAGTGATTcgc  >  1:523346/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                 
AGGATGATTGACCATGACAAATCAACCAACGGAAATTGCCATTGTCGGCGGAGGAATGGTCGGCGGCGCACTGGCGCTGGGGCTGGCACAGCACGGATTTGCGGTAACGGTGATCGAGCACGCAGAACCAGCGCCGTTTGTCGCTGATAGCCAACCGGACGTGCGGATCTCGGCGATCAGCGCGGCTTCGGTATCATTGCTTAAAGGGTTAGGGGTCTGGGATGCAGTACAGGCTATGCGTTGCCATCCTTACCGCAGACTGGAAACGTGGGAGTGGGAAACGGCGCATGTGGTGTTTGACGCCGCTGAACTTAAGCTACCGCTGCTTGGCTATATGGTGGAAAACACTGTCCTGCAACAGGCGTTGTGGCAGGCGCTGGAAGCGCATCCGAAAGTAACGTTACGTGTGCCAGGCTCGCTGATTGCGCTGCATCGCCATGATGATCTTCAGGAGCTGGAGCTGAAAGGCGGTGAAGTGATTCGC  >  NC_000913/695088‑695571

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑