Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,445,251 C→G 9.4% G117R (GGG→CGG) ‡ nicE ← maleate isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,445,2510CG 9.4% 176.9 / 7.8 62G117R (GGG→CGG) ‡nicEmaleate isomerase
Reads supporting (aligned to +/- strand):  ref base C (31/25);  new base G (3/3);  total (34/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.69e-01

GGGGATTTCCAGCGCGCGCCAGTCCTGCACCTCGAAGCCTTCCTCGCGGATGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCATGATCGCCACAAGGCAGGCATAGCCCAGCACATCGACCTTGGCGTCGGACAGC  >  NC_002947/4445105‑4445399
                                                                                                                                                  |                                                                                                                                                    
ggggATTTCCAGCGCGCCCCAGTCCTCCACCTCGAAGCCTTCCTCGCGGATGTAGTTCACCACCAACTCCGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCCCATGCAAGGCCTCGACCAGCGCCCCgg                                                                                                                                                    >  2:125292/1‑149 (MQ=255)
   gATTTCCAGCGCGCGCCAGTCCTGCACCTCGAAGCCTTCCTCGCGGATGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCGGGgcgc                                                                                                                                                 <  3:250605/149‑1 (MQ=255)
     tttCCAGCGCGCGCCAGTCCTGCACCTCGAAGCCTTCCTCGCGGATGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTg                                                                                                                                               <  3:141983/149‑1 (MQ=255)
     tttCCAGCGCGCGCCAGTCCTGCACCTCGAAGCCTTCCTCGCGGATGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTg                                                                                                                                               <  7:169157/149‑1 (MQ=255)
       tCCAGCGCGCGCCAGTCCTGCACCTCGAAGCCTTCCTCGCGGATGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCGGGGCGCTgg                                                                                                                                              <  5:125375/148‑1 (MQ=255)
        ccAGCGCGCGCCAGTCCTGCACCTCGAAGCCTTCCTCGCGGATGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGt                                                                                                                                             <  6:255745/148‑1 (MQ=255)
         cAGCGCGCGCCAGTCCTGCACCTCGAAGCCTTCCTCGCGGATGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCGGGGCGCTGGtcg                                                                                                                                           <  1:255054/149‑3 (MQ=255)
             gcgcgcCAGTCCTGCACCTCGAAGCCTTCCTCGCGGATGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCAc                                                                                                                                       >  2:165746/1‑149 (MQ=255)
               gcgcCAGTCCTGCACCTCGAAGCCTTCCTCGCGGATGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCg                                                                                                                                     <  1:115915/149‑1 (MQ=255)
               gcgcCAGTCCTGCACCTCGAAGCCTTCCTCGCGGATGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCg                                                                                                                                     <  6:269463/149‑1 (MQ=255)
                   cAGTCCTGCACCTCGAAGCCTTCCTCGCGGATGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGc                                                                                                                                 >  8:283912/1‑149 (MQ=255)
                      tCCTGCACCTCGAAGCCTTCCTCGCGGATGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAg                                                                                                                              >  2:192405/1‑149 (MQ=255)
                      tCCTGCACCTCGAAGCCTTCCTCGCGGATGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAg                                                                                                                              <  1:192405/149‑1 (MQ=255)
                          gCACCTCGAAGCCTTCCTCGCGGATGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGca                                                                                                                          >  8:17995/1‑149 (MQ=255)
                          gCACCTCGAAGCCTTCCTCGCGGATGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGc                                                                                                                           <  6:87635/148‑1 (MQ=255)
                          gCACCTCGAAGCCTTCCTCGCGGATGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCCCTGGTGATCACCGGTGCCAGTGc                                                                                                                           >  5:87635/1‑148 (MQ=255)
                                cGAAGCCTTCCTCGCGGATGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTg                                                                                                                    >  7:181604/1‑149 (MQ=255)
                                             gcgGATGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTcgg                                                                                                                >  3:164515/1‑140 (MQ=255)
                                             gcgGATGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTcgg                                                                                                                <  4:164515/140‑1 (MQ=255)
                                                 aTGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGGTGGCctgctg                                                                                                   >  2:136503/1‑149 (MQ=255)
                                                  tGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGGTGGCctgctg                                                                                                   >  5:241266/1‑148 (MQ=255)
                                                   gTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGGTGGCctgctg                                                                                                   >  4:6434/1‑147 (MQ=255)
                                                    tAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAg                                                                                                >  5:343565/1‑149 (MQ=255)
                                                           ccaccaACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCttt                                                                                         >  8:258352/1‑149 (MQ=255)
                                                            caccaACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGgcgc                                                                                            <  8:224375/145‑1 (MQ=255)
                                                            caccaACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGgcgc                                                                                            >  7:224375/1‑145 (MQ=255)
                                                            caccaACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCtttt                                                                                        >  4:43795/1‑149 (MQ=255)
                                                            caccaACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCtttt                                                                                        >  2:283127/1‑149 (MQ=255)
                                                             accaACTCGGTCAGCGGCTTCATGTAAGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCtttt                                                                                        >  1:203827/1‑148 (MQ=255)
                                                                     ggTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGATTgg                                                                                <  1:136503/148‑1 (MQ=255)
                                                                       tCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGTCGGTGGCCTGCTGCAGGCGCTTTTCAGATTGGCgg                                                                             <  4:161198/149‑1 (MQ=255)
                                                                       tCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGATTGGCgg                                                                             <  7:17995/149‑1 (MQ=255)
                                                                           cGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGa                                                                          >  8:349596/1‑148 (MQ=255)
                                                                                   tGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATagcccagc                                                                 >  7:230464/1‑149 (MQ=255)
                                                                                             gacgaTGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCa                                                       >  3:235148/1‑149 (MQ=255)
                                                                                             gacgaTGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCa                                                       >  6:136031/1‑149 (MQ=255)
                                                                                                     cAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCATGATCGcc                                               >  8:129188/1‑149 (MQ=255)
                                                                                                       aTGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCATGATCGCcac                                             <  6:241266/149‑1 (MQ=255)
                                                                                                          cgcTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCATGATCGCCACaa                                           <  8:181604/148‑1 (MQ=255)
                                                                                                            cTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCATGATCGCCACAagg                                         >  1:368626/1‑148 (MQ=255)
                                                                                                                gCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCTGTGCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCATGATCGCCACAaggcagg                                     <  7:217751/148‑1 (MQ=255)
                                                                                                                gCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCATGATCGCCACAaggcaggc                                    <  4:235148/149‑1 (MQ=255)
                                                                                                                   ttCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTCGTGATCACCGGTGCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTAAGATTGGCGGTGATAGCCCAGCCCCATGGCCATGATCGCCACAAggcaggca                                   >  4:330004/1‑147 (MQ=255)
                                                                                                                      atcaCATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCATGATCGCCACAAGGCAGGCATAGc                               <  7:258352/148‑1 (MQ=255)
                                                                                                                       tcaCATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCATGATCGCCACAAGGCAGGCATAGccc                             <  8:230464/149‑1 (MQ=255)
                                                                                                                          caTGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCATGATCGCCACAAGGCAGGCATAGCCCAGc                          <  3:6434/149‑1 (MQ=255)
                                                                                                                           aTGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCATGATCGCCACAAGGCAGGCATAGCCCAGca                         >  5:304233/1‑149 (MQ=255)
                                                                                                                             gCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGGTGGACTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCATGATCGCCACAAGGCAGGCATAGCCCAGcaca                       >  3:47647/1‑149 (MQ=255)
                                                                                                                             gCAAGGCCTCGACCAGCGCCCCGGCGCTGATGATCACCGGTGCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCATGATCGCCACAAGGCAGGCATAGcc                              <  6:204021/142‑1 (MQ=255)
                                                                                                                             gCAAGGCCTCGACCAGCGCCCCGGCGCTGATGATCACCGGTGCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCATGATCGCCACAAGGCAGGCATAGcc                              >  5:204021/1‑142 (MQ=255)
                                                                                                                              cAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCATGATCGCCACAAGGCAGGCATAGCCCAGc                          <  5:296619/145‑1 (MQ=255)
                                                                                                                              cAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCATGATCGCCACAAGGCAGGCATAGCCCAGc                          >  6:296619/1‑145 (MQ=255)
                                                                                                                                     tCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCATGATCGCCACAAGGCAGGCATAGCCCAGCACATCGACCtt               >  5:286957/1‑149 (MQ=255)
                                                                                                                                      tcaCCAGCGCCGGGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCATGATCGCCACAAGGCAGGCATAGCCCAGCACATCGACCTTg              >  8:75368/3‑149 (MQ=255)
                                                                                                                                       gACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCATGATCGCTACAAGGCAGGCATAGCCCa                            >  7:262159/1‑134 (MQ=255)
                                                                                                                                       gACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCATGATCGCTACAAGGCAGGCATAGCCCa                            <  8:262159/134‑1 (MQ=255)
                                                                                                                                        aCCAGCGCCGGGGCGCTGGTGATCACCGGTGCCAGTGCATCACTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCATGATCGCCACAAGGCAGGCATAGCCCAGCACATCGACCTTGGc            >  4:71611/1‑149 (MQ=255)
                                                                                                                                         ccAGCGCCGGGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCATGATCGCCACAAGGCAGGCATAGCCCAGCACATCGACCTTGGCg           >  8:194264/1‑149 (MQ=255)
                                                                                                                                         ccAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCATGATCGCCACAAGGCAGGCATAGCCCAGCACATCGACCTTGGCg           >  2:360790/1‑149 (MQ=255)
                                                                                                                                            gcgcCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCATGATCGCCACTAGGCAGGCATAGCCCAGCACATCGACCTTGGCGTCg        <  6:343565/149‑1 (MQ=255)
                                                                                                                                            gcgcCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCATGATCGCCACAAGGCAGGCATAGCCCAGCACATCGACCTTGGCGTCg        <  2:368626/149‑1 (MQ=255)
                                                                                                                                              gcCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCATGATCGCCACAAGGCAGGCATAGCCCAGCACATCGACCTTGGCGTCggt      <  7:129188/149‑2 (MQ=255)
                                                                                                                                                cccGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCATGATCGCCACAAGGCAGGCATAGCCCAGCACATCGACCTTGGCGTc         <  2:145418/144‑1 (MQ=255)
                                                                                                                                                cccGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCATGATCGCCACAAGGCAGGCATAGCCCAGCACATCGACCTTGGCGTc         >  1:145418/1‑144 (MQ=255)
                                                                                                                                                cccGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCATGATCGCCACAAGGCAGGAATAGCCCAGCACATCGACCTTGGCGTCGGa      >  5:24292/1‑147 (MQ=255)
                                                                                                                                                  cGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCATGATCGCCACAAGGCAGGCATAGCCCAGCACATCGACCTTGGCGTCGGACAGc  <  7:349596/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GGGGATTTCCAGCGCGCGCCAGTCCTGCACCTCGAAGCCTTCCTCGCGGATGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCATGATCGCCACAAGGCAGGCATAGCCCAGCACATCGACCTTGGCGTCGGACAGC  >  NC_002947/4445105‑4445399

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