Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,562,183 G→T 5.1% E57D (GAG→GAT glpD → sn‑glycerol‑3‑phosphate dehydrogenase, aerobic, FAD/NAD(P)‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,562,1830GT 5.1% 270.0 / 3.3 78E57D (GAG→GATglpDsn‑glycerol‑3‑phosphate dehydrogenase, aerobic, FAD/NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base G (38/36);  new base T (2/2);  total (40/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.70e-01

TGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGTTTATCCGTGCTGATGCTGGAGGCGCAGGATCTCGCTTGCGCGACCTCTTCCGCCAGTTCAAAACTCATTCACGGTGGCCTGCGCTACCTTGAGCACTATGAATTCCGCCTGGTCAGCGAGGCGCTGGCTGAACGTGAAGTGCTGCTGAAAATGGCCCCGCATATCGCCTTCCCGATGCGTTTTCGCCTGCCACATCGTCCGCATCTGCGCCCGGCGTGGATGATTCGCATTGGTCTGTT  >  NC_000913/3562035‑3562329
                                                                                                                                                    |                                                                                                                                                  
tgATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGTTTATCCGTGCTGATGCTGGAGGCGCAGGATCTCGCTTGCGCGACCTCTTCCGCCAGTTCAAAACTCATTCACGGTGGCCTGCGCTACCTTGAg                                                                                                                                                    <  6:83471/149‑1 (MQ=255)
tgATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGTTTATCCGTGCTGATGCTGGAGGCGCAGGATCTCGCTTGCGCGACCTCTTCCGCCAGTTCAAAACTCATTCACGGTGGCCTGCGCTACCTTGAg                                                                                                                                                    <  6:43614/149‑1 (MQ=255)
   tAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGTTTATCCGTGCTGATGCTGGAGGCGCAGGATCTCGCTTGCGCGACCTCTTCCGCCAGTTCAAAACTCATTCACGGTGGCCTGCGCTACCTTGAGCAc                                                                                                                                                 <  5:105638/149‑1 (MQ=255)
   tAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGTTTATCCGTGCTGATGCTGGAGGCGCAGGATCTCGCTTGCGCGACCTCTTCCGCCAGTTCAAAACTCATTCACGGTGGCCTGCGCTACCTTGAGCAc                                                                                                                                                 <  2:327287/149‑1 (MQ=255)
     gggggCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGTTTATCCGTGCTGATGCTGGAGGCGCAGGATCTCGCTTGCGCGACCTCTTCCGCCAGTTCAAAACTCATTCACGGTGGCCTGCGCTACCTTGAGCACTa                                                                                                                                               <  5:322660/149‑1 (MQ=255)
      ggggCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGTTTATCCGTGCTGATGCTGGAGGCGCAGGATCTCGCTTGCGCGACCTCTTCCGCCAGTTCAAAACTCATTCACGGTGGCCTGCGCTACCTTGAGCACTAt                                                                                                                                              >  4:372615/1‑149 (MQ=255)
          cggcATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGTTTATCCGTGCTGATGCTGGAGGCGCAGGATCTCGCTTGCGCGACCTCTTCCGCCAGTTCAAAACTCATTCACGGTGGCCTGCGCTACCTTGAGCACTATGAAt                                                                                                                                          <  2:194350/149‑1 (MQ=255)
             cATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGTTTATCCGTGCTGATGCTGGAGGCGCAGGATCTCGCTTGCGCGCCCTCTTCCGCCAGTTCAAAACTCATTCACGGTGGCCTGCGCTACCTTGAGCACTATGAATTcc                                                                                                                                       >  3:129427/1‑149 (MQ=255)
              aTCAATGGTGCTGGTATCGCGGCAGCCGCCGCTGGACGCGGTTTATCCGTGCTGATGCTGGAGGCGCAGGATCTCGCTTGCGCGACCTCTTCCGCCAGTTCAAAACTCATTCACGGTGGCCTGCGCTACCTTGAGCACTATGAATTCCg                                                                                                                                      >  5:84897/1‑149 (MQ=255)
                  aTGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGTTTATCCGTGCTGATGCTGGAGGCGCAGGATCTCGCTTGCGCGACCTCTTCCGCCAGTTCAAAACTCATTCACGGTGGCCTGCGCTACCTTGAGCACTATGAATTCCGCCTg                                                                                                                                  <  4:108503/149‑1 (MQ=255)
                  aTGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGTTTATCCGTGCTGATGCTGGAGGCGCAGGATCTCGCTTGCGCGACCTCTTCCGCCAGTTCAAAACTCATTCACGGTGGCCTGCGCTACCTTGAGCACTATGAATTCCGCCTg                                                                                                                                  <  6:171011/149‑1 (MQ=255)
                      tGCTGGTATCGCGGCAGACGCCGCTGGACGCGGTTTATCCGTGCTGATGCTGGAGGCGCAGGATCTCGCTTGCGCGACCTCTTCCGCCAGTTCAAAACTCATTCACCGTGGCCTGCGCTACCTTGAGCACTATGAATTCCGCCTGGt                                                                                                                                >  6:231073/1‑147 (MQ=255)
                           gTATCGCGGCAGACGCCGCTGGACGCGGTTTATCCGTGCTGATGCTGGAGGCGCAGGATCTCGCTTGCGCGACCTCTTCCGCCAGTTCAAAACTCATTCACGGTGGCCTGCGCTACCTTGAGCACTATGAATTCCGCCTGGTCAGCgcg                                                                                                                         >  7:352375/1‑147 (MQ=255)
                              tCGCGGCAGACGCCGCTGGACGCGGTTTATCCGTGCTGATGCTGGAGGCGCAGGATCTCGCTTGCGCGACCTCTTCCGCCAGTTCAAAACTCATTCACGGTGGCCTGCGCTACCTTGAGCACTATGAATTCCGCCTGGTCAGCGAGgc                                                                                                                       <  7:79615/148‑1 (MQ=255)
                               agcgGCAGACGCCGCTGGACGCGGTTTATCCGTGCTGATGCTGGAGGCGCAGGATCTAGCTTGCGCGACCTCTTCCGCCAGTTCAAAACTCATTAACGGTGGCCTTAGATACCTTGAGCACTATGAATTCCGCATGGTCAGCGAGgcgc                                                                                                                     >  8:288151/2‑149 (MQ=255)
                                 cgGCAGACGCCGCTGGACGCGGTTTATCCGTGCTGATGCTGGAGGCGCAGGATCTCGCTTGCGCGACCTCTTCCGCCAGTTCAAAACTCATTCACGGTGGCCTGCGCTACCTTGAGCACTATGAATTCCGCCTGGTCAGCGAGGCgctg                                                                                                                   >  7:72607/1‑149 (MQ=255)
                                 cgGCAGACGCCGCTGGACGCGGTTTATCCGTGCTGATGCTGGAGGCGCAGAATCTCGCTTGCGCGACCTCTTCCGCCAGTTCAAAACTCATTCACGGTGGCCTGCGCTACCTTGATCACTATGAATTCCGCCTGGTCAGCGAGGCgctg                                                                                                                   >  3:335244/1‑149 (MQ=255)
                                 cgGCAGACGCCGCTGGACGCGGTTTATCCGTGCTGATGCTGGAGGCGCAGAATCTCGCTTGCGCGACCTCTTCCGCCAGTTCAAAACTCATTCACGGTGGCCTGCGCTACCTTGATCACTATGAATTCCGCCTGGTCAGCGAGGCgctg                                                                                                                   >  3:340602/1‑149 (MQ=255)
                                  ggCAGACGCCTCTGGACGCGGTTTATCCGTGCTGATGATGGAGGCGCAGGATCTCGCTTGCGCGACCTCTTCCGCCAGTTCAAAACTCATTCACGGTGGCCTGCGCTACCTTGAGCACTATGAATTCCGCCTGGTCAGCGAGGCtctgg                                                                                                                  <  7:288151/149‑1 (MQ=255)
                                   gCAGACGCCGCTGGACGCGGTTTATCCGTGCTGATGCTGGAGGCGCAGGATCTCGCTTGCGCGACCTCTTCCGCCAGTTCAAAACTCATTCACGGTGGCCTGCGCTACCTTGAGCACTATGAATTCCGCCTGGTCAGCGAGGCgctggc                                                                                                                 <  8:390125/149‑1 (MQ=255)
                                       aCGCCGCTGGACGCGGTTTATCCGTGCTGATGCTGGAGGCGCAGGATCTCGCTTGCGCGACCTCTTCCGCCAGTTCAAAACTCATTCACGGTGGCCTGCGCTACCTTGAGCACTATGAATTCCGCCTGGTCAGCGAGGCGCTGGCTGaa                                                                                                             >  8:37592/1‑149 (MQ=255)
                                        cgccgcTGGACGCGGTTTATCCGTGCTGATGCTGGAGGCGCAGGATCTCGCTTGCGCGACCTCTTCCGCCAGTTCAAAACTCATTCACGGTGGCCTGCGCTACCTTGAGCACTATGAATTCCGCCTGGTCAGCGAGGCGCTGGCTGAAc                                                                                                            <  8:85669/149‑1 (MQ=255)
                                            gcTGGACGCGGTTTATCCGTGCTGATGCTGGAGGCGCAGGATCTCGCTTGCGCGACCTCTTCCGCCAGTTCAAAACTCATTCACGGTGGCCTGCGCTACCTTGAGCACTATGAATTCCGCCTGGTCAGCGAGGCGCTGGCTGAACGTg                                                                                                         >  3:108109/1‑148 (MQ=255)
                                             cTGGACGCGGTTTATCCGTGCTGATGCTGGAGGCGCAGGATCTCGCTTGCGCGACCTCTTCCGCCAGTTCAAAACTCATTCACGGTGGCCTGCGCTACCTTGAGCACTATGAATTCCGCCTGGTCAGCGAGGCGCTGGCTGAACGTGaa                                                                                                       <  3:372615/149‑1 (MQ=255)
                                               ggACGCGGTTTATCCGTGCTGATGCTGGAGGCGCAGGATCTCGCTTGCGCGACCTCTTCCGCCAGTTCAAAACTCATTCACGGTGGCCTGCGCTACCTTGAGCACTATGAATTCCGCCTGGTCAGCGAGGCGCTGGCTGAACGTGAAGt                                                                                                     <  8:352375/149‑1 (MQ=255)
                                                gACGCGGTTTATCCGTGCTGATGCTGGAGGCGCAGGATCTCGCTTGCGCGACCTCTTCCGCCAGTTCAAAACTCATTCACGGTGGCCTGCGCTACCTTGAGCACTATGAATTCCGCCTGGTCAGCGAGGCGCTGGCTGAACGTGAAGtg                                                                                                    <  6:900/149‑1 (MQ=255)
                                                 aCGCGGTTTATCCGTGCTGATGCTGGAGGCGCAGGATCTCGCTTGCGCGACCGCTTCCGCCAGTTCAAAACTCATTCACGGTGGCCTGCGCTACCTTGAGCACTATGAATTCCGCCTGGTCAGCGAGGCGCTGGCTGAACGTGAAGtgc                                                                                                   <  6:374559/149‑1 (MQ=255)
                                                      gTTTATCCGTGCTGATGCTGGAGGCGCAGGATCTCGCTTGCGCGACCTCTTCCGCCAGTTCAAAACTCATTCACGGTGGCCTGCGCTACCTTGAGCACTATGAATTCCGCCTGGTCAGCGAGGCGCTGGCTGAACGTGAAGTgctgct                                                                                               >  8:362087/1‑148 (MQ=255)
                                                           tCCGTGCTGATGCTGGAGGCGCAGTATCTCGCTTGCGCGACCTCTTCCGCCAGTTCAAAACTCATTCACGGTGGCCTGCGCTACCTTGAGCACTATGAATTCCGCCTGGTCAGCGAGGCGCTGGCTGAACGTGAAGTGCTGCTGAAAAt                                                                                         <  6:84897/149‑1 (MQ=255)
                                                           tCCGTGCTGATGCTGGAGGCGCAGGATCTCGCTTGCGCGACCTCTTCCGCCAGTTCAAAACTCATTCACGGTGGCCTGCGCTACCTTGAGCACTATGAATTCCGCCTGGTCAGCGAGGCGCTGGCTGAACGTGAAGTGCTGCTGAAAAt                                                                                         >  5:64985/1‑149 (MQ=255)
                                                           tCCGTGCTGATGCTGGAGGCGCAGGATCTCGCTTGCGCGACCTCTTCCGCCAGTTCAAAACTCATTCACGGTGGCCTGCGCTACCTTGAGCACTATGAATTCCGCCTGGTCAGCGAGGCGCTGGCTGAACGTGAAGTGCTGCTGAAAAt                                                                                         <  8:106346/149‑1 (MQ=255)
                                                           tCCGTGCTGATGCTGGAGGCGCAGGATCTCGCTTGCGCGACCTCTTCCGCCAGTTCAAAACTCATTCACGGTGGCCTGCGCTACCTTGAGCACTATGAATTCCGCCTGGTCAGCGAGGCGCTGGCTGAACGTGAAGTGCTGCTGAAAAt                                                                                         <  8:113222/149‑1 (MQ=255)
                                                           tCCGTGCTGATGCTGGAGGCGCAGGATCTCGCTTGCGCGACCTCTTCCGCCAGTTCAAAACTCATTCACGGTGGCCTGCGCTACCTTGAGCACTATGAATTCCGCCTGGTCAGCGAGGCCCTGGCTGAACGTGAAGTGCTGCTGAAAAt                                                                                         >  6:108101/1‑149 (MQ=255)
                                                              gTGCTGATGCTGGAGGCGCAGGATCTCGCTTGCGCGACCTCTTCCGCCAGTTCAAAACTCATTCACGGTGGCCTGCGCTACCTTGAGCACTATGAATTCCGCCTGGTCAGCGAGGCGCTGGCTGAACGTGAAGTgctgct                                                                                               <  2:194485/140‑1 (MQ=255)
                                                              gTGCTGATGCTGGAGGCGCAGGATCTCGCTTGCGCGACCTCTTCCGCCAGTTCAAAACTCATTCACGGTGGCCTGCGCTACCTTGAGCACTATGAATTCCGCCTGGTCAGCGAGGCGCTGGCTGAACGTGAAGTgctgct                                                                                               >  1:194485/1‑140 (MQ=255)
                                                              gTGCTGATGCTGGAGGCGCAGGATCTCGCTTGCGCGACCTCTTCCGCCAGTTCAAAACTCATTCACGGTGGCCTGCGCTACCTTGAGCACTATGAATTCCGCCTGGTCAGCGAGGCGCTGGCTGAACGTGAAGTGCTGCTGAAAATGGc                                                                                      <  4:300731/149‑1 (MQ=255)
                                                              gTGCTGATGCTGGAGGCGCAGGATCTCGCTTGCGCGACCTCTTCCGCCAGTTCAAAACTCATTCACGGTGGCCTGCGCTACCTTGAGCACTATGAATTCCGCCTGGTCAGCGAGGCGCTGGCTGAACGTGAAGTGCTGCTGAAAATGGc                                                                                      >  3:236958/1‑149 (MQ=255)
                                                                gCTGATGCTGGAGGCGCAGGATCTCGCTTGCGCGACCTCTTCCGCCAGTTCAAAACTCATTCACGGTGGCCTGCGCTACCTTGAGCACTATGAATTCCGCCTGGTCAGCGAGGCGCTGGCTGAACGTGAAGTGCTGCTGAAAATGGccc                                                                                    >  8:103407/1‑149 (MQ=255)
                                                                  tGATGCTGGAGGCGCAGGATCTCGCTTGCGCGACCTCTTCCGCCAGTTCAAAACTCATTCACGGTGGCCTGCGCTACCTTGAGCACTATGAATTCCGCCTGGTCAGCGAGGCGCTGGCTGAACGTGAAGTGCTGCTGAAAATGGCCCCg                                                                                  >  6:175157/1‑149 (MQ=255)
                                                                   gATGCTGGAGGCGCAGGATCTCGCTTGCGCGACCTCTTCCGCCAGTTCAAAACTCATTCACGGTGGCCTGCGCTACCTTGAGCACTATGAATTCCGCCTGGTCAGCGAGGCGCTGGCTGAACGTGAAGTGCTGCTGAAAATGGCCCCg                                                                                  >  1:381715/1‑148 (MQ=255)
                                                                        tGGAGGCGCAGGATCTCGCTTGCGCGACCTCTTCCGCCAGTTCAAAACTCATTCACGGTGGCCTGCGCTACCTTGAGCACTATGAATTCCGCCTGGTCAGCGAGGCGCTGGCTGAACGTGAAGTGCTGCTGAAAATGGCCCCGCATATc                                                                            <  2:381715/149‑1 (MQ=255)
                                                                          gAGGCGCAGGATCTCGCTTGCGCGACCTCTTCCGCCAGTTCAAAACTCATTCACGGTGGCCTGCGCTACCTTGAGCACTCTGAATTCCGCCTGGTCAGCGAGGCGCTGGCTGAACGTGAAGTGCTGCTGAAAATGGCCCCGCATATCGc                                                                          <  7:136488/149‑1 (MQ=255)
                                                                            ggCGCAGGATCTCGCTTGCGCGACCTCTTCCGCCAGTTCAAAACTCATTCACGGTGGCCTGCGCTACCTTGAGCACTATGAATTCCGCCTGGTCAGCGAGGCGCTGGCTGAACGTGAAGTGCTGCTGAAAATGGCCCCGCATATCGCCt                                                                        <  7:37592/149‑1 (MQ=255)
                                                                               gcAGGATCTCGCTTGCGCGACCTCTTCCGCCAGTTCAAAACTCATTCACGGTGGCCTGCGCTACCTTGAGCACTATGAATTCCGCCTGGTCAGCGAGGCGCTGGCTGAACGTGAAGTGCTGCTGAAAATGGCCCCGCATATCGCCTTcc                                                                     >  8:178258/1‑149 (MQ=255)
                                                                               gcAGGATCTCGCTTGCGCGACCTCTTCCGCCAGTTCAAAACTCATTCACGGTGGCCTGCGCTACCTTGAGCACTATGAATTCCGCCTGGTCAGCGAGGCGCTGGCTGAACGTGAAGTGCTGCTGAAAATGGCCCCGCATATCGCCTTcc                                                                     >  1:159329/1‑149 (MQ=255)
                                                                                 aGGATCTCGCTTGCGCGACCTCTTCCGCCAGTTCAAAACTCATTCACGGTGGCCTGCGCTACCTTGAGCACTATGAATTCCGCCTGGTCAGCGAGGCGCTGGCTGAACGTGAAGTGCTGCTGAAAATGGCCCCGCATATCGCCTTccc                                                                    >  7:406154/1‑148 (MQ=255)
                                                                                   gATCTCGCTTGCGCGACCTCTTCCGCCAGTTCAAAACTCATTCACGGTGGCCTGCGCTACCTTGAGCACTATGAATTCCGCCTGGTCAGCGAGGCGCTGGCTGAACGTGAAGTGCTGCTGAAAATGGCCCCGCATATCGCCTTCCCGAt                                                                 >  3:316704/1‑149 (MQ=255)
                                                                                     tctcGCTTGCGCGACCTCTTCCGCCAGTTCAAAACTCATTCACGGTGGCCTGCGCTACCTTGAGCACTGTGAATTCCGCCTGGTCAGCGAGGCGCTGGCTGAACGTGAAGTGCTGCTGAAAATGGCCCCGCATATCGCCTTCCCGATGc                                                               >  1:171916/1‑149 (MQ=255)
                                                                                     tctcGCTTGCGCGACCTCTTCCGCCAGTTCAAAACTCATTCACGGTGGCCTGCGCTACCTTGAGCACTATGAATTCCGCCTGGTCAGCGAGGCGCTGGCTGAACGTGAAGTGCTGCTGAAAATGGCCCCGCATATCGCCTTCCCGATGc                                                               >  7:257282/1‑149 (MQ=255)
                                                                                     tctcGCTTGCGCGACCTCTTCCGCCAGTTCAAAACTCATTCACGGTGGCCTGCGCTACCTTGAGCACTATGAATTCCGCCTGGTCAGCGAGGCGCTGGCTGAACGTGAAGTGCTGCTGAAAATGGCCCCGCATATCGCCTTCCCGATGc                                                               >  3:244807/1‑149 (MQ=255)
                                                                                         gCTTGCGCGACCTCTTCCGCCAGTTCAAAACTCATTCACGGTGGCCTGCGCTACCTTGAGCACTATGAATTCCGCCTGGTCAGCGAGGCGCTGGCTGAACGTGAAGTGCTTCTGAAAATGGCCCCGCATATCGCCTTCCCGATGCGttt                                                           <  1:395435/149‑1 (MQ=255)
                                                                                         gCTTGCGCGACCTCTTCCGCCAGTTCAAAACTCATTCACGGTGGCCTGCGCTACCTTGAGCACTATGAATTCCGCCTGGTCAGCGAGGCGCTGGCTGAACGTGAAGTGCTGCTGAAAATGGCCCCGCATATCGCCTTCCCGATGCGttt                                                           >  7:179011/1‑149 (MQ=255)
                                                                                           ttGCGCGACCTCTTCCGCCAGTTCAAAACTCATTCACGGTGGCCTGCGCTACCTTGAGCACTATGAATTCCGCCTGGTCAGCGAGGCGCTGGCTGAACGTGAAGTGCTGCTGAAAATGGCCCCGCATATCGCCTTCCCGATGCGTTTTc                                                         >  7:111219/1‑149 (MQ=255)
                                                                                             gcgcGACCTCTTCCGCCAGTTCAAAACTCATTCACGGTGGCCTGCGCTACCTTGAGCACTATGAATTCCGCCTGGTCAGCGAGGCGCTGGCTGAACGTGAAGTGCTGCTGAAAATGGCCCCGCATATCGCCTTCCCGATGCGTTTTc                                                         <  5:175157/147‑1 (MQ=255)
                                                                                               gcgACCTCTTCCGCCAGTTCAAAACTCATTCACGGTGGCCTGCGCTACCTTGAGCACTATGAATTCCGCCTGGTCAGCGAGGCgct                                                                                                                    >  3:166929/1‑86 (MQ=255)
                                                                                               gcgACCTCTTCCGCCAGTTCAAAACTCATTCACGGTGGCCTGCGCTACCTTGAGCACTATGAATTCCGCCTGGTCAGCGAGGCgct                                                                                                                    <  4:166929/86‑1 (MQ=255)
                                                                                                 gACCTCTTCCGCCAGTTCAAAACTAATTCACGGTGGCCTGCGCTACCTTGAGCACTATGAATTCCGCCTGGGCAGCGAGGCGCTGGCTGAACGTGAAGTGCTGCTGAAAATTGCCCCGCATATCGCCTTCCCGATTCGTTTTCGCCTGc                                                   <  1:230873/149‑1 (MQ=255)
                                                                                                  aCCTCTTCCGCCAGTTCAAAACTCATTCACGGTGGCCTGCGCTACCTTGAGCACTATGAATTCCGCCTGGTCAGCGAGGCGCTGGCTGAACGTGAAGTGCTGCTGAAAATGGCCCCGCATATCGCCTTCCCGATGCGTTTTCGCCTGc                                                   >  1:62478/1‑148 (MQ=255)
                                                                                                       ttCCGCCAGTTCAAAACTCATTCACGGTGGCCTGCGCTACCTTGAGCACTATGAATTCCGCCTGGTCAGCGAGGCGCTGGCTGAACGTGAAGTGCTGCTGAAAATGGCCCCGCATATCGCCTTCCCGATGCGTTTTCGCCTGCcaca                                               <  5:143097/147‑1 (MQ=255)
                                                                                                       ttCCGCCAGTTCAAAACTCATTCACGGTGGCCTGCGCTACCTTGAGCACTATGAATTCCGCCTGGTCAGCGAGGCGCTGGCTGAACGTGAAGTGCTGCTGAAAATGGCCCCGCATATCGCCTTCCCGATGCGTTTTCGCCTGCcaca                                               <  4:236958/147‑1 (MQ=255)
                                                                                                           gCCAGTTCAAAACTCATTCACGGTGGCCTGCGCTACCTTGATCACTATGAATTCCGCCTGGTCAGCGAGGCGCTGGCTGAACGTGAAGTGCTGCTGAAAATGGCCCCGCATATCGCCTTCCCGATGCGTTTTCGCCTGCCACATCGTcc                                         <  4:335244/149‑1 (MQ=255)
                                                                                                           gCCAGTTCAAAACTCATTCACGGTGGCCTGCGCTACCTTGATCACTATGAATTCCGCCTGGTCAGCGAGGCGCTGGCTGAACGTGAAGTGCTGCTGAAAATGGCCCCGCATATCGCCTTCCCGATGCGTTTTCGCCTGCCACATCGTcc                                         <  4:340602/149‑1 (MQ=255)
                                                                                                                  cAAAACTCATTCACGGTGGCCTGCGCTACCTTGAGCACTATGAATTCCGCCTGGTCAGCGAGGCGCTGGCTGAACGTGAAGTGCTGCTGAAAATGGCCCCGCATATCGCCTTCCCGATGCGTTTTCGCCTGCCACATCGTCCGCATCTg                                  >  8:92003/1‑149 (MQ=255)
                                                                                                                    aaaCTCATTCACGGTGGCCTGCGCTACCTTGAGCACTATGAATTCCGCCTGGTCAGCGAGGCGCTGGCTGAACGTGAAGTGCTGCTGAAAATGGCCCCGCATATCGCCTTCCCGATGCGTTTTCGCCTGCCACATCGTCCGCATCTgcg                                <  8:111219/149‑1 (MQ=255)
                                                                                                                    aaaCTCATTCACGGTGGCCTGCGCTACCTTGAGCACTATGAATTCCGCCTGGTCAGCGAGGCGCTGGCTGAACGTGAAGTGCTGCTGAAAATGGCCCCGCATATCGCCTTCCCGATGCGTTTTCGCCTGCCACATCGTCCGCATCTgcg                                >  8:113729/1‑149 (MQ=255)
                                                                                                                     aaCTCATTCACGGTGGCCTGCGCTACCTTGAGCACTATGAATTCCGCCTGGTCAGCGAGGCGCTGGCTGAACGTGAAGTGCTGCTGAAAATGGCCCCGCATATCGCCTTCCCGATGCGTTTTCGCCTGCCACATCGTCCGCATCTgcgc                               >  8:27362/1‑149 (MQ=255)
                                                                                                                          aTTCACGGTGGCCTGCGCTACCTTGAGCACTATGAATTCCGCCTGGTCAGCGAGGCGCTGGCTGAACGTGAAGTGCTGCTGAAAATGGCCCCGCATATCGCCTTCCCGATGCGTTTTCGCCTGCCACATCGTCCGCATCTGCGCCCgg                           <  8:257282/148‑1 (MQ=255)
                                                                                                                            tCACGGTGGCCTGCGCTACCTTGAGCACTATGAATTCCGCCTGGTCAGCGAGGCGCTGGCTGAACGTGAAGTGCTGCTGAAAATGGCCCCGCATATCGCCTTCCCGATGCGTTTTCGCCTGCCACATCGTCCGCATCTGCGCCCGGCGt                        >  8:93070/1‑149 (MQ=255)
                                                                                                                              aCGGTGGCCTGCGCTACCTTGAGCACTATGAATTCCGCCTGGTCAGCGAGGCGCTGGCTGAACGTGAAGTGCTGCTGAAAATGGCCCCGCATATCGCCTTCCCGATGCGTTTTCGCCTGCCACATCGTCCGCATCTGCGCCCGGCGTg                       >  3:244955/1‑148 (MQ=255)
                                                                                                                              aCGGTGGCCTGCGCTACCTTGAGCACTATGAATTCCGCCTGGTCAGCGAGGCGCTGGCTGAACGTGAAGTGCTGCTGAAAATGGCCCCGCATATCGCCTTCCCGATGCGTTTTCGCCTGCCACATCGTCCGCATCTGCGCCCGGCGTg                       >  7:286441/1‑148 (MQ=255)
                                                                                                                               cGGTGGCCTGCGCTACCTTGAGCACTATGAATTCCGCCTGGTCAGCGAGGCGCTGGCTGAACGTGAAGTGCTGCTGAAAATGGCCCCGCATATCGCCTTCCCGATGCGTTTTCGCCTGCCACATCGTCCGCATCTGCGCCCGGCGTGga                     <  4:108109/149‑1 (MQ=255)
                                                                                                                               cGGTGGCCTGCGCTACCTTGAGCACTATGAATTCCGCCTGGTCAGCGAGGCGCTGGCTGAACGTGAAGTGCTGCTGAAAATGGCCCCGCATATCGCCTTCCCGATGCGTTTTCGCCTGCCACATCGTCCGCATCTGCGCCCGGCGTGga                     <  7:27362/149‑1 (MQ=255)
                                                                                                                                     ccTGCGCTACCTTGAGCACTATGAATTCCGCCTGGTCAGCGAGGCGCTGGCTGAACGTGAAGTGCTGCTGAAAATGGCCCCGCATATCGCCTTCCCGATGCGTTTTCGCCTGCCACATCGTCCGCATCTGCGCCCGGCGTGGATGATTc               >  8:371242/1‑149 (MQ=255)
                                                                                                                                         cgcTACCTTGAGCACTATGAATTCCGCCTGGTCAGCGAGGCGCTGGCTGAACGTGAAGTGCTGCTGAAAATGGCCCCGCATATCGCCTTCCCGATGCGTTTTCGCCTGCCACATCGTCCGCAGCTGCGCCCGGCGTGGATGATTCGCa            <  5:231073/148‑1 (MQ=255)
                                                                                                                                            tACCTTGAGCACTATGAATTCCGCCTGGTCAGCGAGGCGCTGGCTGAACGTGAAGTGCTGCTGAAAATGGCCCCGCATATCGCCTTCCCGATGCGTTTTCGCCTGCCACATCGTCCGCATCTGCGCCCGGCGTGGATGATTCGCa            <  2:281490/145‑1 (MQ=255)
                                                                                                                                            tACCTTGAGCACTATGAATTCCGCCTGGTCAGCGAGGCGCTGGCTGAACGTGAAGTGCTGCTGAAAATGGCCCCGCATATCGCCTTCCCGATGCGTTTTCGCCTGCCACATCGTCCGCATCTGCGCCCGGCGTGGATGATTCGCa            >  1:281490/1‑145 (MQ=255)
                                                                                                                                              ccTTGAGCACTATGAATTCCGCCTGGTCAGCGAGGCGCTGGCTGAACGTGAAGTGCTGCTGAAAATGGCCCCGCATATCGCCTTCCCGATGCGTTTTCGCCTGCCACATCGTCCGCATCTGCGCCCGGCGTGGATGATTCGCATTGGTc      >  7:302162/1‑149 (MQ=255)
                                                                                                                                                ttGAGCACTATGAATTCCGCCTGGTCAGCGAGGCGCTGGCTGAACGTGAAGTGCTGCTGAAAATGGCCCCGCATATCGCCTTCCCGATGCGTTTTCGCCTGCCACATCGTCCGCATCTGCGCCCGGCGTGGATGATTCGCATTGGTCTg    >  4:283492/1‑149 (MQ=255)
                                                                                                                                                  gAGCACTATGAATTCCGCCTGGTCAGCGAGGCGCTGGCTGAACGTGAAGTGCTGCTGAAAATGGCCCCGCATATCGCCTTCCCGATGCGTTTTCGCCTGCCACATCGTCCGCATCTGCGCCCGGCGTGGATGATTCGCATTGGTCTGt   <  2:159329/148‑1 (MQ=255)
                                                                                                                                                   aGCACTATGAATTCCGCCTGGTCAGCGAGGCGCTGGCTGAACGTGAAGTGCTGCTGAAAATGGCCCCGCATATCGCCTTCCCGATGCGTTTTCGCCTGCCACATCGTCCGCATCTGCGCCCGGTGTGGATGATTCGCATTGGTCTGtt  <  4:244807/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGTTTATCCGTGCTGATGCTGGAGGCGCAGGATCTCGCTTGCGCGACCTCTTCCGCCAGTTCAAAACTCATTCACGGTGGCCTGCGCTACCTTGAGCACTATGAATTCCGCCTGGTCAGCGAGGCGCTGGCTGAACGTGAAGTGCTGCTGAAAATGGCCCCGCATATCGCCTTCCCGATGCGTTTTCGCCTGCCACATCGTCCGCATCTGCGCCCGGCGTGGATGATTCGCATTGGTCTGTT  >  NC_000913/3562035‑3562329

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