Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,445,250 C→G 9.7% G117A (GGG→GCG) ‡ nicE ← maleate isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,445,2500CG 9.7% 178.0 / 10.2 62G117A (GGG→GCG) ‡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 = 4.46e-02

GGGGATTTCCAGCGCGCGCCAGTCCTGCACCTCGAAGCCTTCCTCGCGGATGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCATGATCGCCACAAGGCAGGCATAGCCCAGCACATCGACCTTGGCGTCGGA  >  NC_002947/4445105‑4445395
                                                                                                                                                 |                                                                                                                                                 
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                                                                                                            <  4:164515/140‑1 (MQ=255)
                                             gcgGATGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTcgg                                                                                                            >  3:164515/1‑140 (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                                                                                    >  2:283127/1‑149 (MQ=255)
                                                            caccaACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCtttt                                                                                    >  4:43795/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                                                   >  6:136031/1‑149 (MQ=255)
                                                                                             gacgaTGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCa                                                   >  3:235148/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                      >  6:296619/1‑145 (MQ=255)
                                                                                                                              cAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCATGATCGCCACAAGGCAGGCATAGCCCAGc                      <  5:296619/145‑1 (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)
                                                                                                                                                 |                                                                                                                                                 
GGGGATTTCCAGCGCGCGCCAGTCCTGCACCTCGAAGCCTTCCTCGCGGATGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGGTGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCATGATCGCCACAAGGCAGGCATAGCCCAGCACATCGACCTTGGCGTCGGA  >  NC_002947/4445105‑4445395

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