Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,659,730 G→C 5.3% E30Q (GAG→CAG)  nuoF → NADH‑quinone oxidoreductase subunit F

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,659,7300GC 5.3% 246.0 / 3.4 76E30Q (GAG→CAG) nuoFNADH‑quinone oxidoreductase subunit F
Reads supporting (aligned to +/- strand):  ref base G (38/34);  new base C (2/2);  total (40/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.51e-01

CGACACCTTCGGTGACGTGCAGCCGGCTGGCGTTTCCAAACTGCTGGAGGGTTACGTATGACCATTACTTCCTTCGGCCCGGCCAACCGCATCGCGCGCTCGGCCGAAACCCACCCGCTGACCTGGCGCCTGCGTGACGACGGCGAGCCGGTCTGGCTGGCCGAGTACGAATCGAAGAACGGCTACGCCGCAGCCCGCAAGGCGCTGGCGCAAATGTCCGCCGACGACATCGTGCAGAGCGTCAAGGACTCCGGCCTCAAAGGCCGTGGCGGTGCAGGTTTCCCCACTGGCGT  >  NC_002947/4659586‑4659878
                                                                                                                                                |                                                                                                                                                    
cgacACCTTCGGTGACGTGCAGCCGGCTGGCGTTTCCAAACTGCTGGAGGGTTACGTATGACCATTACTTCCTTCGGCCCGGCCAACCGCATCGCGCGCTCGGCCGAAACCCACCCGCTGACCTGGCGCCTGCGTGACGACGGCGAGcc                                                                                                                                                  >  1:161856/1‑149 (MQ=255)
 gacACCTTCGGTGACGTGCAGCCGGCTGGCGTTTCCAAACTGCTGGAGGGTTACGTATGACCATTACTTCCTTCGGCCCGGCCAACCGCATCGCGCGCTCGGCCGAAACCCACCCGCTGACCTGGCGCCTGCGTGACGACGGCGAGCCg                                                                                                                                                 >  7:89384/1‑149 (MQ=255)
 gacACCTTCGGTGACGTGCAGCCGGCTGGCGTTTCCAAACTGCTGGAGGGTTACGTATGACCATTACTTCCTTCGGCCCGGCCAACCGCATCGCGCGCTCGGCCGAAACCCACCCGCTGACCTGGCGCCTGCGTGACGACGGCGAGCCg                                                                                                                                                 >  2:153293/1‑149 (MQ=255)
    acCTTCGGTGACGTGCAGCCGGCTGGCGTTTCCAAACTGCTGGAGGGTTACGTATGACCATTACTTCCTTCGGCCCGGCCAACCGCATCGCGCGCTCGGCCGAAACCCACCCGCTGACCTGGCGCCTGCGTGACGACGGCGAGCCGGTc                                                                                                                                              >  5:219392/1‑149 (MQ=255)
      cTTCGGTGACGTGCAGCCGGCTGGCGTTTCCAAACTGCTGGAGGGTTACGTATGACCATTACTTCCTTCGGCCCGGCCAACCGCATCGCGCGCTCGGCCGAAACCCACCCGCTGACCTGGCGCCTGCGTGACGACGGCGAGCCGGTctg                                                                                                                                            >  3:320573/1‑149 (MQ=255)
      cTTCGGTGACGTGCAGCCGGCTGGCGTTTCCAAACTGCTGGAGGGTTACGTATGACCATTACTTCCTTCGGCCCGGCCAACCGCATCGCGCGCTCGGCCGAAACCCACCCGCTGACCTGGCGCCTGCGTGACGACGGCGAGCCGGTctg                                                                                                                                            >  6:63128/1‑149 (MQ=255)
       ttCGGTGACGTGCAGCCGGCTGGCGTTTCCAAACTGCTGGAGGGTTACGTATGACCATTACTTCCTTCGGCCCGGCCAACCGCATCGCGCGCTCGGCCGAAACCCACCCGCTGACCTGGCGCCTGCGTGACGACGGCGAGCCGGTc                                                                                                                                              <  1:125754/146‑1 (MQ=255)
       ttCGGTGACGTGCAGCCGGCTGGCGTTTCCAAACTGCTGGAGGGTTACGTATGACCATTACTTCCTTCGGCCCGGCCAACCGCATCGCGCGCTCGGCCGAAACCCACCCGCTGACCTGGCGCCTGCGTGACGACGGCGAGCCGGTc                                                                                                                                              >  2:125754/1‑146 (MQ=255)
               cGTGCAGCCGGCTGGCGTTTCCAAACTGCTGGAGGGTTACGTATGACCATTACTTCCTTCGGCCCGGCCAACCGCATCGCGCGCTCGGCCGAAACCCACCCGCTGACCTGGCGCCTGCGTGACGACGGCGAGCCGGTCTGGCTGGCCGa                                                                                                                                   <  1:325221/149‑1 (MQ=255)
               cGTGCAGCCGGCTGGCGTTTCCAAACTGCTGGAGGGTTACGTATGACCATTACTTCCTTCGGCCCGGCCAACCGCATCGCGCGCTCGGCCGAAACCCACCCGCTGACCTGGCGCCTGCGTGACGACGGCGAGCCGGTCTGGCTGGCCGa                                                                                                                                   <  1:153293/149‑1 (MQ=255)
                gTGCAGCCGGCTGGCGTTTCCAAACTGCTGGAGGGTTACGTATGACCATTACTTCCTTCGGCCCGGCCAACCGCATCGCGCGCTCGGCCGAAACCCACCCGCTGACCTGGCGCCTGCGTGACGACGGCGAGCCGGTCTGGCTGGCCGAg                                                                                                                                  <  5:63128/149‑1 (MQ=255)
                  gCAGCCGGCTGGCGTTTCCAAACTGCTGGAGGGTTACGTATGACCATTACTTCCTTCGGCCCGGCCAACCGCATCGCGCGCTCGGCCGAAACCCACCCGCTGACCTGGCGCCTGCGTGACGACGGCGAGCCGGTCTGGCTGGCCGAGTa                                                                                                                                >  2:90153/1‑149 (MQ=255)
                       cGGCTGGCGTTTCCAAACTGCTGGAGGGTTACGTATGACCATTACTTCCTTCGGCCCGGCCAACCGCATCGCGCGCTCGGCCGAAACCCACCCGCTGACCTGGCGCCTGCGTGACGACGGCGAGCCGGTCTGGCTGGCCGAGTACGAAt                                                                                                                           <  7:130479/149‑1 (MQ=255)
                           tGGCGTTTCCAAACTGCTGGAGGGTTACGTATGACCATTACTTCCTTCGGCCCGGCCAACCGCATCGCGCGCTCGGCCGAAACCCACCCGCTGACCTGGCGCCTGCGTGACGACGGCGAGCCGGTCTGGCTGGCCGAGTACGAATCg                                                                                                                         >  7:233403/1‑147 (MQ=255)
                             gCGTTTCCAAACTGCTGGAGGGTTACGTATGACCATTACTTCCTTCGGCCCGGCCAACCGCATCGCGCGCTCGGCCGAAACCCACCCGCTGACCTGGCGCCTGCGTGACGACGGCGAGCCGGTCTGGCTGGCCGAGTACGAATCgaag                                                                                                                      >  7:229938/1‑148 (MQ=255)
                              cGTTTCCAAACTGCTGGAGGGTTACGTATGACCATTACTTCCTTCGGCCCGGCCAACCGCATCGCGCGCTCGGCCGAAACCCACCCGCTGACCTGGCGCCTGCGTGACGACGGCGAGCCGGTctggctgg                                                                                                                                       >  7:336992/1‑130 (MQ=255)
                              cGTTTCCAAACTGCTGGAGGGTTACGTATGACCATTACTTCCTTCGGCCCGGCCAACCGCATCGCGCGCTCGGCCGAAACCCACCCGCTGACCTGGCGCCTGCGTGACGACGGCGAGCCGGTctggctgg                                                                                                                                       <  8:336992/130‑1 (MQ=255)
                                  tCCAAACTGCTGGAGGGTTACGTATGACCATTACTTCCTTCGGCCCGGCCAACCGCATCGCGCGCTCGGCCGAAACCCACCCGCTGACCTGGCGCCTGCGTGACGACGGCGAGCCGGTCTGGCTGGCCGAGTACGAATCGAAGAACgg                                                                                                                 <  2:161856/148‑1 (MQ=255)
                                  tCCAAACTGCTGGAGGGTTACGTATGACCATTACTTCCTTCGGCCCGGCCAACCGCATCGCGCGCTCGGCCGAAACCCACCCGCTGACCTGGCGCCTGCGTGACGACGGCGAGCCGGTCTGGCTGGCCGAGTACGAATCGAAGAACgg                                                                                                                 <  4:77742/148‑1 (MQ=255)
                                     aaaCTGCTGGAGGGTTACGTATGACCATTACTTCCTTCGGCCCGGCCAACCGCATCGCGCGCACGTACGAAACCCACCCCATGACCTGGAGCCTGCGTGACGACGGCGAGCCGGTCTGGATGCCCGAGTACGAATCGAAGAACGGCTAc                                                                                                             >  7:50267/1‑149 (MQ=255)
                                       aCTGCTGGAGGGTTACGTATGACCATTACTTCCTTCGGCCCGGCCAACCGCATCGCGCGCTCGGCCGAAACCCACCCGCTGACCTGGCGCCTGCGTGACGACGGCGAGCCGGTCTGGCTGGCCGAGTACGAATCGAAGAACGGCTAcgc                                                                                                           <  8:89384/149‑1 (MQ=255)
                                           ctgGAGGGTTACGTATGACCATTACTTCCTTCGGCCCGGCCAACCGCATCGCGCGCTCGGCCGAAACCCACCCGCTGACCTGGCGCCTGCGTGACGACGGCGAGCCGGTCTGGCTGGCCGAGTACGAATCGAAGAACGGCTAcgccg                                                                                                         >  2:116449/1‑147 (MQ=255)
                                             ggAGGGTTACGTATGACCATTACTTCCTTCGGCCCGGCCAACCGCATCGCGCGCTCGGCCGAAACCCACCCGCTGACCTGGCGCCTGCGTGACGACGGCCAGCCAGACCGGCTCGCCgtcgtc                                                                                                                               <  1:131238/123‑6 (MQ=255)
                                             ggAGGGTTACGTATGACCATTACTTCCTTCGGCCCGGCCAACCGCATCGCGCGCTCGGCCGAAACCCACCCGCTGACCTGGCGCCTGCGTGACGACGGCCAGCCAGACCGGCTCGCCgtcgtc                                                                                                                               >  2:131238/1‑118 (MQ=255)
                                              gAGGGTTACGTATGACCATTACTTCCTTCGGCCCGGCCAACCGCATCGCGCGCTCGGCCGAAACCCACCCGCTGACCTGGCGCCTGCGTGACGACGGCGAGCCGGTCTGGCTGGCCGAGTACGAATCGAAGAACGGCTACGCCGCAGcc                                                                                                    >  4:30936/1‑149 (MQ=255)
                                                gggTTACGTATGACCATTACTTCCTTCGGCCCGGCCAACCGCATCGCGCGCTCGGCCGAAACCCACCCGCTGACCTGGCGCCTGCGTGACGACGGCGAGCCGGTCTGGCTGGCCGAGTACGAATCGAAGAACGGCTACGCCGCAGCCCg                                                                                                  <  1:90153/149‑1 (MQ=255)
                                                 ggTTACGTATGACCATTACTTCTTTCGGCCCGGCCAACCGCATCGCGCGCTCGGCCGAAACCCACCCGCTGACCTGGCGCCTGCGTGACGACGGCGAGCCGGTCTGGCTGGCCGAGTACGAATCGAAGAACGGCTACGCCGCAGCCCGc                                                                                                 <  3:329538/149‑1 (MQ=255)
                                                   ttACGTATGACCATTACTTCCTTCGGCCCGGCCAACCGCATCGCGCGCTCGGCCGAAACCCACCCGCTGACCTGGCGCCTGCGTGACGACGGCGAGCCGGTCTGGCTGGCCGAGTACGAATCGAAGAACGGCTACGCCGCAGCCCGc                                                                                                 >  1:176261/1‑147 (MQ=255)
                                                   ttACGTATGACCATTACTTCCTTCGGCCCGGCCAACCGCATCGCGCGCTCGGCCGAAACCCACCCGCTGACCTGGCGCCTGCGTGACGACGGCGAGCCGGTCTGGCTGGCCGAGTACGAATCGAAGAACGGCTACGCCGCAGCCCGc                                                                                                 >  5:71769/1‑147 (MQ=255)
                                                    tACGTATGACCATTACTTCCTTCGGCCCGGCCAACCGCATCGCGCGCTCGGCCGAACCCCACCCACTGACCTCGCGCCTGCGTGACGACCGCGAGCCGGTCTGGCTGGCCGAGTACGAATCGAAGAACGGCTACGCCGCAGCCCGCAAg                                                                                              >  7:4924/1‑149 (MQ=255)
                                                       gTATGACCATTACTTCCTTCGGCCCGGCCAACCGCATCGCGCGCTCGGCCGAAACCCACCCGCTGACCTGGCGCCTGCGTGACGACGGCGAGCCGGTCTGGCTGGCCGAGTACGAATCGAAGAACGGCTACGCCGCAGCCCGCAAGgcg                                                                                           >  8:231897/1‑149 (MQ=255)
                                                          tGACCATTACTTCCTTCGGCCCGGCCAACCGCATCGCGCGCTCGGCCGAAACCCACCCGCTGACCTGGCGCCTGCGTGACGACGGCGAGCCGGTCTGGCTGGCCGAGTACGAATCGAAGAACGGCTACGCCGCAGCCCGCAAGGCGCTg                                                                                        <  1:116449/149‑1 (MQ=255)
                                                             ccATTACTTCCTTCGGCCCGGCCAACCGCATCGCGCGCTCGGCCGAAACCCACCCGCTGACCTGGCGCCTGCGTGACGACGGCGAGCCGGTCTGGCTGGCCGAGTACGAATCGAAGAACGGCTACGCCGCAGCCCGCAAGGCGCTGgcg                                                                                     >  6:22487/1‑149 (MQ=255)
                                                              cATTACTTCCTTCGGCCCGGCCAACCGCATCGCGCGCTCGGCCGAAACCCACCCGCTGACCTGGCGCCTGCGTGACGACGGCGAGCCGGTCTGGCTGGCCGAGTACGAATCGAAGAACGGCTACGCCGCAGCCCGCAAGGCGCTGgcgc                                                                                    <  4:320573/149‑1 (MQ=255)
                                                                  aCTTCCTTCGGCCCGGCCAACCGCATCGCGCGCTCGGCCGAAACCCACCCGCTGACCTGGCGCCTGCGTGACGACGGCGAGCCGGTCTGGCTGGCCGAGTACGAATCGAAGAACgg                                                                                                                 <  3:99680/116‑1 (MQ=255)
                                                                  aCTTCCTTCGGCCCGGCCAACCGCATCGCGCGCTCGGCCGAAACCCACCCGCTGACCTGGCGCCTGCGTGACGACGGCGAGCCGGTCTGGCTGGCCGAGTACGAATCGAAGAACgg                                                                                                                 >  4:99680/1‑116 (MQ=255)
                                                                  aCTTCCTTCGGCCCGGCCAACCGCATCGCGCGCTCGGCCGAAACCCACCCGCTGACCTGGCGCCTGCGTGACGACGGCGAGCCGGTCTGGCTGGCCGAGTACGAATCGAAGAACTGCTACGCCGCAGCCCGCAAGTCGCTGGCGCaaa                                                                                 <  8:50267/148‑1 (MQ=255)
                                                                    ttccttcGGCCCGGCCAACCGCATCGCGCGCTCGGCCGAAACCCACCCGCTGACCTGGCGCCTGCGTGACGACGGCGAGCCGGTCTGGCTGGCCGAGTACGAATCgaa                                                                                                                       >  7:212604/1‑108 (MQ=255)
                                                                    ttccttcGGCCCGGCCAACCGCATCGCGCGCTCGGCCGAAACCCACCCGCTGACCTGGCGCCTGCGTGACGACGGCGAGCCGGTCTGGCTGGCCGAGTACGAATCgaa                                                                                                                       <  8:212604/108‑1 (MQ=255)
                                                                       cttcGGCCCGGCCAACCGCATCGCGCGCTCGGCCGAAACCCACCCGCTGACCTGGCGCCTGCGTGACGACGGCGAGCCGGTCTGGCTGGCCGAGTACGAATCGAAGAACGTCTACGCCGCAGCCCGCAAGGCGCTGGCGCAAATGTccg                                                                           <  6:219392/149‑1 (MQ=255)
                                                                       cttcGGCCCGGCCAACCGCATCGCGCGCTCGGCCGAAACCCACCCGCTGACCTGGCGCCTGCGTGACGACGGCGAGCCGGTCTGGCTGGCCGAGTACGAATCGAAGAACGGCTACGCCGCAGCCCGCAAGGCGCTGGCGCAAATGTcc                                                                            <  6:15676/148‑1 (MQ=255)
                                                                          cggcccggccAACCGCATCGCGCGCTCGGCCGAAACCCACCCGCTGACCTGGCGCCTGCGTGACGACGGCGAGCCGGTCTGGCTGGCCGAGTACGAATCGAAGAACGGCTACGCCGCAGCCCGCAAGGCGCTGGCGCAAATGTccgccg                                                                        >  7:61749/1‑149 (MQ=255)
                                                                             cccggccAACCGCATCGCGCGCTCGGCCGAAACCCACCCGCTGACCTGGCGCCTGCGTGACGACGGCGAGCCGGTCTGGCTGGCCGAGTACGAATc                                                                                                                          >  6:1527/1‑96 (MQ=255)
                                                                             cccggccAACCGCATCGCGCGCTCGGCCGAAACCCACCCGCTGACCTGGCGCCTGCGTGACGACGGCGAGCCGGTCTGGCTGGCCGAGTACGAATc                                                                                                                          <  5:1527/96‑1 (MQ=255)
                                                                              ccggccAACCGCATCGCGCGCTCGGCCGAAACCCACCCGCTGACCTGGCGCCTGCGTGACGACGGCGAGCCGGTCTGGCTGGCCGAGTACGAATCGAAGAACGGCTACGCCGCAGCCCGCAAGGCGCTGGCGCAAATGTCCGCcgacga                                                                    >  7:114432/1‑149 (MQ=255)
                                                                                   cAACCGCATCGCGCGCTCGGCCGAAACCCACCCGCTGACCTGGCGCCTGCGTGACGACGGCGAGCCGGTCTGGCTGGCCGAGTACGAATCGAAGAACGGCTACGCCGCAGCCCGCAAGGCGCTGGCGCAAATGTCCGCCGACGACATCg                                                               >  6:201865/1‑149 (MQ=255)
                                                                                    aaCCGCATCGCGCGCTCGGCCGAAACCCACCCGCTGACCTGGCGCCTGCGTGACGACGGCGAGCCGGTCTGGCTGGCCGAGTACGAATCGAAGAACGGCTACGCCGCAGCCCGCAAGGCGCTGGCGCAAATGTCCGCCGACGACATCGt                                                              >  3:347041/1‑149 (MQ=255)
                                                                                    aaCCGCATCGCGCGCTCGGCCGAAACCCACCCGCTGACCTGGCGCCTGCGTGACGACGGCGAGCCGGTCTGGCTGGCCGAGTACGAATCGAAGAACGGCTACGCCGCAGCCCGCAAGGCGCTGGCGCAAATGTCCGCCGACGACATCGt                                                              >  5:136858/1‑149 (MQ=255)
                                                                                         cATCGCGCGCTCGGCCGAAACCCACCCGCTGACCTGGCGCCTGCGTGACGACGGCGAGCCGGTCTGGCTGGCCGAGTACGAATCGAAGAACGGCTACGCCGCAGCCCGCAAGGCGCTGGCGCAAATGTCCGCCGACGACATCGTGCaga                                                         <  7:231897/149‑1 (MQ=255)
                                                                                         cATCGCGCGCTCGGCCGAAACCCACCCGCTGACCTGGCGCCTGCGTGACGACGGCGAGCCGGTCTGGCTGGCCGAGTACGAATCGAAGAACGGCTACGCCGCAGCCCGCAAGGCGCTGGCGCAAATGTCCGCCGACGACATCGTGCa                                                           >  7:321640/1‑147 (MQ=255)
                                                                                                 gcTCGGCCGAAACCCACCCGCTGACCTGGCGCCTGCGTGACGACGGCGAGCCGGTCTGGCTGGCCGAGTACGAATCGAAGAACGGCTACGCCGCAGCCCGCAAGGCGCTGGCGCAAATGTCCGCCGACGACATCGTGCAGAGCGTCAAg                                                 <  6:71769/149‑1 (MQ=255)
                                                                                                 gcTCGGCCGAAACCCACCCACTGACCTGGCGCCTGCGTGACGACGGCGAGCCGGTCTGGCTGGCCGAGTACGAATCGAAGAACGGCTACGCCGCAGCCCGCCAGGCGCTGGCGCAAATGTCCGCCGACGACATCGTGCAGAGCGTCAAg                                                 <  8:4924/149‑1 (MQ=255)
                                                                                                   tCGGCCGAAACCCACCCGCTGACCTGGCGCCTGCGTGACGACGGCGAGCCGGTCTGGCTGGCCGAGTACGAATCGAAGAACGGCTACGCCGCAGCCCGCAAGGCGCTGGCGCAAATGTCCGCCGACGACATCGTGCAGAGCGTCAAgg                                                >  2:350777/1‑148 (MQ=255)
                                                                                                    cGGCCGAAACCCCCCCGCTCCCCGGGCGCCGGCGGGACGCCGGCGCGCGGGTCGGGCGGGCCGATTACGAATCGAACACCGGCTCCGCCGCCGCCCGCAAGGCGCTGGCGCAAATGCCCGCCGACGACATCGTGCAGCGCGTCACGGAc                                              <  8:114432/149‑1 (MQ=255)
                                                                                                      gCCGAAACCCACCCGCTGACCTGGCGCCTGCGTGACGACGGCGAGCCGGTCTGGCTGGCCGAGTACGAATCGAAGAACGGCTACGCCGCAGCCCGCAAGGCGCTGGCGCAAATGTCCGCCGACGACATCGTGCAGAGCGTCAAGGACTc                                            >  5:45989/1‑149 (MQ=255)
                                                                                                      gCCGAAACCCACCCGCTGACCTGGCGCCTGCGTGACGACGGCGAGCCGGTCTGGCTGGCCGAGTACGAATCGAAGAACGGCTACGCCGCAGCCCGCAAGGCGCTGGCGCAAATGTCCGCCGACGACATCGTGCAGAGCGTCAAGGACTc                                            <  5:201865/149‑1 (MQ=255)
                                                                                                        cGAAACCCACCCGCTGACCTGGCGCCTGCGTGACGACGGCGAGCCGGTCTGGCTGGCCGAGTACGAATCGAAGAACGGCTACGCCGCAGCCCGCAAGGCGCTGGCGCAAATGTCCGCCGACGACATCGTGCAGAGCGTCAAGGACTcca                                          <  8:61749/149‑2 (MQ=255)
                                                                                                         gAAACCCACCCGCTGACCTGGCGCCTGCGTGACGACGGCGAGCCGGTCTGGCTGGCCGAGTACGAATCGAAGAACGGCTACGCCGCAGCCCGCAAGGCGCTGGCGCAAATGTCCGCCGACGACATCGTGCAGAGCGTCAAGGCCTCCgg                                         >  8:199301/1‑149 (MQ=255)
                                                                                                             cccacccGCTGACCTGGCGCCTGCGTGACGACGGCGAGCCGGTCTGGCTGGCCGAGTACGAATCGAAGAACGGCTACGCCGCAGCCCGCAAGGCGCTGGCGCAAATGTCCGCCGACGACATCGTGCAGAGCGTCAAGGACTCCGGCCTc                                     <  8:233403/149‑1 (MQ=255)
                                                                                                                acccGCTGACCTGGCGCCTGCGTGACGACGGCGAGCCGGTCTGGCTGGCCGAGTACGAATCGAAGAACGGCTACGCCGCAGCCCGCAAGGCGCTGGCGCAAATGTCCGCCGAGGACATCGTGCag                                                          <  8:106261/125‑1 (MQ=255)
                                                                                                                acccGCTGACCTGGCGCCTGCGTGACGACGGCGAGCCGGTCTGGCTGGCCGAGTACGAATCGAAGAACGGCTACGCCGCAGCCCGCAAGGCGCTGGCGCAAATGTCCGCCGACGACATCGTGCag                                                          >  7:106261/1‑125 (MQ=255)
                                                                                                                acccGCTGACCTGGCGCCTGCGTGACGACGGCGAGCCGGTCTGGCTGGCCGAGTACGAATCGAAGAACGGCTACGCCGCAGCCCGCAAGGCGCTGGCGCAAATGTCCGCCGACGACATCGTGCTGAGCGTCAAGGACTCCGGCCTCa                                    >  6:68162/1‑147 (MQ=255)
                                                                                                                          cTGGCGCCTGCGTGACGACGGCGAGCCGGTCTGGCTGGCCGAGTACGAATCGAAGAACGGCTACGCCGCAGCCCGCAAGGCGCTGGCGCAAATGTCCGCCGACGACATCGTGCAGAGCGTCAAGGACTCCGGCCTTTGAGGCCggagtc                        <  2:176261/149‑6 (MQ=255)
                                                                                                                               gcCTGCGTGACGACGGCGAGCCGGTCTGGCTGGCCGAGTACGAATCGAAGAACGGCTACGCCGCAGCCCGCAAGGCGCTGGCGCAAATGTCCGCCGACGACATCGTGCAGAGCGTCAAGGACTCCGGCCTCAAAGGCCGTGGCGGTGCa                   <  8:321640/149‑1 (MQ=255)
                                                                                                                                ccTGCGTGACGACGGCGAGCCGGTCTGGCTGGCCGAGTACGAATCGAAGAACGGCTACGCCGCAGCCCGCAAGGCGCTGGCGCAAATGTCCGCCGACGACATCGTGCAGAGCGTCAAGGACTCCGGCCTCAAAGGCCGTgg                          <  6:328076/141‑1 (MQ=255)
                                                                                                                                ccTGCGTGACGACGGCGAGCCGGTCTGGCTGGCCGAGTACGAATCGAAGAACGGCTACGCCGCAGCCCGCAAGGCGCTGGCGCAAATGTCCGCCGACGACATCGTGCAGAGCGTCAAGGACTCCGGCCTCAAAGGCCGTgg                          >  5:328076/1‑141 (MQ=255)
                                                                                                                                  tGCGTGACGACGGCGAGCCGGTCTGGCTGGCCGAGTACGAATCGAAGAACGGCTACGCCGCAGCCCGCAAGGCGCTGGCGCAAATGTCCGCCGACGACATCGTGCAGAGCGTCAAGGACTCCGGCCTCAAAGGc                               >  7:14513/1‑134 (MQ=255)
                                                                                                                                  tGCGTGACGACGGCGAGCCGGTCTGGCTGGCCGAGTACGAATCGAAGAACGGCTACGCCGCAGCCCGCAAGGCGCTGGCGCAAATGTCCGCCGACGACATCGTGCAGAGCGTCAAGGACTCCGGCCTCAAAGGc                               <  8:14513/134‑1 (MQ=255)
                                                                                                                                    cGTGACGACGGCGAGCCGGTCTGGCTGGCCGAGTACGAATCGAAGAACGGCTACGCCGCAGCCCGCAAGGCGCTGGCGCAAATGTCCGCCGACGACATCGTGCAGAGCGTCAAGGACTCCGGCCTCAAAGGCCGTGGCGGTGCAGGttt              <  1:350777/149‑1 (MQ=255)
                                                                                                                                       gacgacGGCGAGCCGGTCTGGCTGGCCGAGTACGAATCGAAGAACGGCTACGCCGCAGCCCGCAAGGCGCTGGCGCAAATGTCCGCCGACGACATCGTGCAGAGCGTCAAGGACTCCGGCCTCAAAGGCCGTGGCGGTGCAGGTTTccc           >  7:117091/1‑149 (MQ=255)
                                                                                                                                       gacgacGGCGAGCCGGTCTGGCTGGCCGAGTACGAATCGAAGAACGGCTACGCCGCAGCCCGCAAGGCGCTGGCGCAAATGTCCGCCGACGACATCGTGCAGAGCGTCAAGCACTCCGGCCTCAAAGGCCGTGGCGGTGCAGGTTTccc           >  5:6589/1‑149 (MQ=255)
                                                                                                                                        acgacgGCGAGCCGGTCTGGCTGGCCGAGTACGAATCGAAGAACGGCTACGCCGCAGCCCGCAAGGCGCTGGCGCAAATGTCCGCCGACGACATCGTGCAGAGCGTCAAGGACTCCg                                          >  3:8242/1‑117 (MQ=255)
                                                                                                                                        acgacgGCGAGCCGGTCTGGCTGGCCGAGTACGAATCGAAGAACGGCTACGCCGCAGCCCGCAAGGCGCTGGCGCAAATGTCCGCCGACGACATCGTGCAGAGCGTCAAGGACTCCg                                          <  4:8242/117‑1 (MQ=255)
                                                                                                                                           tcgGCCAGCCAGACCGGCTCGCCGAGTACGAATCGAAGAACGGCTACGCCGCAGCCCGCAAGGCGCTGGCGCAAATGTCCGCCGACGACATCGTGCAGAGCGTCAAg                                                 <  6:217405/106‑1 (MQ=255)
                                                                                                                                           tcgGCCAGCCAGACCGGCTCGCCGAGTACGAATCGAAGAACGGCTACGCCGCAGCCCGCAAGGCGCTGGCGCAAATGTCCGCCGACGACATCGTGCAGAGCGTCAAg                                                 >  5:217405/2‑107 (MQ=255)
                                                                                                                                              gCGAGCCGGTCTGGCTGGCCGAGTACGAATCGAAGAACGGCTACGCCGCAGCCCGCAAGGCGCTGGCGCAAATGTCCGCCGACGACATCGTGCAGAGCGTCAAGGACTCCGGCCTCAAAGGCCGTGGCGGTGCAGGTTTCCCCACTGGc    <  6:136858/149‑1 (MQ=255)
                                                                                                                                                gAGCCGGTCTGGCTGGCCGAGTACGAATCGAAGAACGGCTACGCCGCAGCCCGCAAGGCGCTGGCGCAAATGTCCGCCGACGACATCGTGCAGAGCGTCAAGGACTCCGGCCTCAAAGGCCGTGGCGGTGCAGGTTTCCCCACTGGCGt  >  5:301152/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
CGACACCTTCGGTGACGTGCAGCCGGCTGGCGTTTCCAAACTGCTGGAGGGTTACGTATGACCATTACTTCCTTCGGCCCGGCCAACCGCATCGCGCGCTCGGCCGAAACCCACCCGCTGACCTGGCGCCTGCGTGACGACGGCGAGCCGGTCTGGCTGGCCGAGTACGAATCGAAGAACGGCTACGCCGCAGCCCGCAAGGCGCTGGCGCAAATGTCCGCCGACGACATCGTGCAGAGCGTCAAGGACTCCGGCCTCAAAGGCCGTGGCGGTGCAGGTTTCCCCACTGGCGT  >  NC_002947/4659586‑4659878

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: ‑