Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,594,068 A→C 7.0% intergenic (‑266/+35) PP_1399 ← / ← kgtP hypothetical protein/alpha‑ketoglutarate permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,594,0680AC 7.0% 218.7 / 6.4 71intergenic (‑266/+35)PP_1399/kgtPhypothetical protein/alpha‑ketoglutarate permease
Reads supporting (aligned to +/- strand):  ref base A (32/33);  new base C (3/2);  total (35/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

GCGTATGTCGTGCTTCAGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCG  >  NC_002947/1593924‑1594214
                                                                                                                                                |                                                                                                                                                  
gcgTATGTCGTGCTTCAGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATgg                                                                                                                                                >  8:27266/1‑149 (MQ=255)
gcgTATGTCGTGCTTCAGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATgg                                                                                                                                                <  5:301040/149‑1 (MQ=255)
      gTCGTGCTTCAGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGcc                                                                                                                                          >  6:24304/1‑149 (MQ=255)
           gCTTCAGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCttt                                                                                                                                     >  5:31474/1‑149 (MQ=255)
            cTTCAGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCttt                                                                                                                                     <  4:227493/148‑1 (MQ=255)
                  gCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCGCCATGGAAGGGGGCGCCCCTGTTAGTACt                                                                                                                              >  4:271121/1‑149 (MQ=255)
                   ccGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCtt                                                                                                                                      >  1:124511/1‑140 (MQ=255)
                   ccGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCtt                                                                                                                                      <  2:124511/140‑1 (MQ=255)
                            tggctTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTgccaggccag                                                                                                                    <  4:212245/149‑1 (MQ=255)
                                  gggAGCGTTCGCGCGGGCGCCATCGTTCTGCCAGCCATTGCAGCTTTATAAGGTGGAATTAATCGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTa                                                                                                               <  3:271121/148‑1 (MQ=255)
                                   ggAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACt                                                                                                             >  8:58749/1‑149 (MQ=255)
                                   agaacGTTGGCGCCCGCGCGAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACt                                                                                                             >  2:13917/5‑149 (MQ=255)
                                     aGCGTTCGCGCGGGCGACAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCa                                                                                                           <  5:307703/149‑1 (MQ=255)
                                        gTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTg                                                                                                        <  6:31474/149‑1 (MQ=255)
                                            gcgcgGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGgcgc                                                                                                                                           <  3:197577/110‑1 (MQ=255)
                                            gcgcgGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGgcgc                                                                                                                                           >  4:197577/1‑110 (MQ=255)
                                            gcgcgGGCGCCAACGTTATGCCAGCCAAACCAGCTTCATAACGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGCGGCGCTCCATACAATCGGGCGCCCCTTTTAGTACTGCCACGCCAGCGTTACTAAGTTGTa                                                                                                     >  4:230342/1‑148 (MQ=255)
                                                 ggCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTa                                                                                                               <  7:195059/133‑1 (MQ=255)
                                                 ggCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTa                                                                                                               >  8:195059/1‑133 (MQ=255)
                                                      cAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTa                                                                                                                                   <  6:169143/108‑1 (MQ=255)
                                                      cAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTa                                                                                                                                   >  5:169143/1‑108 (MQ=255)
                                                      cAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGtg                                                                                          <  5:235958/149‑1 (MQ=255)
                                                         cGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGTAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGtgctg                                                                                       <  8:101060/149‑1 (MQ=255)
                                                          gTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGtgctgc                                                                                      >  6:224614/1‑149 (MQ=255)
                                                            tCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTa                                                                                                               >  6:1352/1‑122 (MQ=255)
                                                            tCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTa                                                                                                               <  5:1352/122‑1 (MQ=255)
                                                               gccagccaTTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCCCCATTGCATGGAGCGcc                                                                                                                                          >  3:283272/1‑92 (MQ=38)
                                                               gccagccaTTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCCCCATTGCATGGAGCGcc                                                                                                                                          <  4:283272/92‑1 (MQ=38)
                                                                ccagccaTTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTgg                                                                                  <  1:13917/147‑1 (MQ=255)
                                                                   gccaTTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCt                                                                             <  7:150509/149‑1 (MQ=255)
                                                                       ttGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCa                                                                         <  2:292584/149‑1 (MQ=255)
                                                                       ttGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCa                                                                         >  8:23903/1‑149 (MQ=255)
                                                                       ttGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCa                                                                         <  5:24304/149‑1 (MQ=255)
                                                                       ttGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACAGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCa                                                                         >  1:292584/1‑149 (MQ=255)
                                                                          cAGCTTTATAAGGTGGAATTAATAGGGCGTTAAGTAGGAAAAGTCTGAATAGAAAAGGGGCGCTCGATGCTATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATg                                                                       <  3:230342/148‑1 (MQ=255)
                                                                          cAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATgg                                                                      >  5:224671/1‑149 (MQ=255)
                                                                          cAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATgg                                                                      >  5:229720/1‑149 (MQ=255)
                                                                          cAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATg                                                                       <  7:27266/148‑1 (MQ=255)
                                                                                 aTAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGt                                                                <  7:23903/148‑1 (MQ=255)
                                                                                  tAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTa                                                               >  6:335776/1‑148 (MQ=255)
                                                                                    aGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAAc                                                            <  2:319011/149‑1 (MQ=255)
                                                                                       tGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTAGCGTAATCcag                                                         >  5:224441/1‑149 (MQ=255)
                                                                                         gAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACcag                                                         >  6:293784/1‑147 (MQ=255)
                                                                                            ttAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCagcagc                                                     >  8:175884/1‑148 (MQ=255)
                                                                                               aTAGGGCGTTAAGAAGGAAAAGTCTGAATAGAGAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAg                                                  <  5:213382/148‑1 (MQ=255)
                                                                                                 aGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAg                                               >  3:32184/1‑149 (MQ=255)
                                                                                                  gggCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGa                                                                                             <  4:15626/102‑1 (MQ=255)
                                                                                                  gggCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGa                                                                                             >  3:15626/1‑102 (MQ=255)
                                                                                                   ggCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGc                                             <  5:224614/149‑1 (MQ=255)
                                                                                                     cGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGa                                           <  7:175884/149‑1 (MQ=255)
                                                                                                               gAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGt                                 <  1:300644/149‑1 (MQ=255)
                                                                                                                 aaaGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGa                               <  6:229720/149‑1 (MQ=255)
                                                                                                                 aaaGGCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGa                               <  6:224671/149‑1 (MQ=255)
                                                                                                                     tCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGt                                 >  3:309492/1‑143 (MQ=255)
                                                                                                                     tCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGt                                 <  4:309492/143‑1 (MQ=255)
                                                                                                                     tCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTa                           >  3:249505/1‑149 (MQ=255)
                                                                                                                         aaTAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAg                       <  5:335776/149‑1 (MQ=255)
                                                                                                                           tAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAg                       >  6:283638/1‑147 (MQ=255)
                                                                                                                              aaaaGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCttt                                                                           <  1:354234/92‑1 (MQ=255)
                                                                                                                              aaaaGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCttt                                                                           >  2:354234/1‑92 (MQ=255)
                                                                                                                              aaaaGGGGCGCCCCATTGCATGGAGCGCACCTTTTAGTACTGCCAGGACAGCGTTACTAAGTGGTAATCCGAGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAg                    >  1:171552/1‑147 (MQ=255)
                                                                                                                               aaaGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAg                 <  7:58749/149‑1 (MQ=255)
                                                                                                                                   gggCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGc             >  6:108468/1‑149 (MQ=255)
                                                                                                                                   gggCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGATACTCAGTGGTAATCCGCGAGTGCAGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGc             >  2:351431/1‑149 (MQ=255)
                                                                                                                                    ggCGCCCCATTGCATGGAGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGc                >  7:195100/1‑145 (MQ=255)
                                                                                                                                    ggCGCCCCATTGCATGGAGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGc                <  8:195100/145‑1 (MQ=255)
                                                                                                                                     gcgcTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGc                                             >  3:33354/1‑115 (MQ=255)
                                                                                                                                     gcgcTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGc                                             <  4:33354/115‑1 (MQ=255)
                                                                                                                                      cgcTCCATGCAATGGGGCTCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGctct          <  5:283638/149‑1 (MQ=255)
                                                                                                                                            aTGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATcc    >  2:350719/1‑149 (MQ=255)
                                                                                                                                              gCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCg  >  7:224114/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
GCGTATGTCGTGCTTCAGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCG  >  NC_002947/1593924‑1594214

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: