Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 438,743 C→A 7.3% G479C (GGC→TGC)  dxs ← 1‑deoxyxylulose‑5‑phosphate synthase, thiamine triphosphate‑binding, FAD‑requiring

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913438,7430CA 7.3% 382.1 / 14.4 109G479C (GGC→TGC) dxs1‑deoxyxylulose‑5‑phosphate synthase, thiamine triphosphate‑binding, FAD‑requiring
Reads supporting (aligned to +/- strand):  ref base C (51/50);  new base A (4/4);  total (55/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 112 total aligned reads displayed.

TGGCGTTCAGCGATTCGGCGACTTTCGCCGCTTCTGGCATCAGCGTACCAAAGTTAAGGATCGCCAGTTTCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGTCATAATGACCATTTCCGGTATGCAGCGCAGGTAAGAGAGATCAAAAGCACCCTGATGGGTT  >  NC_000913/438596‑438891
                                                                                                                                                   |                                                                                                                                                    
tGGCGTTCAGCGATTCGGCGACTTTCGCCGCTTCTGGCATCAGCGTACCAAAGTTAAGGATCGCCAGTTTCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCa                                                                                                                                                     >  2:1017565/1‑149 (MQ=255)
tGGCGTTCAGCGATTCGGCGACTTTCGCCGCTTCTGGCATCAGCGTACCAAAGTTAAGGATCGCCAGTTTCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCa                                                                                                                                                     >  2:1017369/1‑149 (MQ=255)
 ggCGTTCAGCGATTCGGCGACTTTCGCCGCTTCTGGCATCAGCGTACCAAAGTTAAGGATCGCCAGTTTCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCAc                                                                                                                                                    >  4:989330/1‑149 (MQ=255)
   cGTTCAGCGATTCGGCGACTTTCGCCGCTTCTGGCATCAGCGTACCAAAGTTAAGGATCGCCAGTTTCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCAcgc                                                                                                                                                  <  4:834886/149‑1 (MQ=255)
    gTTCAGCGATTCGGCGACTTTCGCCGCTTCTGGCATCAGCGTACCAAAGTTAAGGATCGCCAGTTTCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCAcgcg                                                                                                                                                 >  4:640242/1‑149 (MQ=255)
    gTTCAGCGATTCGGCGACTTTCGCCGCTTCTGGCATCAGCGTACCAAAGTTAAGGATCGCCAGTTTCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCAcgcg                                                                                                                                                 >  2:1067770/1‑149 (MQ=255)
    gTTCAGCGATTCGGCGACTTTCGCCGCTTCTGGCATCAGCGTACCAAAGTTAAGGATCGCCAGTTTCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCAcgcg                                                                                                                                                 >  4:265711/1‑149 (MQ=255)
          cGATTCGGCGACTTTCGCCGCTTCTGGCATCAGCGTACCAAAGTTAAGGATCGCCAGTTTCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAgc                                                                                                                                           <  3:972111/149‑1 (MQ=255)
           gATTCGGCGACTTTCGCCGCTTCTGGCATCAGCGTACCAAAGTTAAGGATCGCCAGTTTCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAgcg                                                                                                                                          <  3:723694/149‑1 (MQ=255)
                   gACTTTCGCCGCTTCTGGCATCAGCGTACCAAAGTTAAGGATCGCCAGTTTCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTg                                                                                                                                  >  2:721310/1‑149 (MQ=255)
                   gACTTTCGCCGCTTCTGGCATCAGCGTACCAAAGTTAAGGATCGCCAGTTTCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTg                                                                                                                                  >  4:1056285/1‑149 (MQ=255)
                    aCTTTCGCCGCTTCTGGCATCAGCGTACCAAAGTTAAGGATCGCCAGTTTCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGa                                                                                                                                 >  1:308645/1‑149 (MQ=255)
                      tttCGCCGCTTCTGGCATCAGCGTACCAAAGTTAAGGATCGCCAGTTTCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACg                                                                                                                               >  4:175932/1‑149 (MQ=255)
                      tttCGCCGCTTCTGGCATCAGCGTACCAAAGTTAAGGATCGCCAGTTTCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACg                                                                                                                               >  2:680967/1‑149 (MQ=255)
                      tttCGCCGCTTCTGGCATCAGCGTACCAAAGTTAAGGATCGCCAGTTTCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACg                                                                                                                               >  4:467670/1‑149 (MQ=255)
                       ttCGCCGCTTCTGGCATCAGCGTACCAAAGTTAAGGATCGCCAGTTTCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACgg                                                                                                                              <  2:576397/149‑1 (MQ=255)
                          gccgcTTCTGGCATCAGCGTACCAAAGTTAAGGATCGCCAGTTTCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGcc                                                                                                                           >  1:357875/1‑149 (MQ=255)
                          gccgcTTCTGGCATCAGCGTACCAAAGTTAAGGATCGCCAGTTTCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGcc                                                                                                                           >  1:14751/1‑149 (MQ=255)
                             gcTTCTGGCATCAGCGTACCAAAGTTAAGGATCGCCAGTTTCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATc                                                                                                                        <  2:108777/149‑1 (MQ=255)
                               ttCTGGCATCAGCGTACCAAAGTTAAGGATCGCCAGTTTCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGt                                                                                                                      <  4:205465/149‑1 (MQ=255)
                               ttCTGGCATCAGCGTACCAAAGTTAAGGATCGCCAGTTTCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGt                                                                                                                      <  2:36381/149‑1 (MQ=255)
                                tCTGGCATCAGCGTACCAAAGTTAAGGATCGCCAGTTTCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGtt                                                                                                                     <  4:212577/149‑1 (MQ=255)
                                   ggCATCAGCGTACCAAAGTTAAGGATCGCCAGTTTCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTtata                                                                                                                  <  3:265711/149‑1 (MQ=255)
                                      aTCAGCGTACCAAAGTTAAGGATCGCCAGTTTCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTg                                                                                                               >  4:983844/1‑149 (MQ=255)
                                      aTCAGCGTACCAAAGTTAAGGATCGCCAGTTTCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTg                                                                                                               >  1:1078255/1‑149 (MQ=255)
                                       tCAGCGTACCAAAGTTAAGGATCGCCAGTTTCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGa                                                                                                              <  4:1073760/149‑1 (MQ=255)
                                             tACCAAAGTTAAGGATCGCCAGTTTCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCg                                                                                                        <  1:1033113/149‑1 (MQ=255)
                                               ccAAAGTTAAGGATCGCCAGTTTCTCGCCACGACGCTTCACAAAGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGt                                                                                                      <  1:1017565/149‑1 (MQ=255)
                                                 aaaGTTAAGGATCGCCAGTTTCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGtat                                                                                                    >  3:271825/1‑149 (MQ=255)
                                                   aGTTAAGGATCGCCAGTTTCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATag                                                                                                  <  2:146698/149‑1 (MQ=255)
                                                    gTTAAGGATCGCCAGTTTCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATaga                                                                                                 <  3:527496/149‑1 (MQ=255)
                                                    gTTAAGGATCGCCAGTTTCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATaga                                                                                                 <  2:907174/149‑1 (MQ=255)
                                                     ttAAGGATCGCCAGTTTCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATagag                                                                                                <  1:478199/149‑1 (MQ=255)
                                                          gATCGCCAGTTTCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCt                                                                                           >  1:332755/1‑149 (MQ=255)
                                                             cGCCAGTTTCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGc                                                                                        <  3:989330/149‑1 (MQ=255)
                                                                       tcGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGtt                                                                              <  2:1059233/149‑1 (MQ=255)
                                                                        cGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGttt                                                                             <  2:357875/149‑1 (MQ=255)
                                                                           cACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTg                                                                                                               <  1:668701/112‑1 (MQ=255)
                                                                           cACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTg                                                                                                               >  2:668701/1‑112 (MQ=255)
                                                                                 gCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGtttt                                                                            <  2:183904/141‑1 (MQ=255)
                                                                                 gCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGtttt                                                                            >  1:183904/1‑141 (MQ=255)
                                                                                 gCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTc                                                                    <  4:271825/149‑1 (MQ=255)
                                                                                  cTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTGGCAACGCggtc                                                                                                                                              >  1:976804/1‑72 (MQ=255)
                                                                                  cTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTGGCAACGCggtc                                                                                                                                              <  2:976804/74‑3 (MQ=255)
                                                                                   ttCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCAt                                                                         >  3:283983/1‑142 (MQ=255)
                                                                                   ttCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCgg                                                                  >  3:233013/1‑149 (MQ=255)
                                                                                   ttCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCgg                                                                  <  3:490223/149‑1 (MQ=255)
                                                                                    tCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCggg                                                                 <  1:1067770/149‑1 (MQ=255)
                                                                                     cacaATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCgggg                                                                >  4:668792/1‑149 (MQ=255)
                                                                                      acaATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGt                                                               >  4:1040651/1‑149 (MQ=255)
                                                                                      acaATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGt                                                               >  3:347284/1‑149 (MQ=255)
                                                                                      acaATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCACTCGGGGt                                                               <  3:1056285/149‑1 (MQ=255)
                                                                                        aaTGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGTCa                                                             >  4:325153/1‑149 (MQ=255)
                                                                                          tGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGTCATa                                                           <  4:233013/149‑1 (MQ=255)
                                                                                             cTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGTCATAATg                                                        >  2:42222/1‑149 (MQ=255)
                                                                                              tttGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGTCATaa                                                          >  3:699399/1‑146 (MQ=255)
                                                                                              tttGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGTCATAATGa                                                       <  1:680967/149‑1 (MQ=255)
                                                                                              tttGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGTCATAATGa                                                       >  1:76875/1‑149 (MQ=255)
                                                                                               ttGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGTCATAATGAc                                                      <  3:640242/149‑1 (MQ=255)
                                                                                                  ccAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTtat                                                                                                                   >  4:359493/1‑85 (MQ=255)
                                                                                                  ccAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGTCATAATGACCAt                                                   >  3:965907/1‑149 (MQ=255)
                                                                                                    aaTTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGTCATAATGACCAttt                                                 >  4:182136/1‑149 (MQ=255)
                                                                                                     aTTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGTCATAATGACCATTTc                                                >  4:102735/1‑149 (MQ=255)
                                                                                                        ggTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGTCATAATGACCATTTCCgg                                             <  4:956464/149‑1 (MQ=255)
                                                                                                         gTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCa                                                                          <  3:375259/119‑1 (MQ=255)
                                                                                                         gTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCa                                                                          >  4:375259/1‑119 (MQ=255)
                                                                                                         gTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGTCATAATGACCATTTCCGGt                                            <  3:983844/149‑1 (MQ=255)
                                                                                                         gTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGTCATAATGACCATTTCCGGt                                            >  1:222802/1‑149 (MQ=255)
                                                                                                          tAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGTCATAATGACCATTTCCGGTa                                           <  3:175932/149‑1 (MQ=255)
                                                                                                          tAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGTCATAATGACCATTTCCGGTa                                           >  1:714485/1‑149 (MQ=255)
                                                                                                             tttttCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGa                                                                                                                                 <  4:1138294/60‑1 (MQ=255)
                                                                                                             tttttCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGa                                                                                                                                 >  3:1138294/1‑60 (MQ=255)
                                                                                                             tttttCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGTCATAATGACCATTTCCGGTATGc                                        >  2:150360/1‑149 (MQ=255)
                                                                                                             tttttCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGTCATAATGACCATTTCCGGTATGc                                        >  3:1125569/1‑149 (MQ=255)
                                                                                                                   cAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGTCATAATGACCATTTCCGGTATGCAGCGCa                                  <  4:347284/149‑1 (MQ=255)
                                                                                                                    aGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGTCATAATGACCATTTCCGGTATGCAGCGCAg                                 <  2:14751/149‑1 (MQ=255)
                                                                                                                    aGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGTCATAATGACCATTTCCGGTATGCAGCGCAg                                 <  1:391944/149‑1 (MQ=255)
                                                                                                                    aGCCGCGACAGTTCCACCCCGACCGCGTTGCCACGCGGGGAGCACACAGCTGACGGGCCATCGTTATAGTGATAGCTGGTATGGAGCATCAGACGAAATTCGCGTTCATCGATCGGGGCCATAATGACCCTTTCCGGTAGGCAGCGCAg                                 >  2:515476/1‑149 (MQ=255)
                                                                                                                     gcggcgTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGTCATAATGACCATTTCCGGTATGCAGCGCAgg                                >  3:646239/1‑149 (MQ=255)
                                                                                                                          gTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGTCATAATGACCATTTCCGGTATGCAGCGCAGGTAaga                           <  4:965907/149‑1 (MQ=255)
                                                                                                                          gTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGTCATAATGACCATTTCCGGTATGCAGCGCAGGTAaga                           <  2:205677/149‑1 (MQ=255)
                                                                                                                          gTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGTCATAATGACCATTTCCGGTATGCAGCGCAGGTAaga                           <  2:222802/149‑1 (MQ=255)
                                                                                                                           tCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGTCATAATGACCATTTCCGGTATGCAGCGCAGGTAagag                          >  1:876046/1‑149 (MQ=255)
                                                                                                                             aGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGTCATAATGACCATTTCCGGTATGCAGCGCAGGTAagagag                        >  3:310215/1‑149 (MQ=255)
                                                                                                                               ttCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGTCATAATGACCATTTCCGGTATGCAGCGCAGGTAAGAGAGAt                      >  4:702017/1‑149 (MQ=255)
                                                                                                                                   aCGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATaga                                                                                                 >  4:934623/1‑70 (MQ=255)
                                                                                                                                    cGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGTCATAATGACCATTTCCGGTATGCAGCGCAGGTAAGAGAGATCaaaa                 <  1:1017369/149‑1 (MQ=255)
                                                                                                                                      ctacccGCGTGGCAACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGTCATAATGAc                                                      <  2:75389/106‑1 (MQ=255)
                                                                                                                                       tacccGCGTGGCAACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGtt                                                                              >  4:257154/4‑85 (MQ=255)
                                                                                                                                       tacccGCGTGGCAACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGtt                                                                              <  3:257154/82‑1 (MQ=255)
                                                                                                                                       cGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCg                                                                                                        >  1:665850/1‑59 (MQ=255)
                                                                                                                                       cGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCg                                                                                                        <  2:665850/59‑1 (MQ=255)
                                                                                                                                        acccGCGTGGCAACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGTCATAATGACCATTTCCGGTATGCAGCGCAGGTAAGAGAGAt                      <  1:8239/138‑1 (MQ=255)
                                                                                                                                        acccGCGTGGCAACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGTCATAATGACCATTTCCGGTATGCAGCGCAGGTAAGAGAGAt                      >  2:8239/3‑140 (MQ=255)
                                                                                                                                         aCCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGTCATAATGACCATTTCCGGTATGCAGCGCAGGTAAGAGAGATCAAAAGCAcc            >  3:1108355/1‑149 (MQ=255)
                                                                                                                                              gTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGTCATAATGACCATTTCCGGTATGCAGCGCAGGTAAGAGAGATCAAAAGCACCCTGAt       <  4:646239/149‑1 (MQ=255)
                                                                                                                                               ttGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGTCATAATGACCATTTCCGGTATGCAGCGCAGGTAAGAGAGATCAAAAGCACCCTGATg      <  3:325153/149‑1 (MQ=255)
                                                                                                                                                 gCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGTCATAATGACCATTTCCGGTATGCAGCGCAGGTAAGAGAGATCAAAAGCACCCTGATggg    >  4:994982/1‑149 (MQ=255)
                                                                                                                                                  ccACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGTCATAATGACCATTTCCGGTATGCAGCGCAGGTAAGAGAGATCAAAAGCACCCTGATGGGt   <  3:467670/149‑1 (MQ=255)
                                                                                                                                                   cACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGTCATAATGACCATTTCCGGTATGCAGCGCAGGTAAGAGAGATCAAAAGCACCCTGATGGGtt  <  3:1040651/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TGGCGTTCAGCGATTCGGCGACTTTCGCCGCTTCTGGCATCAGCGTACCAAAGTTAAGGATCGCCAGTTTCTCGCCACGACGCTTCACAATGCCTTTGCCAATTGGTAGTTTTTCCAGCGGCGTCAGTTCCACGCCGACCGCGTTGCCACGCGGGTAGCGCACCGCTGACGGGCCATCGTTATAGTGATAGCCGGTATAGAGCATCTGGCGACATTCGTTTTCATCGCTCGGGGTCATAATGACCATTTCCGGTATGCAGCGCAGGTAAGAGAGATCAAAAGCACCCTGATGGGTT  >  NC_000913/438596‑438891

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