Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,747,081 G→A 8.0% S227S (AGC→AGT) ‡ etfB ← electron transfer flavoprotein subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,747,0810GA 8.0% 178.5 / 4.3 63S227S (AGC→AGT) ‡etfBelectron transfer flavoprotein subunit beta
Reads supporting (aligned to +/- strand):  ref base G (26/31);  new base A (3/2);  total (30/33)
Fisher's exact test for biased strand distribution p-value = 6.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

GCGGCTGCAGCGACAGTGTTCAGGGTGGCCGGGGCTACGGCGCCGGCCTCGTATTCAGCGACAACCAGGATAGTCATTTAGATTACCTTCGCTTCGTTCTTCAGCTTCTCGACCAGTTCGGCCACCGATTTGACCTTGATGCCCGCGCTGCGGGCAGCCGGGGCTTCAACCTTCACGGTCTTGTTGGTGGAGGCGAGGGAAACGCCCAGTGCGTCTGGAGTAACAGTCTCCAGCG‑GCTTCTTCTTGGCCTTCATGATGTTCGGCAGCGACGCGTAGCGTGGCTCGT  >  NC_002947/4746935‑4747220
                                                                                                                                                  |                                                                                                                                            
gcggcTGCAGCGACAGTGTTCAGGGTGGCCGGGGCTACGGCGCCGGCCTCGTATTCAGCGACAACCAGGATAGTCATTTAGATTACCTTCGCTTCGTTCTTCAGCTTCTCGACCAGTTCGGCCACCGATTTGACCTTGATGCCCGCGCt                                                                                                                                            >  2:328526/1‑149 (MQ=255)
       cAGCGACAGTGTTCAGGGTGGCCGGGGCTACGGCGCCGGCCTCGTATTCAGCGACAACCAGGATAGTCATTTAGATTACCTTCGCTTCGTTCTTCAGCTTCTCGACCAGTTCGGCCACCGATTTGACCTTGATGCCCGCGCTGCggg                                                                                                                                       >  2:327227/1‑147 (MQ=255)
            aCAGTGTTCAGGGTGGCCGGGGCTACGGCGCCGGCCTCGTATTCAGCGACAACCAGGATAGTCATTTAGATTACCTTCGCTTCGTTCTTCAGCTTCTCGACCAGTTCGGCCACCGATTTGACCTTGATGCCCGCGCTGCGGGCAGCCgg                                                                                                                                >  5:314028/1‑149 (MQ=255)
            aCAGTGTTCAGGGTGGCCGGGGCTACGGCGCCGGCCTCGTATTCAGCGACAACCAGGATAGTCATTTAGATTACCTTCGCTTCGTTCTTCAGCTTCTCGACCAGTTCGGCCACCGATTTGACCTTGATGCCCGCGCTGCGGGCAGCCgg                                                                                                                                >  6:312753/1‑149 (MQ=255)
            aCAGTGTTCAGGGTGGCAGTGGCTACGGCGCAGGCCTCGTATTCAGCGACAACCAGGATAGTCATTTAGATTACCTTCGCTGCGTTCTTCAGCTTCTCGACCAGTTCGGCCACCGATTTGACCTTGATGCCCGCGCTGCGGGCAGCCgg                                                                                                                                <  5:116284/149‑1 (MQ=255)
                  ttCAGGGTGGCCGGGGCTACGGCGCCGGCCTCGTATTCAGCGACAACCAGGATAGTCATTTAGATTACCTTCGCTTCGTTCTTCAGCTTCTCGACCAGTTCGGCCACCGATTTGACCTTGATGCCCGCGCTGCGGGCAGCCGGGGCtt                                                                                                                           <  5:107152/148‑1 (MQ=255)
                    cAGGGTGGCCGGGGCTACGGCGCCGGCCTCGTATTCAGCGACAACCAGGATAGTCATTTAGATTACCTTCGCTTCGTTCTTCAGCTTCTCGACCAGTTCGGCCACCGATTTGACCTTGATGCCCGCGCTGCGGGCAGCCGGGGCTTCaa                                                                                                                        >  6:231834/1‑149 (MQ=255)
                    cAGGGTGGCCGGGGCTACGGCGCCGGCCTCGTATTCAGCGACAACCAGGATAGTCATTTAGATTACCTTCGCTTCGTTCTTCAGCTTCTCGACCAGTTCGGCCACCGATTTGACCTTGATGCCCGCGCTGCGGGCAGCCGGGGCTTCa                                                                                                                         <  5:312753/148‑1 (MQ=255)
                     aGGGTGGCCGGGGCTACGGCGCCGGCCTCGTATTCAGCGACAACCAGGATAGTCATTTAGATTACCTTCGCTTCGTTCTTCAGCTTCTCGACCAGTTCGGCCACCGATTTGACCTTGATGCCCGCAGCGCGGGCatca                                                                                                                                  >  8:88625/1‑135 (MQ=255)
                     aGGGTGGCCGGGGCTACGGCGCCGGCCTCGTATTCAGCGACAACCAGGATAGTCATTTAGATTACCTTCGCTTCGTTCTTCAGCTTCTCGACCAGTTCGGCCACCGATTTGACCTTGATGCCCGCAGCGCGGGCatca                                                                                                                                  <  7:88625/138‑4 (MQ=255)
                      gggTGGCCGGGGCTACGGCGCCGGCCTCGTATTCAGCGACAACCAGGATAGTCATTTAGATTACCTTCGCTTCGTTCTTCAGCTTCTCGACCAGTTCGGCCACCGATTTGACCTTGATGCCCGCGCTGCGGGCAGCCgg                                                                                                                                >  8:122266/1‑139 (MQ=255)
                      gggTGGCCGGGGCTACGGCGCCGGCCTCGTATTCAGCGACAACCAGGATAGTCATTTAGATTACCTTCGCTTCGTTCTTCAGCTTCTCGACCAGTTCGGCCACCGATTTGACCTTGATGCCCGCGCTGCGGGCAGCCgg                                                                                                                                <  7:122266/139‑1 (MQ=255)
                        gTGGCCGGGGCTACGGCGCCGGCCTCGTATTCAGCGACAACCAGGATAGTCATTTAGATTACCTTCGCTTCGTTCTTCAGCTTCTCGACCAGTTCGGCCACCGATTTGACCTTGATGCCCGCGCTGCGGGCAGCCGGGGCTTCAACCtt                                                                                                                    >  5:103066/1‑149 (MQ=255)
                        gTGGCCGGGGCTACGGCGCCGGCCTCGTATTCAGCGACAACCAGGATAGTCATTTAGATTACCTTCGCTTCGTTCTTCAGCTTCTCGACCAGTTCGGCCACCGATTTGACCTTGATGCCCGCGCTGCGGGCAGCCGGGGCTTCAACCtt                                                                                                                    >  7:174143/1‑149 (MQ=255)
                        gTGGCCGGGGCTACGGCGCCGGCCTCGTATTCAGCGACAACCAGGATAGTCATTTAGATTACCTTCGCTTCGTTCTTCAGCTTCTCGACCAGTTCGGCCACCGATTTGACCTTGATGCCCGCGCTGCGGGCAGCCGGGGCTTCAACCt                                                                                                                     >  4:44542/1‑148 (MQ=255)
                                ggCTACGGCGCCGGCCTCGTATTCAGCGACAACCAGGATAGTCATTTAGATTACCTTCGCTTCGTTCTTCAGCTTCTCGACCAGTTCGGCCACCGATTTGACCTTGATGCCCGCGCTGCGGGCAGCCGGGGCTTCAACCTTCACGGTCt                                                                                                            <  2:146890/149‑1 (MQ=255)
                                ggCTACGGCGCCGGCCTCGTATTCAGCGACAACCAGGATAGTCATTTAGATTACCTTCGCTTCGTTCTTCAGCTTCTCGACCAGTTCGGCCACCGATTTGACCTTGATGCCCGCGCTGCGGGCAGCCGGGGCTTCAACCTTCACGGTCt                                                                                                            <  1:151567/149‑1 (MQ=255)
                                     cGGCGCCGGCCTCGTATTCAGCGACAACCAGGATAGTCATTTAGATTACCTTCGCTTCGTTCTTCAGCTTCTCGACCAGTTCGGCCACCGATTTGACCTTGATGCCCGCGCTGCGGGCAGCCGGGGCTTCAACCTTCACGGTCttgt                                                                                                         <  2:182275/147‑1 (MQ=255)
                                      ggCGCCGGCCTCGTATTCAGCGACAACCAGGATAGTCATTTAGATTACCTTCGCTTCGTTCTTCAGCTTCTCGACCAGTTCGGCCACCGATTTGACCTTGATGCCCGCGCTGCGGGCAGCCGGGGCTTCAACCTTCACGGTCTTGTtgg                                                                                                      <  1:327227/149‑1 (MQ=255)
                                         gcCGGCCTCGTATTCAGCGACAACCAGGATAGTCATTTAGATTACCTTCGCTTCGTTCTTCAGCTTCTCGACCAGTTCGGCCACCGATTTGACCTTGATGCCCGCGCTGCGGGCAGCCGGGGCTTCAACCTTCACGGTCTTGTtggtg                                                                                                    >  8:304243/1‑148 (MQ=255)
                                          ccGGCCTCGTATTCAGCGACAACCAGGATAGTCATTTAGATTACCTTCGCTTCGTTCTTCAGCTTCTCGACCAGTTCGGCCACCGATTTGACCTTGATGCCCGCGCTGCGGGCAGCCGGGGCTTCAACCTTCACGGTCTTGTTGGTGGa                                                                                                  <  3:59993/149‑1 (MQ=255)
                                           cGGCCTCGTATTCAGCGACAACCAGGATAGTCATTTAGATTACCTTCGCTTCGTTCTTCAGCTTCTCGACCAGTTCGGCCACCGATTTGACCTTGATGCCCGCGCTGCGGGCAGCCGGGGCTTCAACCTTCACGGTCTTGTTGGTGGAg                                                                                                 <  7:178122/149‑1 (MQ=255)
                                             gCCTCGTATTCAGCGACAACCAGGATAGTCATTTAGATTACCTTCGCTTCGTTCTTCAGCTTCTCGACCAGTTCGGCCACCGATTTGACCTTGATGCCCGCGCTGCGGGCAGCCGGGGCTTCAACCTTCACGGTCTTGTTGGTGGAGGc                                                                                               <  3:44542/149‑1 (MQ=255)
                                               cTCGTATTCAGCGACAACCAGGATAGTCATTTAGATTACCTTCGCTTCGTTCTTCAGCTTCTCGACCAGTTCGGCCACCGATTTGACCTTGATGCCCGCGCTGCGGGCAGCCGGGGCTTCAACCTTCACGGTCTTGTTGGTGGAGGCGa                                                                                             <  6:314028/149‑1 (MQ=255)
                                               cTCGTATTCAGCGACAACCAGGATAGTCATTTAGATTACCTTCGCTTCGTTCTTCAGCTTCTCGACCAGTTCGGCCACCGATTTGACCTTGATGCCCGCGCTGCGGGCAGCCGGGGCTTCAACCTTCACGGTCTTGTTGGTGGAGGCGa                                                                                             <  7:304243/149‑1 (MQ=255)
                                                        aGCGACAACCAGGATAGTCATTTAGATTACCTTCGCTTCGTTCTTCAGCTTCTCGACCAGTTCGGCCACCGATTTGACCTTGATGCCCGCGCTGCGGGCAGCCGGGGCTTCAACCTTCACGGTCTTGTTGGTGGAGGCGAGGGAAACGc                                                                                    >  8:192904/1‑149 (MQ=255)
                                                           gACAACCAGGATAGTCATTTAGATTACCTTCGCTTCGTTCTTCAGCTTCTCGACCAGTTCGGCCACCGATTTGACCTTGATGCCCGCGCTGCGGGCAGCCGGGGCTTCAACCTTCACGGTCTTGTTGGTGGAGGCGAGGGAAACGCCCa                                                                                 <  3:299950/149‑1 (MQ=255)
                                                           gACAACCAGGATAGTCATTTAGATTACCTTCGCTTCGTTCTTCAGCTTCTCGACCAGTTCGGCCACCGATTTGACCTTGATGCCCGCGCTGCGGGCAGCCGGGGCTTCAACCTTCACGGTCTTGTTGGTGGAGGCGAGGGAAACGCCCa                                                                                 <  5:231834/149‑1 (MQ=255)
                                                            aCAACCAGGATAGTCATTTAGATTACCTTCGCTTCGTTCTTCAGCTTCTCGACCAGTTCGGCCACCGATTTGACCTTGATGCCCGCGCTGCGGGCAGCCGGGGCTTCAACCTTCACGGTCTTGTTGGTGGAGGCGAGGGAAACGCCCAg                                                                                >  3:158708/1‑149 (MQ=255)
                                                                ccAGGATAGTCATTTAGATTACCTTCGCTTCGTTCTTCAGCTTCTCGACCAGTTCGGCCACCGATTTGACCTTGATGCCCGCGCTGCGGGCAGCCGGGGCTTCAACCTTCACGGTCTTGTTGGTGGAGGCGAGGGAAACGCCCAGTGCg                                                                            >  4:117621/1‑149 (MQ=255)
                                                                    gATAGTCATTTAGATTACCTTCGCTTCGTTCTTCAGCTTCTCGACCAGTTCGGCCACCGATTTGACCTTGATGCCCGCGCTGCGGGCAGCCGGGGCTTCAACCTTCACGGTCTTGTTGGTGGAGGCGAGGGAAACGCCCAGTGCGTCTg                                                                        <  4:158708/149‑1 (MQ=255)
                                                                           aTTTAGATTACCTTCGCTTCGTTCTTCAGCTTCTCGACCAGTTCGGCCACCGATTTGACCTTGATGCCCGCGCTGCGGGCAGCCGGGGCTTCAACCTTCACGGTCTTGTTGGTGGAGGCGAGGGAAACGCCCAGTGCGTCTGGAGTaa                                                                  <  1:328526/148‑1 (MQ=255)
                                                                            tttAGATTACCTTCGCTTCGTTCTTCAGCTTCTCGACCAGTTCGGCCACCGATTTGACCTTGATGCCCGCGCTGCGGGCAGCCGGGGCTTCAACCTTCACGGTCTTGTTGGTGGAGGCGAGGGAAACGCCCAGTGCGTCTGGa                                                                      <  3:21162/143‑1 (MQ=255)
                                                                            tttAGATTACCTTCGCTTCGTTCTTCAGCTTCTCGACCAGTTCGGCCACCGATTTGACCTTGATGCCCGCGCTGCGGGCAGCCGGGGCTTCAACCTTCACGGTCTTGTTGGTGGAGGCGAGGGAAACGCCCAGTGCGTCTGGa                                                                      >  4:21162/1‑143 (MQ=255)
                                                                            tttAGATTACCTTCGCTTCGTTCTTCAGCTTCTCGACCAGTTCGGCCACCGATTTGACCTTGATGCCCGCGCTGCGGGCAGCCGGGGCTTCAACCTTCACGGTCTTGTTGGTGGAGGCGAGGGAAACGCCCAGTGCGTCTGGAGTaa                                                                  >  6:305656/1‑147 (MQ=255)
                                                                             ttAGATTACCTTCGCTTCGTTCTTCAGCTTCTCGACCAGTTCGGCCACCGATTTGACCTTGATGCCCGCGCTGCGGGCAGCCGGGGCTTCAACCTTCACGGTCTTGTTGGTGGAGGCGAGGGAAACGCCCAGTGCGTCTgg                                                                       >  5:211466/1‑141 (MQ=255)
                                                                             ttAGATTACCTTCGCTTCGTTCTTCAGCTTCTCGACCAGTTCGGCCACCGATTTGACCTTGATGCCCGCGCTGCGGGCAGCCGGGGCTTCAACCTTCACGGTCTTGTTGGTGGAGGCGAGGGAAACGCCCAGTGCGTCTgg                                                                       >  5:206690/1‑141 (MQ=255)
                                                                             ttAGATTACCTTCGCTTCGTTCTTCAGCTTCTCGACCAGTTCGGCCACCGATTTGACCTTGATGCCCGCGCTGCGGGCAGCCGGGGCTTCAACCTTCACGGTCTTGTTGGTGGAGGCGAGGGAAACGCCCAGTGCGTCTgg                                                                       <  6:206690/141‑1 (MQ=255)
                                                                             ttAGATTACCTTCGCTTCGTTCTTCAGCTTCTCGACCAGTTCGGCCACCGATTTGACCTTGATGCCCGCGCTGCGGGCAGCCGGGGCTTCAACCTTCACGGTCTTGTTGGTGGAGGCGAGGGAAACGCCCAGTGCGTCTgg                                                                       <  6:211466/141‑1 (MQ=255)
                                                                                     ccTTCGCTTCGTTCTTCAGCTTCTCGACCAGTTCGGCCACCGATTTGACCTTGATGCCCGCGCTGCGGGCAGCCGGGGCTTCAACCTTCACGGTCTTGTTGGTGGAGGCGAGGGAAACGCCCAGTGCGTCTGGAGTAACAGTCTCCAGc                                                       <  6:103066/149‑1 (MQ=255)
                                                                                      cttcgcttcgTTCTTCAGCTTCTCGACCAGTTCGGCCACCGATTTGACCTTGATGCCCGCGCTGCGGGCAGCCGGGGCTTCAACCTTCACGGTCTTGTTGGTGGAGGCGAGGGAAACGCCCAGTGCGTCTGGAg                                                                     <  3:91211/134‑1 (MQ=255)
                                                                                      cttcgcttcgTTCTTCAGCTTCTCGACCAGTTCGGCCACCGATTTGACCTTGATGCCCGCGCTGCGGGCAGCCGGGGCTTCAACCTTCACGGTCTTGTTGGTGGAGGCGAGGGAAACGCCCAGTGCGTCTGGAg                                                                     >  4:91211/1‑134 (MQ=255)
                                                                                                ttcttcAGCTTCTCGACCAGTTCGGCCACCGATTTGACCTTGATGCCCGCGCTGCGGGCAGCCGGGGCTTCAACCTTCACGGTCTTGTTGGTGGAGGCGAGGGAAACGCCCAGTGCGTCTGGAGTAACAGTCTCCAGCG‑Gcttcttctt                                           <  3:117621/149‑1 (MQ=255)
                                                                                                 tcttcAGCTTCTCGACCAGTTCGGCCACCGATTTGACCTTGATGCCCGCGCTGCGGGCAGCCGGGGCTTCAACCTTCACGGTCTTGTTGGTGGAGGCGAGGGAAACGCCCAGTGCGTCTGGAGTAAc                                                                 <  7:79401/127‑1 (MQ=255)
                                                                                                 tcttcAGCTTCTCGACCAGTTCGGCCACCGATTTGACCTTGATGCCCGCGCTGCGGGCAGCCGGGGCTTCAACCTTCACGGTCTTGTTGGTGGAGGCGAGGGAAACGCCCAGTGCGTCTGGAGTAAc                                                                 >  8:79401/1‑127 (MQ=255)
                                                                                                       gCTTCTCGACCAGTTCGGCCACCGATTTGACCTTGATGCCCGCGCTGCGGGCAGCCGGGGCTTCAACCTTCACGGTCTTGTTGGTGGAGGCGAGGGAAACGCCCAGTGCGTCTGGAGTAACAGTCTCCAGCG‑GCTTCTTCTTGGCCTTc                                    <  7:192904/149‑1 (MQ=255)
                                                                                                           ctcGACCAGTTCGGCCACCGATTTGACCTTGATGCCCGCGCTGCGGGCAGCCGGGGCTTCAACCTTCACGGTCTTGTTGGTGGAGGCGAGGGAAACGCCCAGTGCGTCTGGAGTAACAGTCTCCAGCG‑GCTTCTTCTTGGCCTTCatga                                >  1:34210/1‑149 (MQ=255)
                                                                                                             cGACCAGTTCGGCCACCGATTTGACCTTGATGCCCGCGCTGCGGGCAGCCGGGGCTTCAACCTTCACGGTCTTGTTGGTGGAGGCGAGGGAAACGCCCAGTGCGTCTGGAGTAACAGTCTCCAGCG‑GCTTCTTCTTGGCCTTCatgatg                              >  7:91340/1‑149 (MQ=255)
                                                                                                              gACCAGTTCGGCCCCCGATTTGACCTTGATGCCCGCGCTGCGGGCAGCCGGGGCTTCAACCTTCACGGTCTTGTTGGTGGAGGCGAGGGAAACGCCCAGTGCGTCTGGAGTAACAGTCTCCAGCG‑GCTTCTTCTTGGCCTTCatga                                <  4:126634/146‑1 (MQ=255)
                                                                                                                  aGTTCGGCCACGGATTTGACCTTGATGCCCGCGCTGCGGGCAGCCGGGGCTTCAAGCTTCACGGTGTTGTTGGTGGAGGCGTGGGAAACGCCCAGTGCGTCTGGAGTAACAGTCTCCAGCG‑GCTTCTTCTTGGCCTTCATGATGTTCgg                         <  7:183474/149‑1 (MQ=255)
                                                                                                                  aGTTCGGCCACCGATTTGACCTTGATGCCCGCGCTGCGGGCAGCCGGGGCTTCAACCTTCACGGTCTTGTTGGTGGAGGCGAGGGAAACGCCCAGTGCGTCTGGAGTAACAGTCTCCAGCG‑GCTTCTTCTTGGCCTTCATGATGTTCgg                         <  8:91340/149‑1 (MQ=255)
                                                                                                                       ggCCACCGATTTGACCTTGATGCCCGCGCTGCGGGCAGCCGGGGCTTCAACCTTCACGGTCTTGTTGGTGGAGGCGAGGGAAACGCCCAGTGCGTCTGGAGTAACAGTCTCCAGCG‑GCTTCTTCTTGGCCTTCATGATGTTCGGCAGc                     <  5:305656/148‑1 (MQ=255)
                                                                                                                         ccACCGATTTGACCTTGATGCCCGCGCTGCGGGCAGCCGGGGCTTCAACCTTCACGGTCTTGTTGGTGGAGGCGAGGGAAACGCCCAGTGCGTCTGGAGTAACAGTCTCCAGCG‑GCTTCTTCTTGGCCTTCATGATGTTCGGCAGCGAc                  >  8:68433/1‑149 (MQ=255)
                                                                                                                          cACCGATTTGACCTTGATGCCCGCGCTGCGGGCAGCCGGGGCTTCAACCTTCACGGTCTTGTTGGTGGAGGCGAGGGAAACGCCCAGTGCGTCTGGAGTAACAGTCTCCAGCG‑GCTTCTTCTTGGCCTTCATGATGTTCGGCAGCGAcg                 >  8:14071/1‑149 (MQ=255)
                                                                                                                           aCCGATTTGACCTTGATGCCCGCGCTGCGGGCAGCCGGGGCTTCAACCTTCACGGTCTTGTTGGTGGAGGCGAGGGAAACGCCCAGTGCGTCTGGAGTAACAGTCTCCAGCG‑Gcttcttctt                                           >  4:25001/1‑122 (MQ=255)
                                                                                                                           aCCGATTTGACCTTGATGCCCGCGCTGCGGGCAGCCGGGGCTTCAACCTTCACGGTCTTGTTGGTGGAGGCGAGGGAAACGCCCAGTGCGTCTGGAGTAACAGTCTCCAGCG‑Gcttcttctt                                           <  3:25001/122‑1 (MQ=255)
                                                                                                                               aTTTGACCTTGATGCCCGCGCTGCGGGCAGCCGGGGCTTCAACCTTCACGGTCTTGTTGGTGGAGGCGAGGGAAACGCCCAGTGCGTCTGGAGTAACAGTCTCCAGCG‑GCTTCTTCTTGGCCTTCATGATGTTCGGCAGCGACGCGTAg            >  5:61032/1‑149 (MQ=255)
                                                                                                                                    aCCTTGATGCCCGCGCTGCGGGCAGCCGGGGCTTCAACCTTCACGGTCTTGTTGGTGGAGGCGAGGGAAACGCCCAGTGCGTCTGGAGTAACAGTCTCCAGCG‑GCTTCTTCTTGGCCTTCATGATGTTCGGCAGCGACGCGTAGCGTgg       <  7:68433/149‑1 (MQ=255)
                                                                                                                                     cccggctGCCCGCAGCGCGGGCATCCGGGGCTTCAACCTTCACGGTCTTGTTGGTGGAGGCGAGGGAAACGCCCAGTGCGTCTGGAGTAACAGTCTCCAGCG‑GCTTCTTCTTGGCCTTCATGATGTTCGGCAGCGACGCGTAGCGTGGc      >  2:258887/7‑149 (MQ=255)
                                                                                                                                         gctGCCCGCAGCGCGGGCATCCGGGGCTTCAACCTTCACGGTCTTGTTGGTGGAGGCGAGGGAAACGCCCAGTGCGTCTGGAGTAACAGTCTCCAGCGCGCTTCTTCTTGGCCTTCATGATGTTCGGCAGCGACGCGTAGCGTGGCTCg   <  1:258887/147‑1 (MQ=255)
                                                                                                                                         gATGCCCGCCCTGCGGGCAGCCGGGGCTTCAACCTTCACGGTCTTGTTGGTGGAGGCGAGGGAAACGCCCACTGCGTCTGGAGTAACAGTCTCCAGCG‑GCTTCTTCTTCGCCTTCATGAAGTTCCGAAGCGACGCGTAGCGTGGCTCGt  >  4:152106/1‑149 (MQ=255)
                                                                                                                                          ctGCCCGCAGCGCGGGCAGCCGGGGCTTCAACCTTCACGGTCTTGTTGGTGGAGGCGAGGGAAACGCCCAGTGCGTCTGGAGTAACAGTCTCCAGCG‑GCTTCTTCTTGGCCTTCATGATGTTCgg                         >  5:46814/2‑125 (MQ=255)
                                                                                                                                             cccGCGCTGCGGGCAGCCGGGGCTTCAACCTTCACGGTCTTGTTGGTGGAGGCGAGGGAAACGCCCAGTGCGTCTGGAGTAACAGTCTCCAGCG‑GCTTCTTCTTGGCCTTCatgatg                              >  6:206622/1‑117 (MQ=255)
                                                                                                                                                  |                                                                                                                                            
GCGGCTGCAGCGACAGTGTTCAGGGTGGCCGGGGCTACGGCGCCGGCCTCGTATTCAGCGACAACCAGGATAGTCATTTAGATTACCTTCGCTTCGTTCTTCAGCTTCTCGACCAGTTCGGCCACCGATTTGACCTTGATGCCCGCGCTGCGGGCAGCCGGGGCTTCAACCTTCACGGTCTTGTTGGTGGAGGCGAGGGAAACGCCCAGTGCGTCTGGAGTAACAGTCTCCAGCG‑GCTTCTTCTTGGCCTTCATGATGTTCGGCAGCGACGCGTAGCGTGGCTCGT  >  NC_002947/4746935‑4747220

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