Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,666,511 G→A 5.2% V93M (GTG→ATG)  nuoL → NADH‑quinone oxidoreductase subunit L

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,666,5110GA 5.2% 293.7 / 3.4 95V93M (GTG→ATG) nuoLNADH‑quinone oxidoreductase subunit L
Reads supporting (aligned to +/- strand):  ref base G (41/49);  new base A (3/2);  total (44/51)
Fisher's exact test for biased strand distribution p-value = 6.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.74e-01

ACGGCCGCCTACGTCATCTGGCAGTTCAACGTCGCCCCGCCTGAGGGCGGCGCGTACAGCCAGCTGCTGTGGCAGTGGATGTCGGTGGACGGCTTCGCGCCGAACTTCACCCTGTACCTGGACGGCCTGTCGGTCACCATGCTCGGCGTGGTGACCGGTGTCGGCTTCCTGATCCACCTGTTCGCATCCTGGTACATGCGTGGCGAAGCCGGCTACTCGCGCTTCTTCTCGTACACCAACCTGTTCATCGCCAGCATGCTGTTCCTGATCCTTGGCGATAACCTGCTGTTCATC  >  NC_002947/4666364‑4666657
                                                                                                                                                   |                                                                                                                                                   
aCGGCCGCCTACGTCATCTGGCAGTTCAACGTCGCCCCGCCTGAGGGCGGCGCGTACAGCCAGCTGCTGTGGCAGTGGATGTCGGTGGACGGCTTCGCGCCGAACTTCACCCTGTACCTGGACGGCCTGTCGGTCACCATGCTCGGCgt                                                                                                                                                    >  3:329551/1‑149 (MQ=255)
 cGGCCGCCTACGTCATCTGGCAGTTCAACGTCGCCCCGCCTGAGGGCGGCGCGTACAGCCAGCTGCTGTGGCAGTGGATGTCGGTGGACGGCTTCGCGCCGAACTTCACCCTGTACCTGGACGGCCTGTCGGTCACCATGCTCGGCgtg                                                                                                                                                   <  5:304524/149‑1 (MQ=255)
   gccgccTACGTCATCTGGCAGTTCAACGTCGCCCCGCCTGAGGGCGGCGCGTACAGCCAGCTGCTGTGGCAGTGGATGTCGGTGGACGGCTTCGCGCCGAACTTCACCCTGTACCTGGACGGCCTGTCGGTCACCATGCTCGGCgtggt                                                                                                                                                 <  6:292647/149‑1 (MQ=255)
   gccgccTACGTCATCTGGCAGTTCAACGTCGCCCCGCCTGAGGGCGGCGCGTACAGCCAGCTGCTGTGGCAGTGGATGTCGGTGGACGGCTTCGCGCCGAACTTCACCCTGTACCTGGACGGCCTGTCGGTCACCATGCTCGGCgtggt                                                                                                                                                 <  6:425/149‑1 (MQ=255)
     cgccTACGTCATCTGGCAGTTCAACGTCGCCCCGCCTGAGTGCGTCGCGTACAGCCAGCTGCTGTGGCAGTGGATGTCGGTGGTCGGCTTCGCGCCGAACTTCACCCTGTACCTGGACGGCCTGTCGGTCACCACGCCGAGCATGGTGa                                                                                                                                               <  7:74879/149‑1 (MQ=255)
      gccTACGTCATCTGGCAGTTCAACGTCGCCCCGCCTGAGGGCGGCGCGTACAGCCAGCTGCTGTGGCAGTGGATGTCGGTGGACGGCTTCGCGCCGAACTTCACCCTGTACCTGGACGGCCTGTCGGTCACCATGCTCGGCGTGGTGAc                                                                                                                                              >  1:271773/1‑149 (MQ=255)
      gccTACGTCATCTGGCAGTTCAACGTCGCCCCGCCTGAGGGCGGCGCGTACAGCCAGCTGCTGTGGCAGTGGATGTCGGTGGACGGCTTCGCGCCGAACTTCACCCTGTACCTGGACGGCCTGTCGGTCACCATGCTCGGCGTGGTGAc                                                                                                                                              <  6:293988/149‑1 (MQ=255)
        cTACGTCATCTGGCAGTTCAACGTCGCCCCGCCTGAGGGCGGCGCGTACAGCCAGCTGCTGTGGCAGTGGATGTCGGTGGACGGCTTCGCGCCGAACTTCACCCTGTACCTGGACGGCCTGTCGGTCACCATGCTCGGCGTGGTGACCg                                                                                                                                            <  4:75336/149‑1 (MQ=255)
        cTACGTCATCTGGCAGTTCAACGTCGCCCCGCCTGAGGGCGGCGCGTACAGCCAGCTGCTGTGGCAGTGGATGTCGGTGGACGGCTTCGCGCCGAACTTCACCCTGTACCTGGACGGCCTGTCGGTCACCATGCTCGGCGTGGTGACCg                                                                                                                                            <  4:369518/149‑1 (MQ=255)
        cTACGTCATCTGGCAGTTCAACGTCGCCCCGCCTGAGGGCGGCGCGTACAGCCAGCTGCTGTGGCAGTGGATGTCGGTGGACGGCTTCGCGCCGAACTTCACCCTGTACCTGGACGGCCTGTCGGTCACCATGCTCGGCGTGGTGACCg                                                                                                                                            <  7:138344/149‑1 (MQ=255)
           cGTCATCTGGCAGTTCAACGTCGCCGCGCCTGAGGGCGGCGCGTACAGCCAGCTGCTGTGGCAGTGGATGTCGGTGGACGGCTTCGCGCCGAACTTCACCCTGTACCTGGACGGCCTGTCGGTCACCATGCTCGGCGTGGTGACCGgtg                                                                                                                                         >  3:44948/1‑149 (MQ=255)
           cGTCATCTGGCAGTTCAACGTCGCCGCGCCTGAGGGCGGCGCGTACAGCCAGCTGCTGTGGCAGTGGATGTCGGTGGACGGCTTCGCGCCGAACTTCACCCTGTACCTGGACGGCCTGCCGGTCAGCACGCCGAGCATGGTGACCgaca                                                                                                                                         <  6:215315/149‑4 (MQ=255)
           cGTCATCTGGCAGTTCAACGTCGCCCCGCCTGAGGGCGGCGCGTACAGCCAGCTGCTGTGGCAGTGGATGTCGGTGGACGGCTTCGCGCCGAACTTCACCCTGTACCTGGACGGCCTGTCGGTCACCATGCTCGGCGTGGTGACCGgtg                                                                                                                                         <  8:83656/149‑1 (MQ=255)
              cATCTGGCAGTTCAACGTCGCCCCGCCTGAGGGCGGCGCGTACAGCCAGCTGCTGTGGCAGTGGATGTCGGTGGACGGCTTCGCGCCGAACTTCACCCTGTACCTGGACGGCCTGTCGGTCACCATGCTCGGCGTGGTGACCGGTGTCg                                                                                                                                      >  8:73604/1‑149 (MQ=255)
                tCTGGCAGTTCAACGTCGCCCCGCCTGAGGGCGGCGCGTACAGCCAGCTGCTGTGGCAGTGGATGTCGGTGGACGGCTTCGCGCCGAACTTCACCCTGTACCTGGACGGCCTGTCGGTCACCATGCTCGGCGTGGTGACCGGTGTCGGc                                                                                                                                    >  1:369569/1‑149 (MQ=255)
                  tGGCAGTTCAACGTCGCCCCGCCTGAGGGCGGCGCGTACAGCCAGCTGCTGTGGCAGTGGATGTCGGTGGACGGCTTCGCGCCGAACTTCACCCTGTACCTGGACGGCCTGTCGGTCACCATGCTCGGCGTGGTGACCGGTGTCGGCtt                                                                                                                                  >  7:160443/1‑149 (MQ=255)
                   ggCAGTTCAACGTCGCCCCGCCTGAGGGCGGCGCGTACAGCCAGCTGCTGTGGCAGTGGATGTCGGTGGACGGCTTCGCGCCGAACTTCACCCTGTACCTGGACGGCCTGTCGGTCACCATGCTCGGCGTGGTGACCGGTGTCGGCTTc                                                                                                                                 >  7:75383/1‑149 (MQ=255)
                           aaCGTCGCCCCGCCTGAGGGCGGCGCGTACAGCCAGCTGCTGTGGCAGTGGATGTCGGTGGACGGCTTCGCGCCGAACTTCACCCTGTACCTGGACGGCCTGTCGGTCACCATGCTCGGCGTGGTGACCGGTGTc                                                                                                                                       <  3:90008/135‑1 (MQ=255)
                           aaCGTCGCCCCGCCTGAGGGCGGCGCGTACAGCCAGCTGCTGTGGCAGTGGATGTCGGTGGACGGCTTCGCGCCGAACTTCACCCTGTACCTGGACGGCCTGTCGGTCACCATGCTCGGCGTGGTGACCGGTGTc                                                                                                                                       >  4:90008/1‑135 (MQ=255)
                           aaCGTCGCCCCGCCTGAGGGCGGCGCGTACAGCCAGCTGCTGTGGCAGTGGATGTCGGTGGACGGCTTCGCGCCGAACTTCACCCTGTACCTGGACGGCCTGTCGGTCACCATGCTCGGCGTGGTGACCGGTGTCGGCTTCCTGATcc                                                                                                                         <  8:75383/148‑1 (MQ=255)
                              gTCGCCCCGCCTGAGGGCGGCGCGTACAGCCAGCTGCTGTGGCAGTGGATGTCGGTGGACGGCTTCGCGCCGAACTTCACCCTGTACCTGGACGGCCTGTCGGTCACCATGCTCGGCGTGGTGACCGGTGTCGGCTTCCTGATCCACCt                                                                                                                     <  4:361125/149‑1 (MQ=255)
                              gTCGCCCCGCCTGAGGGCGGCGCGTACAGCCAGCTGCTGTGGCAGTGGATGTCGGTGGACGGCTTCGCGCCGAACTTCACCCTGTACCTGGACGGCCTGTCGGTCACCATGCTCGGCGTGGTGACCGGTGTCGGCTTCCTGATCCACCt                                                                                                                     <  7:73604/149‑1 (MQ=255)
                               tCGCCCCGCCTGAGGGCGGCGCGTACAGCCAGCTGCTGTGGCAGTGGATGTCGGTGGACGGCTTCGCGCCGAACTTCACCCTGTACCTGGACGGCCTGTCGGTCACCATGCTCGGCGTGGTGACCGGTGTCGGCTTCCAGATCCCCCTg                                                                                                                    >  5:105599/1‑149 (MQ=255)
                                 gCCCCGCCTGAGGGCGGCGCGTACAGCCAGCTGCTGTGGCAGTGGATGTCGGTGGACGGCTTCGCGCCGAACTTCACCCTGTACCTGGACGGCCTGTCGGTCACCATGCTCGGCGTGGTGACCGGAGTCGGCTTCCTGATCCACCTGtt                                                                                                                  >  8:39935/1‑149 (MQ=255)
                                  ccccGCCTGAGGGCGGCGCGTACAGCCAGCTGCTGTGGCAGTGGATGTCGGTGGACGGCTTCGCGCCGAACTTCACCCTGTACCTGGACGGCCTGTCGGTCACCATGCTCGGCGTGGTGACCGGTGTCGGCTTCCTGATCCACCTGTTc                                                                                                                 <  5:181257/149‑1 (MQ=255)
                                  ccccGCCTGAGGGCGGCGCGTACAGCCAGCTGCTGTGGCAGTGGATGTCGGTGGACGGCTTCGCGCCGAACTTCACCCTGTACCTGGACGGCCTGTCGGTCACCATGCTCGGCGTGGTGACCGGTGTCGGCTTCCTGATCCACCTGTTc                                                                                                                 <  8:319046/149‑1 (MQ=255)
                                    ccGCCTGAGGGCGGCGCGTACAGCCAGCTGCTGTGGCAGTGGATGTCGGTGGACGGCTTCGCGCCGAACTTCACCCTGTACCTGGACGGCCTGTCGGTCACCATGCTCGGCGTGGTGACCGGTGTCGGCTTCCTGATCCACCTGTTCg                                                                                                                >  2:3190/1‑148 (MQ=255)
                                            gggCGGCGCGTACAGCCAGCTGCTGTGGCAGTGGATGTCGGTGGACGGCTTCGCGCCGAACTTCACCCTGTACCTGGACGGCCTGTCGGTCACCATGCTCTGCGTGGTGAGCGGTGTCGGCTTCCTGATCCACCTGTTCGCATCCTGGt                                                                                                       <  7:39935/149‑1 (MQ=255)
                                            gggCGGCGCGTACAGCCAGCTGCTGTGGCAGTGGATGTCGGTGGACGGCTTCGCGCCGAACTTCACCCTGTACCTGGACGGCCTGTCGGTCACCATGCTCGGCGTGGTGACCGGTGTCGGCTTCCTGATCCACCTGTTCGCATCCTGGt                                                                                                       <  8:160443/149‑1 (MQ=255)
                                             ggcggcGCGTACAGCCAGCTGCTGTGGCAGTGGATGTCGGTGGACGGCTTCGCGCCGAACTTCACCCTGTACCTGGACGGCCTGTCGGTCACCATGCTCGGCGTGGTGACCGGTGTCGGCTTCCTGATCCACCTGTTCGCATCCTGGTa                                                                                                      <  2:310247/149‑1 (MQ=255)
                                              gcggcgCGTACAGCCAGCTGCTGTGGCAGTGGATGTCGGTGGACGGCTTCGCGCCGAACTTCACCCTGTACCTGGACGGCCTGTCGGTCACCATGCTCGGCGTGGTGACCGGTGTCGGCTTCCTGATCCACCTGTTCGCATCCTGGTAc                                                                                                     <  6:102673/149‑1 (MQ=255)
                                              gcggcgCGTACAGCCAGCTGCTGTGGCAGTGGATGTCGGTGGACGGCTTCGCGCCGAACTTCACCCTGTACCTGGACGGCCTGTCGGTCACCATGCTCGGCGTGGTGACCGGTGTCGGCTTCCTGATCCACCTGTTCGCATCCTGGTAc                                                                                                     <  5:334483/149‑1 (MQ=255)
                                                ggcgCGTACAGCCAGCTGCTGTGGCAGTGGATGTCGGTGGACGGCTTCGCGCCGAACTTCACCCTGTACCTGGACGGCCTGTCGGTCACCATGCTCGGCGTGGTGACCGGTGTCgg                                                                                                                                     <  2:283077/116‑1 (MQ=255)
                                                ggcgCGTACAGCCAGCTGCTGTGGCAGTGGATGTCGGTGGACGGCTTCGCGCCGAACTTCACCCTGTACCTGGACGGCCTGTCGGTCACCATGCTCGGCGTGGTGACCGGTGTCgg                                                                                                                                     >  1:283077/1‑116 (MQ=255)
                                                    cgTACAGCCAGCTGCTGTGGCAGTGGATGTCGGTGGACGGCTTCGCGCCGAACTTCACCCTGTACCTGGACGGCCTGTCGGTCACCATGCTCGGCGTGGTGACCGGTGTCGGCTTCCTGATCCACCTGTTCGCATCCTGGTACATGCGt                                                                                               >  1:33647/1‑149 (MQ=255)
                                                    cgTACAGCCAGCTGCTGTGGCAGTGGATGTCGGTGGACGGCTTCGCGCCGAACTTCACCCTGTACCTGGACGGCCTGTCGGTCACCATGCTCGGCGTGGTGACCGGTGTCGGCTTCCTGATCCACCTGTTCGCATCCTGGTACATGCGt                                                                                               <  1:250620/149‑1 (MQ=255)
                                                      tACAGCCAGCTGCTGTGGCAGTGGATGTCGGTGGACGGCTTCGCGCCGAACTTCACCCTGTACCTGGACGGCCTGTCGGTCACCATGCTCGGCGTGGTGACCGGTGTCGGCTTCCTGATCCACCTGTTCGCATCCTGGTACATGCGTgg                                                                                             <  2:240408/149‑1 (MQ=255)
                                                      tACAGCCAGCTGCTGTGGCAGTGGATGTCGGTGGACGGCTTCGCGCCGAACTTCACCCTGTACCTGGACGGCCTGTCGGTCACCATGCTCGGCGTGGTGACCGGTGTCGGCTTCCTGATCCACCTGTTCGCATCCTGGTACATGCGTgg                                                                                             <  7:219068/149‑1 (MQ=255)
                                                          gccagcTGCTGTGGCAGTGGATGTCGGTGGACGGCTTCGCGCCGAACTTCACCCTGTACCTGGACGGCCTGTCGGTCACCATGCTCGGCGTGGTGACCGGTGTCGGCTTCCTGATCCACCTGTTCGCATCCTGGTACATGCGTGGCGaa                                                                                         <  2:341718/149‑1 (MQ=255)
                                                           ccagcTGCTGTGGCAGTGGATGTCGGTGGACGGCTTCGCGCCGAACTTCACCCTGTACCTGGACGGCCTGTCGGTCACCATGCTCGGCGTGGTGACCGGTGTCGGCTTCCTGATCCACCTGTTCGCATCCTGGTACATGCGTGGCGAAg                                                                                        >  3:47488/1‑149 (MQ=255)
                                                           ccagcTGCTGTGGCAGTGGATGTCGGTGGACGGCTTCGCGCCGAACTTCACCCTGTACCTGGACGGCCTGTCGGTCACCATGCTCGGCGTGGTGACCGGTGTCGGCTTCCTGATCCACCTGTTCGCATCCTGGTACATGCGTGGCGAAg                                                                                        >  5:102857/1‑149 (MQ=255)
                                                            cagcTGCTGTGGCAGTGGATGTCGGTGGACGGCTTCGCGCCGAACTTCACCCTGTACCTGGACGGCCTGTCGGTCACCATGCTCGGCGTGGTGACCGGTGTCGGCTTCCTGATCCACCTGTTCGCATCCTGGTACATGCGTGGCGAAGc                                                                                       <  3:238142/149‑1 (MQ=255)
                                                               ctgctgTGGCAGTGGATGTCGGTGGACGGCTTCGCGCCGAACTTCACCCTGTACCTGGACGGCCTGTCGGTCACCATGCTCGGCGTGGTGACCGGTGTCGGCTTCCTGATCCACCTGTTCGCATCCTGGTACATGCGTGGCGAAGCCgg                                                                                    <  1:182453/149‑1 (MQ=255)
                                                               ctgctgTGGCAGTGGATGTCGGTGGACGGCTTCGCGCCGAACTTCACCCTGTACCTGGACGGCCTGTCGGTCACCATGCTCGGCGTGGTGACCGGTGTCGGCTTCCTGATCCACCTGTTCGCATCCTGGTACATGCGTGGCGAAGCCgg                                                                                    <  4:47488/149‑1 (MQ=255)
                                                               ctgctgTGGCAGTGGATGTCGGTGGACGGCTTCGCGCCGAACTTCACCCTGTACCTGGACGGCCTGTCGGTCACCATGCTCGGCGTGGTGACCGGTGTCGGCTTCCTGATCCACCTGTTCGCATCCTGGTACATGCGTGGCGAAGCCgg                                                                                    <  4:44948/149‑1 (MQ=255)
                                                                tgctgTGGCAGTGGATGTCGGTGGACGGCTTCGCGCCGAACTTCACCCTGTACCTGGACGGCCTGTCGGTCACCATGCTCGGCGTGGTGACCGGTGTCGGCTTCCTGATCCACCTGTTCGCATCCTGGTACATGCGTgg                                                                                             >  4:56140/1‑139 (MQ=255)
                                                                tgctgTGGCAGTGGATGTCGGTGGACGGCTTCGCGCCGAACTTCACCCTGTACCTGGACGGCCTGTCGGTCACCATGCTCGGCGTGGTGACCGGTGTCGGCTTCCTGATCCACCTGTTCGCATCCTGGTACATGCGTgg                                                                                             <  3:56140/139‑1 (MQ=255)
                                                                tgctgTGGCAGTGGATGTCGGTGGACGGCTTCGCGCCGAACTTCACCCTGTACCTGGACGGCCTGTCGGTCACCATGCTCGGCGTGGTGACCGGTGTCGGCTTCCTGATCCACCTGTTCGCATCCTGGTACATGCGTGGCGAAGCCGGc                                                                                   >  7:308731/1‑149 (MQ=255)
                                                                       gCAGTGGATGTCGGTGGACGGCTTCGCGCCGAGCTTCACCCTGTACCTGGACGGCCTGTCGGTCACCATGCTCGGCGTGGTGACCGGTGTCGGCTTCCTGATCCACCTGTTCGCATCCTGGTACATGCGTGGCGAAGCCGGCTACTcgc                                                                            >  7:24163/1‑149 (MQ=255)
                                                                        cAGTGGATGTCGGTGGACGGCTTCGCGCCGAACTTCACCCTGTACCTGGACGGCCTGTCGGTCACCATGCTCGGCGTGGTGACCGGTGTCGGCTTCCTGATCCACCTGTTCGCATCCTGGTACATGCGTGGCGAAGCCGGCTACTcgcg                                                                           >  8:78555/1‑149 (MQ=255)
                                                                          gTGGATGTCGGTGGACGGCTTCGCGCCGAACTTCACCCTGTACCTGGACGGCCTGTCGGTCACCATGCTCGGCGTGGTGACCGGTGTCGGCTTCCTGATCCAc                                                                                                                       <  4:114385/103‑1 (MQ=255)
                                                                          gTGGATGTCGGTGGACGGCTTCGCGCCGAACTTCACCCTGTACCTGGACGGCCTGTCGGTCACCATGCTCGGCGTGGTGACCGGTGTAGGCTTCCTGATCCAc                                                                                                                       >  3:114385/1‑103 (MQ=255)
                                                                               tGTCGGTGGACGGCTTCGCGCCGAACTTCACCCTGTACCTGGACGGCCTGTCGGTCACCATGCTCGGCGTGGTGACCGGTGTCGGCTTCCTGATCCACCTGTTCGCATCCTGGTACATGCGTGGCGAAGCCGGCTACTCGCGCttcttc                                                                    <  5:383193/149‑1 (MQ=255)
                                                                                gTCGGTGGACGGCTTCGCGCCGAACTTCACCCTGTACCTGGACGGCCTGTCGGTCACCATGCTCGGCGTGGAATCCGGTGTCGGCTTCCAGGATCCACCTGTTCGCATCCTGGTCCATGCGTGGCGAAGCCGGCTACTCGCGCttcttc                                                                    >  1:316012/1‑149 (MQ=255)
                                                                                       gACGGCTTCGCGCCGAACTTCACCCTGTACCTGGACGGCCTGTCGGTCACCATGCTCGGCGTGGTGACCGGTGTCGGCTTCCTGATCCACCTGTTCGCATCCTGGTACATGCGTGGCGAAGCCGGCTACTCGCGCTTCTTCTCGTacac                                                            <  4:329551/149‑1 (MQ=255)
                                                                                          ggCTTCGCGCCGAACTTCACCCTGTACCTGGACGGCCTGTCGGTCACTATGCTCGGCGTGGTGACCGGTGTCGGCTTCCTGATCCACCTGTTCGCATCCTGGTACATGCGTGGCGAAGCCGGCTACTCGCGCTTCTTCTCGTACACCaa                                                         >  8:35146/1‑149 (MQ=255)
                                                                                            cTTCGCGCCGAACTTCACCCTGTACCTGGACGGCCTGTCGGTCACCATGCTCGGCGTGGTGACCGGTGTCGGCTTCCTGATCCACCTGTTCGCATCCTGGTACATGCGTGGCGAAGCCGGCTACTCGCGCTTCTTCTCGTACACCAAcc                                                       >  4:188531/1‑149 (MQ=255)
                                                                                                  gcCGAACTTCACCCTGTACCTGGACGGCCTGTCGGTCACCATGCTCGGCGTGGTGACCGGTGTCGGCTTCCTGATCCACCTGTTCGCATCCTGGTACATGCGTGGCGAAGCCGGCTACTCGCGCTTCTTCTCGTACACCAACCTGTTc                                                  >  6:62637/1‑148 (MQ=255)
                                                                                                  gcCGAACTTCACCCTGTACCTGGACGGCCTGTCGGTCACCATGCTCGGCGTGGTGACCGGTGTCGGCTTCCTGATCCACCTGTTCGCATCCTGGTACATGCGTGGCGAAGCCGGCTACTCGCGCTTCTTCTCGTACACCAACCTGTTCa                                                 >  8:6867/1‑149 (MQ=255)
                                                                                                    cGAACTTCACCCTGTACCTGGACGGCCTGTCGGTCACCATGCTCGGCGTGGTGACCGGTGTCGGCTTCCTGATCCACCTGTTCGCATCCTGGTACATGCGTGGCGAAGCCGGCTACTCGCGCTTCTTCTCGTACACCAACCTGTTCATc                                               <  1:3190/149‑1 (MQ=255)
                                                                                                      aaCTTCACCCTGTACCTGGACGGCCTGTCGGTCACCATGCTCGGCGTGGTGACCGGTGTCGGCTTCCTGATCCACCTGTTCGCATCCTGGTACATGCGTGGCGAAGCCGGCTACTCGCGCTTCTTCTCGTACACCAACCTGTTCATCGc                                             <  7:78555/149‑1 (MQ=255)
                                                                                                       aCTTCACCCTGTACCTGGACGGCCTGTCGGTCACCATGCTCGGCGTGGTGACCGGTGTCGGCTTCCTGATCCACCTGTTCGCATCCTGGTACATGCGTGGCGAAGCCGGCTACTCGCGCTTCTTCTCGTACACCAACCTGTTCATCGcc                                            <  1:154254/149‑1 (MQ=255)
                                                                                                         ttCACCCTGTACCTGGACGGCCTGTCGGTCACCATGCTCGGCGTGGTGACCGGTGTCGGCTTCCTGATCCACCTGTTCGCATCCTGGTACATGCGTGGCGAAGCCGGCTACTCGCGCTTCTTCTCGTACACCAACCTGTTCATCGCCAg                                          >  3:302856/1‑149 (MQ=255)
                                                                                                           cACCCTGTACCTGGACGGCCTGTCGGTCACCATGCTCGGCGTGGTGACCGGTGTCGGCTTCCTGATCCACCTGTTCGCATCCTGGTACATGCGTg                                                                                              >  8:95542/1‑95 (MQ=255)
                                                                                                           cACCCTGTACCTGGACGGCCTGTCGGTCACCATGCTCGGCGTGGTGACCGGTGTCGGCTTCCTGATCCACCTGTTCGCATCCTGGTACATGCGTg                                                                                              <  7:95542/95‑1 (MQ=255)
                                                                                                           cACCCTGTACCTGGACGGCCTGTCGGTCACCATGCTCGGCGTGGTGACCGGTGTCGGCTTCCTGATCCACCTGTTCGCATCCTGGTACATGCGTGGCGAAGCCGGCTACTCGCGCTTCTTCTCGTACACCAACCTGTTCATCGCCAGc                                         <  5:62637/148‑1 (MQ=255)
                                                                                                           cACCCTGTACCTGGACGGCCTGTCGGTCACCATGCTCGGCGTGGTGACCGGTGTCGGCTTCCTGATCCACCTGTTCGCATCCTGGTACATGCGTGGCGAAGCCGGCTACTCGCGCTTCTTCTCGTACACCAACCTGTTCATCGCCAGCa                                        >  5:231732/1‑149 (MQ=255)
                                                                                                            aCCCTGTACCTGGACGGCCTGTCGGTCACCATGCTCGGCGTGGTGACCGGTGTCGGCTTCCTGATCCACCTGTTCGCATCCTGGTACATGCGTGGCGAAGCCGGCTACTCGCGCTTCTTCTCGTACACCAACCTGTTCATCGCCAGCAt                                       >  7:82334/1‑149 (MQ=255)
                                                                                                                  tACCTGGACGGCCTGTCGGTCACCATGCTCGGCGTGGTGACCGGTGTCGGCTTCCTGATCCACCTGTTCGCATCCTGGTACATGCGTGGCGAAGCCGGCTACTCGCGCTTCTTCTCGTACACCAACCTGTTCATCGCCAGCATGCTg                                   <  8:82334/147‑1 (MQ=255)
                                                                                                                   aCCTGGACGGCCTGTCGGTCACCATGATCGGCGTGGTGACCGGTGTCGGCTTCCTGATCCACCTGTTCGCATCATGGTACAAGCGTGGCGAAGCCGGCTACTCGAGCATCTTCTCGTACACCAACCTGATCAACGCCAGCatgatgat                                 >  4:36211/1‑146 (MQ=255)
                                                                                                                    ccTGGTCGGCCTGTCGTTCACCATGCTCGGCGTTGTGAGCGGTGTCGGCTTCCTGATCCACCTGTTCGCATCCTGGTACATGCGTGGCGAAGCCGGCTACTCGCGCTTCTTCTCGTACACCAACCTGTTCATCGCCAGCATGCTGTTc                                <  3:36211/148‑1 (MQ=255)
                                                                                                                    ccTGGACGGCCTGTCGGTCACCATGCTCGGCGTGGTGACCGGTGTCGGCTTCCTGATCCACCTGTTCGCATCCTGGTACATGCGTGGCGAAGCCGGCTACTCGCGCTTCTTCTCGTACACCAACCTGTTCATCGCCAGCATGCTGTTcc                               <  2:33647/149‑1 (MQ=255)
                                                                                                                    ccTGGACGGCCTGTCGGTCACCATGCTCGGCGTGGTGACCGGTGTCGGCTTCCTGATCCACCTGTTCGCATCCTGGTACATGCGTGGCGAAGCCGGCTACTCGCGCTTCTTCTCGTACACCAACCTGTTCATCGCCAGCATGCTGTTcc                               >  4:264562/1‑149 (MQ=255)
                                                                                                                     cTGGACGGCCTGTCGGTCACCATGCTCGGCGTGGTGACCGGTGTCGGCTTCCTGATCCACCTGTTCGCATCCTGGTACATGCGTGGCGAAGCCGGCTACTCGCGCTTCTTCTCGTACACCAACCTGTTCATCGCCAGCATGCTGTTCCt                              >  4:134189/1‑149 (MQ=255)
                                                                                                                     cTGGACGGCCTGTCGGTCACCATGCTCGGCGTGGTGACCGGTGTCGGCTTCCTGATCCACCTGTTCGCATCCTGGTACATGCGTGGCGAAGCCGGCTACTCGCGCTTCTTCTCGTACACCAACCTGTTCATCGCCAGCATGCTGTTCCt                              >  5:41259/1‑149 (MQ=255)
                                                                                                                     cTGGACGGCCTGTCGGTCACCATGCTCGGCGTGGTGACCGGTGTCGGCTTCCTGATCCACCTGTTCGCATCCTGGTACATGCGTGGCGAAGCCGGCTACTCGCGCTTCTTCTCGTACACCAACCTGTTCATCGCCAGCATGCTGTTCCt                              <  4:302856/149‑1 (MQ=255)
                                                                                                                        gACGGCCTGTCGGTCACCATGCTCGGCGTGGTGACCGGTGTCGGCTTCCTGATCCACCTGTTCGCATCCTGGTACATGCGTGGCGAAGCCGGCTACTCGCGCTTCTTCTCGTACACCAACCTGTTCATCGCCAGCATGCTGTTCCTGa                            <  6:231732/148‑1 (MQ=255)
                                                                                                                         aCGGCCTGTCGGTCACCATGCTCGGCGTGGTGACCGGTGTCGGCTTCCTGATCCACCTGTTCGCATCCTGGTACATGCGTGGCGAAGCCGGCTACTCGCGcttc                                                                       >  5:91008/1‑104 (MQ=255)
                                                                                                                         aCGGCCTGTCGGTCACCATGCTCGGCGTGGTGACCGGTGTCGGCTTCCTGATCCACCTGTTCGCATCCTGGTACATGCGTGGCGAAGCCGGCTACTCGCGcttc                                                                       <  6:91008/104‑1 (MQ=255)
                                                                                                                             gacaccGGTCACCACGCCGAGCATGGTGACCGGTGTCGGCTTCCTGATCCACCTGTTCGCATCCTGGTACATGCGTGGCGAAGCCGGCTACTCGCGCTTCTTCTCGTACACCAACCTGTTCATCGCCAGCATGCTGTTCCTGATCCTTg                      >  6:47313/6‑149 (MQ=255)
                                                                                                                              cTGTCGGTCACCATGCTCGGCGTGGTGACCGGTGTCGGCTTCCTGATCCACCTGTTCGCATCCTGGTACATGCGTGGCGAAGCCGGCTACTCGCGCTTCTTCTCGTACACCAACCTGTTCATCGCCAGCATGCTGTTCCTGATCCTTgg                     <  3:134189/149‑1 (MQ=255)
                                                                                                                              cTGTCGGTCACCATGCTCGGCGTGGTGACCGGTGTCGGCTTCCTGATCCACCTGTTCGCATCCTGGTACATGCGTGGCGAAGCCGGCTACTCGCGCTTCTTCTCGTACACCAACCTGTTCATCGCCAGCATGCTGTTCCTGATCCTTgg                     >  6:68965/1‑149 (MQ=255)
                                                                                                                               tGTCGGTCACCATGCTCGGCGTGGTGACCGGTGTCGGCTTCCTGATCCACCTGTTCGCATCCTGGTACATGCGTGGCGAAGCCGGCTACTCGCGCTTCTTCTCGTACACCAACCTGTTCATCGCCAGCATGCTGTTCCTGATCCTTGGc                    >  5:60312/1‑149 (MQ=255)
                                                                                                                                 tCGGTCACCATGCTCGGCGTGGTGACCGGTGTCGGCTTCCTGATCCACCTGTTCGCATCCTGGTACATGCGTGGCGAAGCCGGCTACTCGCGCTTCTTCTCGTACACCAACCTGTTCATCGCCAGCATGCTGTTCCTGATCCTTGGCGa                  >  2:202147/1‑149 (MQ=255)
                                                                                                                                  cGGTCACCACGCCGAGCATGGTGACCGGTGTCGGCTTCCTGATCCACCTGTTCGCATCCTGGTACATGCGTGGCGAAGCCGGCTACTCGCGCTTCTTCTCGTACACCAACCTGTTCATCGCCAGCATGCTGTTCCTGATCCTTGGCGa                  >  8:229215/1‑148 (MQ=255)
                                                                                                                                   ggTCACCACGCCGAGCATGGTGACCGGTGTCGGCTTCCAGATCCACCTCTTCGCATCCTGGTACATGCGTGGCGAAGCCGGCTACTCGCGCTTCTTCTCGTACACCAACCTGTTCATCGCCAGCATGCTGTTCCTGATACTTGGCGAt                 >  5:7090/1‑148 (MQ=255)
                                                                                                                                    gTCACCATGCTCGGCGTGGTGACCGGTGTCGGCTTCCTGATCCACCTGTTCGCATCCTGGTACATGCGTGGCGAAGCCGGCTACTCGCGCTTCTTCGCGTACACCACCCTGTTCATCGCCAGCATGCTGTTCCTGATCCTTGGCGATaa               >  7:176552/1‑149 (MQ=255)
                                                                                                                                     tCACCATGCTCGGCGTGGTGACCGGTGTCGGCTTCCTGATCCACCTGTTCGCATCCTGGTACATGCGTGGCGAAGCCGGCTACTCGCGCTTCTTCTCGTACACCAACCTGTTCATCGCCAGCATGCTGTTCCTGATCCTTGGCGATAAc              <  2:316012/149‑1 (MQ=255)
                                                                                                                                     tCACCATGCTCGGCGTGGTGACCGGTGTCGGCTTCCTGATCCACCTGTTCGCATCCTGGTACATGCGTGGCGAAGCCGGCTACTCGCGCTTCTTCTCGTACACCAACCTGTTCATCGCCAGCATGCTGTTCCTGATCCTTGGCGATAAc              <  2:271773/149‑1 (MQ=255)
                                                                                                                                        ccATGCTCGGCGTGGTGACCGGTGTCGGCTTCCTGATCCACCTGTTCGCATCCTGGTACATGCGTGGCGAAGCCGGCTACTCGCGCTTCTTCTCGTACACCAACCTGTTCATCGCCAGCATg                                      <  4:52802/122‑1 (MQ=255)
                                                                                                                                        ccATGCTCGGCGTGGTGACCGGTGTCGGCTTCCTGATCCACCTGTTCGCATCCTGGTACATGCGTGGCGAAGCCGGCTACTCGCGCTTCTTCTCGTACACCAACCTGTTCATCGCCAGCATg                                      >  3:52802/1‑122 (MQ=255)
                                                                                                                                         cATGCTCGTCGTGGTGACCGGTGTCGGCTTCATGATCCACCTGTTCGCATCCTTGTACATGCGTGGCGAAGCCGGCTACTAGCGCTTCTTCTCGTACACCAACCTGTTCATCGCCAGCATGCTGTTCCTGATCCttt                      <  7:340314/137‑2 (MQ=255)
                                                                                                                                         cATGCTCGGCGTGGTGACCGGTGTCGGCTTCCTGATCCACCTGTTCGCATCCTGGTACATGCGTGGCGAAGCCGGCTACTCGCGCTTCTTCTCGTACACCAACCTGTTCATCGCCAGCATGCTGTTCCTGATCCTTGGCGATAACctgc          >  3:167831/1‑149 (MQ=255)
                                                                                                                                         cATGCTCGGCGTGGTGAACGGTGTCGGCTTACTGATCCACCTGTTCGCATCATGGTACATGCGTGGCGAAGCCGGCTAATCGCGCTTCTTCTCGTACACCAACCTCTTAATCGCCAGCAAGCTGTTCCTGATCCTTGGAGAtcgcaaga          >  8:340314/1‑142 (MQ=255)
                                                                                                                                          aTGCTCGGCGTGGTGACCGGTGTCGGCTTCCTGATCCACCTGTTCGCATCCTGGTACATGCGTGGCGAAGCCGGCTACTCGCGCTTCTTCTCGTACACCAACCTGTTCATCGCCAGCATGCTGTTCCTGATCCTTGGCGATAACctgct         <  6:60312/149‑1 (MQ=255)
                                                                                                                                            gCTCGGCGTGGTGACCGGTGTCGGCTTCCTGATACACCTCTTCGCATACTGGTAAATGCGTGGCGAAGCAGGCTACTAGCGCTTCTTATCGTACACAAAACTGTTCATCGCCAGCATACTGTAACTGATCCTAGGAGATAACCTTCTGt       >  2:223036/1‑149 (MQ=255)
                                                                                                                                                 gCGTGGTGACCGGTGTCGGCTTCCTGATCCACCTGTTCGCATCCTGGTACATGCGTGGCGAAGCCGGCTACTCGCGCTTCTTCTCGTACACCAACCTGTTCATCGCCAGCATGCTGTTCCTGATCCTTGGCGATAACCTGCTGTTCATc  <  8:308731/149‑1 (MQ=255)
                                                                                                                                                 gCGTGGTGACCGGTGTCGGCTTCCTGATCCACCTGTTCGCATCCTGGTACATGCGTGGCGAAGCCGGCTACTCGCGCTTCTTCTCGTACACCAACCTGTTCATCGCCAGCATGCTGTTCCTGATCCTTGGCGATAACCTGCTGTTCATc  <  2:369569/149‑1 (MQ=255)
                                                                                                                                                 gCGTGGTGACCGGTGTCGGCTTCCTGATCCACCTGTTCGCATCCTGGTACATGCGTGGCGAAGCCGGCTACTCGCGCTTCTTCTCGTACACCAACCTGTTCATCGCCAGCATGCTGTTCCTGATCCTTGGCGATAACCTGCTGTTCATc  <  6:102857/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ACGGCCGCCTACGTCATCTGGCAGTTCAACGTCGCCCCGCCTGAGGGCGGCGCGTACAGCCAGCTGCTGTGGCAGTGGATGTCGGTGGACGGCTTCGCGCCGAACTTCACCCTGTACCTGGACGGCCTGTCGGTCACCATGCTCGGCGTGGTGACCGGTGTCGGCTTCCTGATCCACCTGTTCGCATCCTGGTACATGCGTGGCGAAGCCGGCTACTCGCGCTTCTTCTCGTACACCAACCTGTTCATCGCCAGCATGCTGTTCCTGATCCTTGGCGATAACCTGCTGTTCATC  >  NC_002947/4666364‑4666657

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