Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 82,857 G→T 8.2% P91Q (CCG→CAG)  qor‑I ← quinone oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294782,8570GT 8.2% 178.0 / 5.4 61P91Q (CCG→CAG) qor‑Iquinone oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (24/32);  new base T (2/3);  total (26/35)
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.46e-01

CTTCACCTCATAGGTCTGCTTCAGCAGGTACTGTACGGTCAAACCCTTGAGCATCACCGCTGCCGCCTGCTCGAAGCTGATGGTGTCCGGCAGCTTGACCAGGTTGGCTTCCGGTAGCGTATGCACCTCGCTGTAGGCGCCCAGCGGGCCGCCGGCATGTGCCACACGGTCTCCTACCTTCAGCCGGGACACGCCTTCGCCCACGGCCTCTACCACACCTGCCGCTTCGGTACCCAGGCCAGAAGGCAGGGACGGCGGCGCATACAGCCCGCTTCGGAAGTAGGTATCGAT  >  NC_002947/82712‑83002
                                                                                                                                                 |                                                                                                                                                  
cTTCACCTCATAGGTCTGCTTCAGCAGGTACTGTACGGTCAAACCCTTGAGCATCACCGCTGCCGCCTGCTCGAAGCTGATGGTGTCCGGCAGCTTGACCAGGTTGGCTTCCGGTAGCGTATGCACCTCGCTGTAGGCGCCCAGCGGgc                                                                                                                                                 >  7:289183/1‑149 (MQ=255)
cTTCACCTCATAGGTCTGCTTCAGCAGGTACTGTACGGTCAAACCCTTGAGCATCACCGCTGCCGCCTGCTCGAAGCTGATGGTGTCCGGCAGCTTGACCAGGTTGGCTTCCGGTAGCGTATGCACCTCGCTGTAGGCGCCCAGCGGgc                                                                                                                                                 <  7:4361/149‑1 (MQ=255)
 ttCACCTCATAGGTCTGCTTCAGCAGGTACTGTACGGTCAAACCCTTGAGCATCACCGCTGCCGCCTGCTCGAAGCTGATGGTGTCCGGCAGCTTGACCAGGTTGGCTTCCGGTAGCGTATGCACCTCGCTGTAGGCGCCCAGCGGgcc                                                                                                                                                >  1:174968/1‑149 (MQ=255)
      cTCATAGGTCTGCTTCAGCAGGTACTGTACGGTCAAACCCTTGAGCATCACCGCTGCCGCCTGCTCGAAGCTGATGGTGTCCGGCAGCTTGACCAGGTTGGCTTCCGGTAGCGTATGCACCTCGCTGTAGGCGGCCCGCTGggcgccta                                                                                                                                           <  1:336063/149‑3 (MQ=255)
            ggTCTGATTCAGCAGGTACTGTACGGTCAAACCCTTGAGCATCCCCGCTGCCCCCTGCTCGAAGCTGATGGTGGCCGGCAGCGTGCCCAGGTTGGCTTCCGGTAGCGTATGCACCTCGCTCTAGGCGCCCAGCGCGCCGCCGGCAtgtg                                                                                                                                    >  2:336409/1‑149 (MQ=255)
              tCTGCTTCAGCAGGTACTGTACGGTCAAACCCTTGAGCATCACCGCTGCCGCCTGCTCGAAGCTGATGGTGTCCGGCAGCTTGACCAGGTTGGCTTCCGGTAGCGTATGCACCTCGCTGTAGGCGCCCAGCGGGCCGCCGGCATGTGcc                                                                                                                                  >  6:204779/1‑149 (MQ=255)
              tCTGCTTCAGCAGGTACTGTACGGTCAAACCCTTGAGCATCACCGCTGCCGCCTGCTCGAAGCTGATGGTGTCCGGCAGCTTGACCAGGTTGGCTTCCGGTAGCGTATGCACCTCGCTGTAGGCGCCCAGCGGGCCGCCGGCATGTGc                                                                                                                                   <  7:317816/148‑1 (MQ=255)
                  cTTCAGCAGGTACTGTACGGTCAAACCCTTGAGCATCACCGCTGCCGCCTGCTCGAAGCTGATGGTGTCCGGCAGCTTGACCAGGTTGGCTTCCGGTAGCGTATGCACCTCGCTGTAGGCGCCCAGCGGGCCGCCGGCATGTGCcaca                                                                                                                               <  3:186508/148‑1 (MQ=255)
                  cTTCAGCAGGTACTGTACGGTCAAACCCTTGAGCATCACCGCTGCCGCCTGCTCGAAGCTGATGGTGTCCGGCAGCTTGACCAGGTTGGCTTCCGGTAGCGTATGCACCTCGCTGTAGGCGCCCAGCGGGCCGCCGGCATGTGCcac                                                                                                                                <  5:336376/147‑1 (MQ=255)
                  cTTCAGCAGGTACTGTACGGTCAAACCCTTGAGCATCACCGCTGCCGCCTGCTCGAAGCTGATGGTGTCCGGCAGCTTGACCAGGTTGGCTTCCGGTAGCGTATGCACCTCGCTGTAGGCGCCCAGCGGGCCGCCGGCATGTGCCacac                                                                                                                              >  1:46174/1‑149 (MQ=255)
                     cagcagGTACTGTACGGTCAAACCCTTGAGCATCACCGCTGCCGCCTGCTCGAAGCTGATGGTGTCCGGCAGCTTGACCAGGTTGGCTTCCGGTAGCGTATGCACCTCGCTGTAGGCGCCCAGCGGGCCGCCGGCATGTGCCACACg                                                                                                                             <  1:63083/147‑1 (MQ=255)
                      agcagGTACTGTACGGTCAAACCCTTGAGCATCACCGCTGCCGCCTGCTCGAAGCTGATGGTGTCCGGCAGCTTGACCAGGTTGGCTTCCGGTAGCGTATGCACCTCGCTGTAGGCGCCCAGCGGGCCGCCGGCATGTGc                                                                                                                                   >  3:109321/1‑140 (MQ=255)
                      agcagGTACTGTACGGTCAAACCCTTGAGCATCACCGCTGCCGCCTGCTCGAAGCTGATGGTGTCCGGCAGCTTGACCAGGTTGGCTTCCGGTAGCGTATGCACCTCGCTGTAGGCGCCCAGCGGGCCGCCGGCATGTGc                                                                                                                                   <  4:109321/140‑1 (MQ=255)
                         agGTACTGTACGGTCAAACCCTTGATCATCACCGCTGCCGCCTGCTCGAAGCTGATGGTGTCCGGCAGCTTGACCAGGTTGGCTTCCGGTAGCGTATGCACCTCGCTGTAGGCGCCCAGCGGGCCGCCGGCATGTGCCACACGGTCTcc                                                                                                                       <  4:256423/149‑1 (MQ=255)
                           gTACTGTACGGTCAAACCCTTGAGCATCACCGCTGCCGCCTGCTCGAAGCTGATGGTGTCCGGCAGCTTGACCAGGTTGGCTTCCGGTAGCGTATGCACCACGCTGTAGGCGGCCCGCTGggcgcctaca                                                                                                                                         >  7:20307/1‑126 (MQ=255)
                           gTACTGTACGGTCAAACCCTTGAGCATCACCGCTGCCGCCTGCTCGAAGCTGATGGTGTCCGGCAGCTTGACCAGGTTGGCTTCCGGTAGCGTATGCACCACGCTGTAGGCGGCCCGCTGggcgcctaca                                                                                                                                         <  8:20307/130‑5 (MQ=255)
                              cTGTACGGTCAAACCCTTGAGCATCACCGCTGCCGCCTGCTCGAAGCTGATGGTGTCCGGCAGCTTGACCAGGTTGGCTTCCGGTAGCGTATGCACCTCGCTGTAGGCGCCCAGCGGGCCGCCGGCATGTGCCACACGGTCTCCTACCt                                                                                                                  <  5:71998/149‑1 (MQ=255)
                                gTACGGTCAAACCCTTGAGCATCACCGCTGCCGCCTGCTCGAAGCTGATGGTGTCCGGCAGCTTGACCAGGTTGGCTTCCGGTAGCGTATGCACCTCGCTGTAGGCGCCCAGCGGGCCGCCGGCATGTGCCACACGGTCTCCTACCtt                                                                                                                 >  3:274411/1‑148 (MQ=255)
                                 tACGGTCAAACCCTTGAGCATCACCGCTGCCGCCTGCTCGAAGCTGATGGTGTCCGGCAGCTTGACCAGGTTGGCTTCCGGTAGCGTATGCACCTCGCTGTAGGCGCCCAGCg                                                                                                                                                    >  6:188922/1‑113 (MQ=255)
                                   cGGTCAAACCCTTGAGCATCACCGCTGCCGCCTGCTCGAAGCTGATGGTGTCCGGCAGCTTGACCAGGTTGGCTTCCGGTAGCGTATGCACCTCGCTGTAGGCGCCCAGCGGGCCGCCGGCATGTGCCACACGGTCTCCTACCTTCAg                                                                                                              <  4:274411/148‑1 (MQ=255)
                                     ggccaaCCCTTGTGCATCACCGCTGCCGCCTGCTCGATGCTGATGGTGTCCTGCAGCTTGACCAGGTTGGCTTCCGGTAGCGTATGCACCTCGCTGTAGGCGCCCAGCGGGCCGCCGGCATGTGCCACACGGTCTCCTACCTTCAGCCg                                                                                                           <  6:257538/145‑1 (MQ=255)
                                      tCAAACCCTTGAGCATCACCGCTGCCGCCTGCTCGAAGCTGATGGTGTCCGGCAGCTTGACCAGGTTGGCTTCCGGTAGCGTATGCACCTCGCTGTAGGCGCCCAGCGGGCCGCCGGCATGTGCCACACGGTCTCCTACCTTCAGCCgg                                                                                                          >  3:385968/1‑149 (MQ=255)
                                          aCCCTTGAGCATCACCGCTGCCGCCTGCTCGAAGCTGATGGTGTCCGGCAGCTTGACCAGGTTGGCTTCCGGTAGCGTATGCACCTCGCTGTAGGCGCCCAGCGGGCCGCCGGCATGTGCCACACGGTCTCCTACCTTCAGCCGGGaca                                                                                                      <  2:304816/149‑1 (MQ=255)
                                                     tCACCGCTGCCGCCTGCTCGAAGCTGATGGTGTCCGGCAGCTTGACCAGGTTGGCTTCCGGTAGCGTATGCACCTCGCTGTAGGCGCCCAGCGGGCCGCCGGCAtgtg                                                                                                                                    <  5:76897/108‑1 (MQ=255)
                                                     tCACCGCTGCCGCCTGCTCGAAGCTGATGGTGTCCGGCAGCTTGACCAGGTTGGCTTCCGGTAGCGTATGCACCTCGCTGTAGGCGCCCAGCGGGCCGCCGGCAtgtg                                                                                                                                    >  6:76897/1‑108 (MQ=255)
                                                       aCCGCTGCCGCCTGCTTGAAGCTGATGGTGTCCGGCAGCTTGACCAGGTTGGCTTCCGGTAGCGTATGCACCTCGCTGTAGGCGCCCAGCGGGCCGCCGGCATGTTGCCACACGGTCTCCTACCTTCAGCCGGGACACGCCTTCGCCCa                                                                                          <  1:81235/149‑1 (MQ=255)
                                                       aCCGCTGCCGCCTGCTCGAAGCTGATGGTGTCCGGCAGCTTGACCAGGTTGGCTTCCGGTAGCGTATGCACCTCGCTGTAGGCGCCCAGCGGGCCGCCGGCATGTGCCACACGGTCTCCTACCTTCAGCCGGGACACGCCTTCGCCCa                                                                                          <  8:129959/148‑1 (MQ=255)
                                                              ccgccTGCTCGAAGCTGATGGTGTCCGGCAGCTTGACCAGGTTGGCTTCCGGTAGCGTATGCACCTCGCTGTAGGCGCCCAGCGGGCCGCCGGCATGTGCCACACGGTCTCCTACCTTCAGCCGGGACACGCCTTCGCCCACGGCctc                                                                                   <  4:385968/148‑1 (MQ=255)
                                                                        gAAGCTGATGGTGTCCGGCAGCTTGACCAGGTTGGCTTCCGGTAGCGTATGCACCTCGCTGTAGGCGCCCAGCGGGCCGCCGGCATGTGCCAAACGGTCTCCTACCTTCAGCCGGGACACGCCTTCGCCCACGGCCTCTACCACACCTg                                                                        >  1:230778/1‑149 (MQ=255)
                                                                          aGCTGATGGTGTCCGGCAGCTTGACCAGGTTGGCTTCCGGTAGCGTATGCACCTCGCTGTAGGCGCCCAGCGGGCCGCCGGCATGTGCCACACGGGCTCCTACCTTCAGCCGGGACACGCCTTCGCCCACGGCCTCTACCACACCTGcc                                                                      <  8:252421/149‑1 (MQ=255)
                                                                          aGCTGATGGTGTCCGGAAGCTTGACCAGGTTGGCTTCCGGTAGCGTATGCACCTCGCTGTAGGCGCCCAGCGGGCCGCCGGCATGTGCCACACGGTCTCCTACCTTCAGCCGGGACACGCCTTCGCCCACGGCCTCTACCACACCTGcc                                                                      <  8:252418/149‑1 (MQ=255)
                                                                             tGATGGTGTCCGGCAGATTGACCAGGTTGGCTTCCGGTAGCTTATGCACCTCGCTGTAGGCGCCCAGCGGGACGCCGTCATGTGCGACACTGTCTCCTACCTTCAGCCGGTAAACGCCTTc                                                                                               <  8:284505/121‑1 (MQ=255)
                                                                              gATGGTGTCCGGCAGCTTGACCAGGTTGGCTTCCGGTAGCGTATGCACCTCGCTGTAGGCGCCCAGCGGGCCGCCGGCATGTGCCACACGGTCTCCTACCTTCAGCCGGGACACGCCTTCGCCCACGGCCTCTACCACACCTGCCGCtt                                                                  >  4:155283/1‑149 (MQ=255)
                                                                               aTGGTGTCCGGCAGCTTGACCAGGTTGGCTTCCGGTAGCGTATGCACCTCGCTGTAGGCGCCCAGCGGGCCGCCGGCATGTGCCACACGGTCTCCTACCTTCAGCCGGGACACGCCTTCGCCCACGGCCTCTACCACACCTGCCGCtt                                                                  <  2:174968/148‑1 (MQ=255)
                                                                                   tgtCCGGCAGCTTGACCAGGTTGGCTTCCGGTAGCGTATGCACCTCGCTGTAGGCGCCCAGCGGGCCGCCGGCATGTGCCACACGGTCTCCTACCTTCAGCCGGGACACGCCTTCGCCCACGGCCTCTACCACACCTGCCGCTTCGGTa                                                             <  5:204779/149‑1 (MQ=255)
                                                                                    gtCCGGCAGCTTGACCAGGTTGGCTTCCGGTAGCGTATGCACCTCGCTGTAGGCGCCCAGCGGGCCGCCGGCATGTGCCACACGGTCTCCTACCTTCAGCCGGGACACGCCTTCGCCCACGGCCTCTACcaca                                                                            <  1:178125/133‑1 (MQ=255)
                                                                                    gtCCGGCAGCTTGACCAGGTTGGCTTCCGGTAGCGTATGCACCTCGCTGTAGGCGCCCAGCGGGCCGCCGGCATGTGCCACACGGTCTCCTACCTTCAGCCGGGACACGCCTTCGCCCACGGCCTCTACcaca                                                                            >  2:178125/1‑133 (MQ=255)
                                                                                             cTTGACCAGGTTGGCTTCCGGTAGCGTATGCACCTCGCTGTAGGCGCCCAGCGGGCCGCCGGCATGTGCCACACGGTCTCCTACCTTCAGCCGGGACACGCCTTCGCCCACGGCCTCTACCACACCTGCCGCTTCGGTACCCAGGc                                                      >  5:208453/1‑146 (MQ=255)
                                                                                             cTTGACCAGGTTGGCTTCCGGTAGCGTATGCACCTCGCTGTAGGCGCCCAGCGGGCCGCCGGCATGTGCCACACGGTCTCCTACCTTCAGCCGGGACACGCCTTCGCCCACGGCCTCTACCACACCTGCCGCTTCGGTACCCAGGc                                                      >  4:86981/1‑146 (MQ=255)
                                                                                               tGACCAGGTTGGCTTCCGGTAGCGTATGCACCTCGCTGTAGGCGCCCAGCGGGCCGCCCGCATGTGCCACACGCTATCCTACATTCAGCCGGGACACGCCTTCCCCCACGCCCTCTAACACACCTGCCGCTGCGCTACCCAGGCCac                                                   >  8:274095/1‑146 (MQ=255)
                                                                                                       ttGGCTTCCGGTAGCGTATGCACCTCGCTGTAGGCGCCCAGCGGGCCGCCGGCATGTGCCACACGGTCTCCTACCTTCAGCCGGGACACGCCTTCGCCCACGGCCTCTACCACACCTGCCGCTTCGGTACCCAGGCCAGAAGGCAGGGa                                         <  2:230778/149‑1 (MQ=255)
                                                                                                            ttCCGGTAGCGTATGCACCTCGCTGTAGGCGCCCAGCGGGCCGCCGGCATGTGCCACACGGTCTCCTACCTTCAGCCGGGACACGCCTTCGCCCACGGCCTCTACCACACCTGCCGCTTCGGTACCCAGGCCAGAAGGCAGGGAcggcg                                    >  8:202763/1‑149 (MQ=255)
                                                                                                               cGGTAGCGTATGCACCTCGCTGTAGGCGCCCAGCGGGCCGCCGGCATGTGCCACACGGTCTCCTACCTTCAGCCGGGACACGCCTTCGCCCACGGCCTCTACCACACCTGCCGCTTCGGTACCCAGGCCAGAAGGCAGGGACggcggc                                  <  7:202763/148‑1 (MQ=255)
                                                                                                               cGGTAGCGTATGCACCTCGCTGTAGGCGCCCAGCGGGCCGCCGGCATGTGCCACACGGTCTCCTACCTTCAGCCGGGACACGCCTTCGCCCACGGCCTCTACCACACCTGCCGCTTCGGTACCCAGGCCAGAAGGCAGGGACggcggc                                  <  2:46174/148‑1 (MQ=255)
                                                                                                                ggTAGCGTATGCACCTCGCTGTAGGCGCCCAGCGGGCCGCCGGCATGTGCCACACGGTCTCCTACCTTCAGCCGGGACACGCCTTCGCCCACGGCCTCTACCACACCTGCCGCTTCGGTAccc                                                          >  4:286236/1‑123 (MQ=255)
                                                                                                                ggTAGCGTATGCACCTCGCTGTAGGCGCCCAGCGGGCCGCCGGCATGTGCCACACGGTCTCCTACCTTCAGCCGGGACACGCCTTCGCCCACGGCCTCTACCACACCTGCCGCTTCGGTAccc                                                          <  3:286236/123‑1 (MQ=255)
                                                                                                                           cACCCCGCTGTAGGCGCCCAGCGGGCCGCCTGCATGTGCCACACGGTCTCCTACCTTCAGCCGGGACACGCCTTCGCCCACGGCCTCTACCACACCTGCCGCTTCGGTACCCAGGCCAGAAGGCAGGGACGGCGGCGCATACAGCCCg                      <  2:384382/148‑1 (MQ=255)
                                                                                                                            aCCTCGCTGTAGGCGCCCAGCGGGCCGCCGGCATGTGCCACACGGTCTCCTACCTTCAGCCGGGACACGCCTTCGCCCACGGCCTCTACCACACCTGCCGCTTCGGTACCCAGGCCAGAAGGCAGGGACGGCGGCGCATACAGcgaggt                    >  6:326481/1‑144 (MQ=255)
                                                                                                                             ccTCGCTGTAGGCGCCCAGCGGGCCGCCGGCATGTGCCACACGGTCTCCTACCTTCAGCCGGGACACGCCTTCGCCCACGGCCTCTACCACACCTGCCGCTTCGGTACCCAGGCCAGAAGGCAGGGACGGCGGCGCATACAGCCCGCtt                   <  6:208453/149‑1 (MQ=255)
                                                                                                                             ccTCGCTGTAGGCGCCCAGCGGGCCGCCGGCATGTGCCACACGGTCTCCTACCTTCAGCCGGGACACGCCTTCGCCCACGGCCTCTACCACACCTGCCGCTTCGGTACCCAGGCCAGAAGGCAGGGACGGCGGCGCATACAGCCCGCtt                   >  4:143355/1‑149 (MQ=255)
                                                                                                                              cTCGCTGTAGGCGCCCAGCGGGCCGCCGGCATGTGCCACACGGTCTCCTACCTTCAGCCGGGACACGCCTTCGCCCACGGCCTCTACCACACCTGCCGCTTCGGTACCCAGGCCAGAAGGCAGGGACGGCGGCGCATACAGCCCGCTTc                  >  1:19943/1‑149 (MQ=255)
                                                                                                                              cTCGCTGTAGGCGCCCAGCGGGCCGCCGGCATGTGCCACACGGTCTCCTACCTTCAGCCGGGACACGCCTTCGCCCACGGCCTCTACCACACCTGCCGCTTCGGTACCCAGGCCAGAAGGCAGGGACGGCGGCGCATACAGCCCGCTTc                  >  6:274945/1‑149 (MQ=255)
                                                                                                                              cTCGCTGTAGGCGCCCAGCGGGCCGCCGGCATGTGCCACACGGTCTCCTACCTTCAGCCGGGACACGCCTTCGCCCACGGCCTCTACCACACCTGCCGCTTCGGTACCCAGGCCAGAAGGCAGGGACGGCGGCGCATACAGCCCGCTTc                  >  3:354740/1‑149 (MQ=255)
                                                                                                                                 gCTGTAGGCGCCCAGCGGGCCGCCGGCATGTGCCACACGGTCTCCTACCTTCAGCCGGGACACGCCTTCGCCCACGGCCTCTACCACACCTGCCGCTTCGGTACCCAGGCCAGAAGGCAGGGACGGCGGCGCATACAGCCCGCTTCgg                <  7:78881/148‑1 (MQ=255)
                                                                                                                                 gCTGTAGGCGCCCAGCGGGCCGCCGGCATGTGCCACACGGTCTCCTACCTTCAGCCGGGACACGCCTTCGCCCACGGCCTCTACCACACCTGCCGCTTCGGTACCCAGGCCAGAAGGCAGGGACGGCGGCGCATACAGCCCGCTTCgg                <  2:19943/148‑1 (MQ=255)
                                                                                                                                  cTGTAGGCGCCCAGCGGGCCGCCGGCATGTGCCACACGGTCTCCTACCTTCAGCCGGGACACGCCTTCGCCCACGGCCTCTACCACACCTGCCGCTTCGGTACCCAGGCCAGAAGGCAGGGACGGCGGCGCATACAGCCCGCTTCGGaa              <  8:289183/149‑1 (MQ=255)
                                                                                                                                   tGTAGGCGCCCAGCGGGCCGCCGGCATGTGCCACACGGTCTCCTACCTTCAGCCGGGACACGCCTTCGCCCACGGCCTCTACCACACCTGCCGCTTCGGTACCCAGGCCAGAAGGCAGGGACGGCGGCGCATACAGCCCGCTTCGGAAg             >  5:106013/1‑149 (MQ=255)
                                                                                                                                     tAGGCGCCCAGCGGGCCGCCGGCATGTGCCACACGGTCTCCTACCTTCAGCCGGGACACGCCTTCGCCCACGGCCTCTACCACACCTGCCGCTTCGGTACCCAGGCCAGAAGGCAGGGACGGCGGCGCATACAGCCCGCTTCGGAAGt            <  3:86981/148‑1 (MQ=255)
                                                                                                                                      cggcggcCCGCTGGGCGCCGGCATGTGCCACACGGTCTCCTACCTTCAGCCGGGACACGCCTTCGCCCACGGCCTCTACCACACCTGCCGCTTCGGTACCCAGGCCAGAAGGCAGGGACGGCGGCGCATACAg                          >  2:41530/2‑133 (MQ=255)
                                                                                                                                      cggcggcCCGCTGGGCGCCGGCATGTGCCACACGGTCTCCTACCTTCAGCCGGGACACGCCTTCGCCCACGGCCTCTACCACACCTGCCGCTTCGGTACCCAGGCCAGAAGGCAGGGACGGCGGCGCATACAg                          <  1:41530/132‑1 (MQ=255)
                                                                                                                                        gcgcCCAGCGGGCCGCCGGCATGTGCCACACGGTCTCCTACCTTCAGCCGGGACACGCCTTCGCCCACGGCCTCTACCACACCTGCCGCTTCGGTACCCAGGCCAGAAGGCAGGGACGGCGGCGCATACAGCCCGCTTCGGAAGTAgg         >  4:282926/1‑148 (MQ=255)
                                                                                                                                              aGCGGGCCGCCGGCATGTGCCACACGGTCTCCTACCTTCAGCCGGGACACGCCTTCGCCCACGGCCTCTACCACACCTGCCGCTTCGGTACCCAGGCCAGAAGGCAGGGACGGCGTCGCATACAGCCCGCTTCGGAAGTAGGTATCGAt  <  4:354740/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
CTTCACCTCATAGGTCTGCTTCAGCAGGTACTGTACGGTCAAACCCTTGAGCATCACCGCTGCCGCCTGCTCGAAGCTGATGGTGTCCGGCAGCTTGACCAGGTTGGCTTCCGGTAGCGTATGCACCTCGCTGTAGGCGCCCAGCGGGCCGCCGGCATGTGCCACACGGTCTCCTACCTTCAGCCGGGACACGCCTTCGCCCACGGCCTCTACCACACCTGCCGCTTCGGTACCCAGGCCAGAAGGCAGGGACGGCGGCGCATACAGCCCGCTTCGGAAGTAGGTATCGAT  >  NC_002947/82712‑83002

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: