Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,244,285 C→A 5.1% R66L (CGC→CTC)  hmgA ← homogentisate 1,2‑dioxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,244,2850CA 5.1% 257.9 / 2.2 78R66L (CGC→CTC) hmgAhomogentisate 1,2‑dioxygenase
Reads supporting (aligned to +/- strand):  ref base C (36/38);  new base A (2/2);  total (38/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.62e-01

GCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACAGCAGCTCGGCATACAGGCCATACGGGGCCTTCTGCGGGGAGTTCTGCCCAACAGGCAGCGCCCCGGGCAATGCTTCGCTGGCGAATTCGTT  >  NC_002947/5244139‑5244433
                                                                                                                                                  |                                                                                                                                                    
gCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCgg                                                                                                                                                    <  1:239400/149‑1 (MQ=255)
 ccAACCTTCGATGAAATCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGt                                                                                                                                                   <  5:179017/149‑1 (MQ=255)
   aaCCTTCGATGAAATCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTa                                                                                                                                                  >  3:86526/1‑148 (MQ=255)
       ttCGATGAAATCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCGCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGcc                                                                                                                                             <  6:243998/149‑1 (MQ=255)
         cGATGAAATCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGcc                                                                                                                                             <  7:57544/147‑1 (MQ=255)
          gATTAAATGGGTCTGTGCAGCCGGAATCGGGTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAg                                                                                                                                           <  7:45228/148‑1 (MQ=255)
          gATGAAATCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAgg                                                                                                                                          >  3:110833/1‑149 (MQ=255)
          gATGAAATCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAg                                                                                                                                           <  4:293329/148‑1 (MQ=255)
              aactgggtcggtTCAGCAGGAATCGGCTGCGGGCGCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACg                                                                                                                                      <  7:148419/144‑1 (MQ=255)
              aaaTCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTAc                                                                                                                                       <  2:122434/148‑1 (MQ=255)
                 tcggtcggTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGgcg                                                                                                                                   >  7:178613/1‑149 (MQ=255)
                    gtcggtTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATACCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAg                                                                                                                                >  1:1368/1‑149 (MQ=255)
                          tCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGTCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCt                                                                                                                          >  6:199330/1‑149 (MQ=255)
                             gcagGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGAGGTACAGCCAGGTACGGCGCAGCTCGCTGCg                                                                                                                       >  4:189828/1‑149 (MQ=255)
                              cagGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGGTGCTGCGACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGCAGCGTGGGTGCAAGGCGGCAGGGCGAATGCTGTTCAGCCTGGTACGGCGCAGAGCGCTGCgg                                                                                                                      >  4:26411/1‑149 (MQ=255)
                                ggAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGc                                                                                                                    <  5:200188/149‑1 (MQ=255)
                                 gAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGc                                                                                                                    >  3:107174/1‑148 (MQ=255)
                                  aaTCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGTGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCa                                                                                                                  >  3:227748/1‑149 (MQ=255)
                                    tCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATg                                                                                                                >  1:145770/1‑149 (MQ=255)
                                      ggCTGCGGGCTGCAGCGCAGGCGGTTGGGGTTGATGCCAGCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGGGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGt                                                                                                              <  6:157515/149‑1 (MQ=255)
                                      ggCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGt                                                                                                              >  7:194193/1‑149 (MQ=255)
                                       gCTGCGGGCTCCCGCGCAGGCGGTTGGGGTGGAGGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGGGCGCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCTGGTCGCTGCGGGCCATGGTg                                                                                                             <  7:254225/149‑1 (MQ=255)
                                        cTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGa                                                                                                            <  5:228795/149‑1 (MQ=255)
                                            gggCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGc                                                                                                        >  2:224215/1‑149 (MQ=255)
                                            gggCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGc                                                                                                        >  4:211471/1‑149 (MQ=255)
                                               cTCCAGCGCAGTCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCTAGGCGCTGGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGGGGGCCATGGTGAACTCCGt                                                                                                     <  4:240152/149‑1 (MQ=255)
                                               cTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGt                                                                                                     >  3:94032/1‑149 (MQ=255)
                                                  cAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGc                                                                                                   <  4:110833/148‑1 (MQ=255)
                                                  cAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGAGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGc                                                                                                   <  4:154860/148‑1 (MQ=255)
                                                   aGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGTGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGccc                                                                                                 <  4:227748/149‑1 (MQ=255)
                                                   aGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCACTGCGGGCCATGGTGAACGCCGTGccc                                                                                                 <  8:178613/149‑1 (MQ=255)
                                                   aGCGCAGGCGGTTGGGGTTGATACCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGccc                                                                                                 <  2:1368/149‑1 (MQ=255)
                                                             gTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGcc                                                                                                       <  4:188673/133‑1 (MQ=255)
                                                             gTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGcc                                                                                                       >  3:188673/1‑133 (MQ=255)
                                                                  tgTTGATGCCTCGCTGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCTGTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACTGCGAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACAGCAGCTCGGCa                                                                                  <  4:305520/148‑1 (MQ=255)
                                                                  ggTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACAGCAGCTCGGCa                                                                                  <  2:145770/149‑1 (MQ=255)
                                                                   gTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGt                                                                                                              <  6:190856/120‑1 (MQ=255)
                                                                   gTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGt                                                                                                              >  5:190856/1‑120 (MQ=255)
                                                                   gTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACAGCAGCTCGGCa                                                                                  >  3:165428/1‑148 (MQ=255)
                                                                   gTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACAGCAGCTCGGCa                                                                                  >  6:65866/1‑148 (MQ=255)
                                                                   gTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACAGCAGCTCGGCa                                                                                  >  6:315729/1‑148 (MQ=255)
                                                                   gTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACAGCAGCTCGGCAt                                                                                 <  4:86526/149‑1 (MQ=255)
                                                                    ttGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACAGCAGCTCGGCATa                                                                                >  1:130940/1‑149 (MQ=255)
                                                                    ttGATGACACCCAGTGGCACGCCGAGCGGCTGGCGCGCCAGGCGATCGAAGCGTGGGTGCAAGGCACAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACAGCAGCTCGGCa                                                                                  >  3:305520/1‑147 (MQ=255)
                                                                      gATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACAGCAGCTCGGCATa                                                                                >  6:324756/1‑147 (MQ=255)
                                                                         gCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCg                                                                                                                                                     >  1:59217/1‑75 (MQ=255)
                                                                         gCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCg                                                                                                                                                     <  2:59217/75‑1 (MQ=255)
                                                                          ccACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACAGCAGCTCGGCATACAGGcc                                                                          <  5:19734/149‑1 (MQ=255)
                                                                          ccACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACAGCAGCTCGGCATACAGGc                                                                           >  8:36539/1‑148 (MQ=255)
                                                                            aCCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATgg                                                                                                               >  7:238245/1‑110 (MQ=255)
                                                                            aCCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATgg                                                                                                               <  8:238245/110‑1 (MQ=255)
                                                                              ccAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACAGCAGCTCGGCATACAGGCCATAc                                                                      >  5:34366/1‑149 (MQ=255)
                                                                                 gTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGAGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACAGCAGCTCGGCATACAGGCCATACggg                                                                   <  3:189828/149‑1 (MQ=255)
                                                                                   ggCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACAGCAGCTCGGCATACAGGCCATACGGGGc                                                                 <  4:94032/149‑1 (MQ=255)
                                                                                           gAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACAGCAGCTCGGCATACAGGCCATACGGGGCCTTCTGCg                                                         <  3:211471/149‑1 (MQ=255)
                                                                                             gCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGGACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACAGCAGCTCGGCATACAGGCCATACGGGGCCTTCTGCggg                                                       <  3:81240/149‑1 (MQ=255)
                                                                                                  tGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGa                                                                                               <  2:208204/104‑1 (MQ=255)
                                                                                                  tGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGa                                                                                               >  1:208204/1‑104 (MQ=255)
                                                                                                    gcgcgcCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACAGCAGCTCGGCATACAGGCCATACGGGGCCTTCTGCGGGGAGTTCt                                                <  4:165428/149‑1 (MQ=255)
                                                                                                    gcgcgcCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACAGCAGCTCGGCATACAGGCCATACGGGGCCTTCTGCGGGGAGTTCt                                                <  4:107174/149‑1 (MQ=255)
                                                                                                      gcgcCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACAGCAGCTCGGCATACAGGCCATACGGGGCCTTCTGCGGGGAGTTCTGc                                              <  1:224215/149‑1 (MQ=255)
                                                                                                       cgcCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACAGCAGCTCGGCATACAGGCCATACGGGGCCTTCTGCGGGGAGTTCTGcc                                             <  3:41552/149‑1 (MQ=255)
                                                                                                            ggCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACAGCAGCTCGGCATACAGGCCATACGGGGCCTTCTGCGGGGAGTTCTGCCCaa                                          >  4:173932/1‑147 (MQ=255)
                                                                                                              cgcTCGAAGCCTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACAGCAGCTCGGCATACAGGCCATACGGGGCCTTCTGCGGGGAGTTCTGCCCAACAgg                                      <  6:65393/149‑1 (MQ=255)
                                                                                                                 tCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACAGCAGCTCGGCATACAGGCCATACGGGGCCTTCTGCGGGGAGTTCTGCCCAACAGGCAg                                   >  4:268698/1‑149 (MQ=255)
                                                                                                                         tGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACAGCAGCTCGGCATACAGGCCATACgg                                                                    >  8:107769/1‑108 (MQ=255)
                                                                                                                          gggTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACAGCAGCTCGGCATACAGGCCATACGGGGCCTTCTGCGGGGAGTTCTGCCCAACAGGCAGCGCCCCggg                          >  3:32547/1‑149 (MQ=255)
                                                                                                                           ggTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACAGCAGCTCGGCATACAGGCCATACGGGGCCTTCTGCGGGGAGTTCTGCCCAACAGGCAGCGCCCCggg                          >  6:220987/1‑148 (MQ=255)
                                                                                                                           ggTGCAAGGCAGAAGGGCGAATGAGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACAGCAGCTCgg                                                                                    <  8:49442/90‑1 (MQ=255)
                                                                                                                           ggTGCAAGGCAGAAGGGCGAATGAGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACAGCAGCTCgg                                                                                    >  7:49442/1‑90 (MQ=255)
                                                                                                                              gCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACAGCAGCTCGGCATACAGGCCATACGGGGCCTTCTGCGGGGAGTTCTGCCCAACAGGCAGCGCCCCGGGCAAt                      >  1:162650/1‑149 (MQ=255)
                                                                                                                               cAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACAGCAGCTCGGCATACAGGCCATACGGGGCCTTCTGCGGGGAGTTCTGCCCAACAGGCAGCGCCCCGGGCAAt                      >  5:210585/1‑148 (MQ=255)
                                                                                                                                 aGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACAGCAGCTCGGCATACAGGCCATACGGGGCCTTCTGCGGGGa                                                     <  3:33451/115‑1 (MQ=255)
                                                                                                                                 aGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACAGCAGCTCGGCATACAGGCCATACGGGGCCTTCTGCGGGGa                                                     >  4:33451/1‑115 (MQ=255)
                                                                                                                                    cAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACAGCAGCTCGGCATACAGGCCATACGGGGCCTTCTGCGGGGAGTTCTGCCCAACAGGCAGCGCCCCGGGCAATGCTTCg                >  4:175921/1‑149 (MQ=255)
                                                                                                                                     aGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACAGCAGCTCGGCATACAGGCCATACGGGGCCTTCTGCGGGGAGTTCTGCCCAACAGGCAGCGCCCCGGGCAATGCTTCGc               >  3:36786/1‑149 (MQ=255)
                                                                                                                                        aGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACAGCAGCTCGGCATACAGGCCATACGGGGCCTTCTGCGGGGAGTTCTGCCCAACAGGCAGCGCCCCGGGCAATGCTTCGCTgg            <  3:194366/149‑1 (MQ=255)
                                                                                                                                               aTGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACAGCAGCTCGGCATACAGGCCATACGGGGCCTTCTGCGGGGAGTTCTGCCCAACAGGCAGCGCCCCGGGCAATGCTTCGCTGGCGAATTc     <  5:199330/149‑1 (MQ=255)
                                                                                                                                                  cGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACAGCAGCTCGGCATACAGGCCATACGGGGCCTTCTGCGGGGAGTTCTGCCCAACAGGCAGCGCCCCGGGCAATGCTTCGCTGGCGAATTCGtt  <  6:210585/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACAGCAGCTCGGCATACAGGCCATACGGGGCCTTCTGCGGGGAGTTCTGCCCAACAGGCAGCGCCCCGGGCAATGCTTCGCTGGCGAATTCGTT  >  NC_002947/5244139‑5244433

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: