Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,594,067 C→G 13.7% intergenic (‑265/+36) 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,0670CG13.7% 196.7 / 21.6 73intergenic (‑265/+36)PP_1399/kgtPhypothetical protein/alpha‑ketoglutarate permease
Reads supporting (aligned to +/- strand):  ref base C (30/33);  new base G (7/3);  total (37/36)
Fisher's exact test for biased strand distribution p-value = 3.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

TGTCGTGCTTCAGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTC  >  NC_002947/1593929‑1594207
                                                                                                                                          |                                                                                                                                              
tGTCGTGCTTCAGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGgcgc                                                                                                                                      <  8:321013/149‑1 (MQ=255)
  tCGTGCTTCAGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCACCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGcc‑‑c                                                                                                                                  >  2:42464/1‑149 (MQ=255)
    gTGCTTCAGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGcc‑‑cc                                                                                                                                 <  7:31007/148‑1 (MQ=255)
        ttCAGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCttt                                                                                                                              >  5:188293/1‑147 (MQ=255)
        ttCAGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTa                                                                                                                            >  7:159985/1‑149 (MQ=255)
        ttCAGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTa                                                                                                                            >  7:53658/1‑149 (MQ=255)
           aGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTa                                                                                                                         >  8:226930/1‑149 (MQ=255)
           aGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCCCCATTGCATGGAgcgc                                                                                                                                      <  6:16229/138‑1 (MQ=255)
           aGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCCCCATTGCATGGAgcgc                                                                                                                                      >  5:16229/1‑138 (MQ=255)
            ggCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTAc                                                                                                                        >  6:251775/1‑149 (MQ=255)
             gCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACt                                                                                                                       <  5:206344/149‑1 (MQ=255)
                gCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTgcc                                                                                                                    <  8:159985/149‑1 (MQ=255)
                gCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTgc                                                                                                                     >  2:73971/1‑148 (MQ=255)
                  aaGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGATTTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTgccag                                                                                                                  >  3:335385/1‑149 (MQ=255)
                      ctggctTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTgccaggcca                                                                                                              <  5:251775/149‑1 (MQ=255)
                       tggctTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTgccaggccag                                                                                                             >  2:291956/1‑149 (MQ=255)
                        ggctTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGc                                                                                                            <  3:332785/149‑1 (MQ=255)
                          ctTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGt                                                                                                          >  5:206049/1‑149 (MQ=255)
                              ggAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAATAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTAc                                                                                                       <  6:247535/148‑1 (MQ=255)
                              agaacGTTGGCGCCCGCGCGAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACt                                                                                                      >  7:210316/5‑149 (MQ=255)
                               gAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACt                                                                                                      >  7:243556/1‑148 (MQ=255)
                                aGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCa                                                                                                    <  5:200381/149‑1 (MQ=255)
                                aGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCa                                                                                                    <  1:42464/149‑1 (MQ=255)
                                   gTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTg                                                                                                 <  1:73971/149‑1 (MQ=255)
                                    ttCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTgg                                                                                                <  8:53658/149‑1 (MQ=255)
                                    ttCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTgg                                                                                                >  8:152285/1‑149 (MQ=255)
                                    ttCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTgg                                                                                                >  8:157274/1‑149 (MQ=255)
                                                ccAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGt                                                                                    <  7:157274/149‑1 (MQ=255)
                                                ccAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGt                                                                                    <  7:152285/149‑1 (MQ=255)
                                                 cAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCa                                                                                                    >  6:43908/1‑132 (MQ=255)
                                                 cAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCa                                                                                                    <  5:43908/132‑1 (MQ=255)
                                                 cAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTg                                                                                                 >  7:344609/1‑135 (MQ=255)
                                                 cAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTg                                                                                                 <  8:344609/135‑1 (MQ=255)
                                                 cAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGtg                                                                                   >  2:49019/1‑149 (MQ=255)
                                                  aaCGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCt                                                                                                                                >  2:60735/1‑103 (MQ=255)
                                                  aaCGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCt                                                                                                                                <  1:60735/103‑1 (MQ=255)
                                                    cGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGtgctg                                                                                <  7:226930/149‑1 (MQ=255)
                                                          gccagccaTTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGgt                                                                          >  6:294583/1‑149 (MQ=255)
                                                           ccagccaTTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGgtg                                                                         <  4:335385/149‑1 (MQ=255)
                                                           ccagccaTTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGgt                                                                          <  8:249542/148‑1 (MQ=255)
                                                           ccagccaTTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGgt                                                                          <  5:294583/148‑1 (MQ=255)
                                                               ccaTTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCTCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTc                                                                       <  1:291956/149‑1 (MQ=255)
                                                                     cAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATg                                                                <  6:206049/148‑1 (MQ=255)
                                                                           gATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGtg                                                                                   <  8:167101/123‑1 (MQ=255)
                                                                           gATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGtg                                                                                   >  7:167101/1‑123 (MQ=255)
                                                                             tAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTa                                                        >  3:140078/1‑148 (MQ=255)
                                                                                 gTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAAc                                                     >  3:264308/1‑147 (MQ=255)
                                                                                    gAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACcagca                                                <  8:243556/149‑1 (MQ=255)
                                                                                          aTAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTg                                                                            <  4:315383/115‑1 (MQ=255)
                                                                                                cGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGa                                    <  6:188293/149‑1 (MQ=255)
                                                                                                cGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGa                                    >  7:13996/1‑149 (MQ=255)
                                                                                                  ttAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATa                                  <  3:310983/149‑1 (MQ=255)
                                                                                                       aagGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAg                                           >  1:114520/1‑135 (MQ=255)
                                                                                                       aagGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAg                                           <  2:114520/135‑1 (MQ=255)
                                                                                                         ggAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCgg                           <  4:140078/149‑1 (MQ=255)
                                                                                                           aaaaGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGt                          <  4:109668/148‑1 (MQ=255)
                                                                                                            cctctCTGAATCGAACAGGGCCGCTCCATGCAAGGGGGCGCC‑‑CCTTTTAGTACTGCCACGCCAGCGCTACTCCGACGCAATCCGCGAGTGCCGCTTGGTGTCTTTCCTGGTCGCGTAAACCAGCAGCGAGCAGGCGACACAGGCGGTGa                        <  8:318957/145‑1 (MQ=255)
                                                                                                             aaGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGAc                       <  1:49019/149‑1 (MQ=255)
                                                                                                                  tGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGa                        >  7:318957/1‑143 (MQ=255)
                                                                                                                  tGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTAc                   >  6:325012/1‑148 (MQ=255)
                                                                                                                    aaTAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCGGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTGCCAg                >  1:7218/1‑149 (MQ=255)
                                                                                                                    aaTAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTAc                   <  5:325012/146‑1 (MQ=255)
                                                                                                                       ctaaaaGGGGCGCCCCATTGCATGGAGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTAc                   <  6:180545/141‑1 (MQ=39)
                                                                                                                       ctaaaaGGGGCGCCCCATTGCATGGAGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTAc                   >  5:180545/3‑143 (MQ=39)
                                                                                                                         aaaaGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAg             >  2:201115/1‑147 (MQ=255)
                                                                                                                          aaaGGGGCGCCCCATTGCATGGAGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAg          >  4:171885/1‑149 (MQ=255)
                                                                                                                            aGGGGCGCCCCATTGCATGGAGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGcc        >  4:24378/1‑149 (MQ=255)
                                                                                                                             ggggCGCCCCATTGCATGGAGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCg       <  3:24378/149‑1 (MQ=255)
                                                                                                                               ggCGCCCCATTGCATGGAGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGct     >  7:325065/1‑149 (MQ=255)
                                                                                                                               ggCGCCCCATTGCATGGAGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGct     >  3:42649/1‑149 (MQ=255)
                                                                                                                               ggCGCCCCATTGCATGGAGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGct     >  5:18034/1‑149 (MQ=255)
                                                                                                                                gcgcTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCATCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGctc    >  7:142030/1‑149 (MQ=255)
                                                                                                                                   cTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGTCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCtctc  >  4:361339/1‑148 (MQ=255)
                                                                                                                                          |                                                                                                                                              
TGTCGTGCTTCAGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCC‑‑CCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTC  >  NC_002947/1593929‑1594207

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: