Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,975,827 C→G 5.2% T253T (ACC→ACG pheA → bifunctional chorismate mutase/prephenate dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,975,8270CG 5.2% 281.0 / 3.6 87T253T (ACC→ACGpheAbifunctional chorismate mutase/prephenate dehydratase
Reads supporting (aligned to +/- strand):  ref base C (39/43);  new base G (2/2);  total (41/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.38e-01

GTGCCGAAAGTGGCTGGACGCTCACTACCCGAACGTGGAGCGCGTGGCCGTTTCGAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTG  >  NC_002947/1975683‑1975974
                                                                                                                                                 |                                                                                                                                                    
gTGCCGAAAGTGGCTGGACGCTCACTACCCGAACGTGGAGCGCGTGGCCGTTTCGAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGCTCGCCGGCGATATGGCGGCCATCCTGTACGGCCTGACCCGTc                                                                                                                                                  >  4:64394/1‑149 (MQ=255)
   ccGAAAGTGGCTGGAAGCTCACTACCCGAACGTGGAGCTCGTGGCCGTTTCGAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGTCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTg                                                                                                                                                <  5:377114/148‑1 (MQ=255)
     gAAAGTGGCTGGACGCTCACTACCCGAACGTGGAGCGCGTGGCCGTTTCGAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGcc                                                                                                                                             >  3:210182/1‑149 (MQ=255)
         gTGGCTGGACGCTCACTACCCGAACGTGGAGCGCGTGGCCGTTTCGAGCAATGCCGAAGCGGTCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCgag                                                                                                                                          >  7:240770/1‑148 (MQ=255)
            gCTGGACGCTCACTACCCGAACGTGGAGCGCGTGGCCGTTTCGAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGagaaga                                                                                                                                      <  3:240498/149‑1 (MQ=255)
            gCTGGACGCTCACTACCCGAACGTGGAGCGCGTGGCCGTTTCGAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCAGACGGGTCAGGCCgtacag                                                                                                                                       <  6:296423/148‑6 (MQ=255)
             cTGGACGCTCACTACCCGAACGTGGAGCGCGTGGCCGTTTCGAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCgag                                                                                                                                          >  4:33616/1‑144 (MQ=255)
             cTGGACGCTCACTACCCGAACGTGGAGCGCGTGGCCGTTTCGAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGAt                                                                                                                                    <  2:117527/149‑1 (MQ=255)
                 aCGCTCACTACCCGAACGTGGAGCGCGTGGCCGTTTCGAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGaa                                                                                                                                <  6:205378/149‑1 (MQ=255)
                 aCGCTCACTACCCGAACGTGGAGCGCGTGGCCGTTTCGAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGaa                                                                                                                                <  1:327474/149‑1 (MQ=255)
                     tCACTACCCGAACGTGGAGCGCGTGGCCGTTTCGAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGAcc                                                                                                                            >  3:335845/1‑149 (MQ=255)
                       aCTACCCGAACGTGGAGCGCGTGGCCGTTTCGAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGc                                                                                                                          >  7:290522/1‑149 (MQ=255)
                        cTACCCGAACGTGGAGCGCGTGGCCGTTTCGAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCAGACGGGTCAGGcc                                                                                                                                             >  2:284647/1‑130 (MQ=255)
                        cTACCCGAACGTGGAGCGCGTGGCCGTTTCGAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCAGACGGGTCAGGcc                                                                                                                                             <  1:284647/130‑1 (MQ=255)
                         tACCCGAACGTGGAGCGCGTGGCCGTTTCGAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGccc                                                                                                                        >  3:2666/1‑149 (MQ=255)
                           cccGAACGTGGAGCGCGTGGCCGTTTCGAGCAATGCCGAAGCGGTCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGccc                                                                                                                        <  8:240770/147‑1 (MQ=255)
                           cccGAACGTGGAGCGCGTGGCCGTTTCGAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCgg                                                                                                                      >  8:298008/1‑149 (MQ=255)
                           cccGAACGTGGAGCGCGTGGCCGTTTCGAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCg                                                                                                                       >  4:198586/1‑148 (MQ=255)
                            ccGAACGTGGAGCGCGTGGCCGTTTCGAGCAATGCCGAAGCGGCCAAGCGTGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGa                                                                                                                     <  8:212163/149‑1 (MQ=255)
                            ccGAACGTGGAGCGCGTGGCCGTTTCGAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGa                                                                                                                     >  6:359956/1‑149 (MQ=255)
                             cGAACGTGGAGCGCGTGGCCGTTTCGAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGa                                                                                                                     >  6:191432/1‑148 (MQ=255)
                             cGAACGTGGAGCGCGTGGCCGTTTCGAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGAc                                                                                                                    <  7:97192/149‑1 (MQ=255)
                             cGAACGTGGAGCGCGTGGCCGTTTCGAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGAc                                                                                                                    <  2:73708/149‑1 (MQ=255)
                             cGAACGTGGAGCGCGTGGCCGTTTCGAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGAGATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGa                                                                                                                     >  8:63628/1‑148 (MQ=255)
                              gAACGTGGAGCGCGTGGCCGTTTCGAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGCAGACCGCCCGGAc                                                                                                                    >  5:196825/1‑148 (MQ=255)
                               aaCGTGGAGCGCGTGGCCGTTTCGAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACaa                                                                                                                  <  3:95454/149‑1 (MQ=255)
                                aCGTGGAGCGCGTGGCCGTTTCGAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAAc                                                                                                                 <  7:173506/149‑1 (MQ=255)
                                aCGTGGAGCGCGTGGCCGTTTCGAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAAc                                                                                                                 <  7:93784/149‑1 (MQ=255)
                                     tagagcgTGGCCGTTTCGTGCAATGCCGATGCGGCCAAGCGGTGCAAGTGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGTCCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCa                                                                                                             <  5:42939/144‑1 (MQ=255)
                                       gcgcGTGGCCGTTTCGAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCAcgc                                                                                                          <  4:306402/149‑1 (MQ=255)
                                        cgcgTGGCCGTTTCGAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCAcgcg                                                                                                         <  8:290522/149‑1 (MQ=255)
                                         gcttGGCCGCTTCGGCATTGCTCGAAACGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACgcgc                                                                                                        >  3:18067/4‑149 (MQ=255)
                                           gTGGCCGTTTCGAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCtt                                                                                                      >  1:286766/1‑149 (MQ=255)
                                               ccGTTTCGAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCtg                                                                                                  <  7:298008/149‑1 (MQ=255)
                                               ccGTTTCGAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCt                                                                                                   >  6:20853/1‑148 (MQ=255)
                                               ccGTTTCGAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCt                                                                                                   >  6:211169/1‑148 (MQ=255)
                                                cGTTTCGAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCtga                                                                                                 >  5:322253/1‑149 (MQ=255)
                                                cGTTTCGAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCtga                                                                                                 >  5:177234/1‑149 (MQ=255)
                                                cGTTTCGAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCtga                                                                                                 <  2:330745/149‑1 (MQ=255)
                                                cGTTTCGAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGACTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCtga                                                                                                 >  5:276337/1‑149 (MQ=255)
                                                  tttCGAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCtgatg                                                                                               <  4:335845/149‑1 (MQ=255)
                                                      gAGCAATGCCGAAGCGGTCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATc                                                                                            <  5:42642/148‑1 (MQ=255)
                                                      gAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATc                                                                                            <  5:280568/148‑1 (MQ=255)
                                                      gAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCg                                                                                           <  2:286766/149‑1 (MQ=255)
                                                        gCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGGCGAGAAGATCGAAGACCGCCCGGTCAACTCCACGCGCTTCCTGATGATCGGt                                                                                         <  2:226986/149‑1 (MQ=255)
                                                          aaTGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTaa                                                                                       >  3:251042/1‑149 (MQ=255)
                                                           aTGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAAc                                                                                      >  3:368246/1‑149 (MQ=255)
                                                             gCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAAc                                                                                      >  8:312443/1‑147 (MQ=255)
                                                                  aaGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGaa                                                                                <  2:289824/149‑1 (MQ=255)
                                                                  aaGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGaa                                                                                <  5:134493/149‑1 (MQ=255)
                                                                          aaGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTAccgcc                                                                        <  6:177234/149‑1 (MQ=255)
                                                                                   aaGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGcga                                                               <  3:198586/149‑1 (MQ=255)
                                                                                   aaGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGcga                                                               <  5:274743/149‑1 (MQ=255)
                                                                                    aGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGcgac                                                              <  8:298752/149‑1 (MQ=255)
                                                                                      ggTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGcgacga                                                            <  4:251042/149‑1 (MQ=255)
                                                                                        tGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACa                                                          <  3:33616/149‑1 (MQ=255)
                                                                                          aGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCgacgac                                                           >  1:49202/1‑146 (MQ=255)
                                                                                          aGTGGAACTCGGCGGCGAGCGCGTGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGAGCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACGGGCGACGACaa                                                         <  3:64394/148‑1 (MQ=255)
                                                                                              gAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCt                                                    >  1:176982/1‑149 (MQ=255)
                                                                                               aaCTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTc                                                   >  7:232560/1‑149 (MQ=255)
                                                                                                 cTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGa                                                 <  2:176982/149‑1 (MQ=255)
                                                                                                          gATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGttt                                        <  6:276337/149‑1 (MQ=255)
                                                                                                              gCCGGCGATATGGCGGCCAACCTGTACGGCCTGACTCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGa                                     <  5:170298/148‑1 (MQ=255)
                                                                                                               ccGGCGATATGGCGGCCAACCTGTACGGGCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTcc                                                                                                    >  1:216163/1‑84 (MQ=255)
                                                                                                               ccGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTcc                                                                                                    <  2:216163/84‑1 (MQ=255)
                                                                                                               ccGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGAt                                    >  7:116743/1‑148 (MQ=255)
                                                                                                               ccGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGAt                                    >  4:370868/1‑148 (MQ=255)
                                                                                                                cGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATTATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGa                                  <  4:2666/149‑1 (MQ=255)
                                                                                                                  gCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGc                                <  6:322253/149‑1 (MQ=255)
                                                                                                                    gATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGcaa                              >  8:272381/1‑149 (MQ=255)
                                                                                                                        tGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGcaacaa                           <  6:196825/148‑1 (MQ=255)
                                                                                                                        tGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGcaacaa                           <  5:359956/148‑1 (MQ=255)
                                                                                                                         ggcggcCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGc                         >  7:105615/1‑149 (MQ=255)
                                                                                                                         ggcggcCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGc                         >  2:56912/1‑149 (MQ=255)
                                                                                                                            ggcCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAg                      <  5:211169/149‑1 (MQ=255)
                                                                                                                               cAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGAGTCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGt                    >  1:247385/1‑149 (MQ=255)
                                                                                                                               cAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTg                   >  7:72758/1‑149 (MQ=255)
                                                                                                                               cAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCCATGAGCAACAAGCCAGGTg                   >  2:204466/1‑149 (MQ=255)
                                                                                                                                  ccTGTACGGCCTGACCCGTCTGGCCGAGTAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGt                >  1:212478/1‑149 (MQ=255)
                                                                                                                                  ccTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGt                >  5:286253/1‑149 (MQ=255)
                                                                                                                                   cTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGtt               >  2:359006/1‑149 (MQ=255)
                                                                                                                                   cTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGtt               >  8:118693/1‑149 (MQ=255)
                                                                                                                                     tctcGGCCAGACGGGTCAGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTg              >  6:154084/4‑148 (MQ=255)
                                                                                                                                      tACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGc             <  6:286253/148‑1 (MQ=255)
                                                                                                                                           ccTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGt       <  5:20853/149‑1 (MQ=255)
                                                                                                                                           ccTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGt       <  4:368246/149‑1 (MQ=255)
                                                                                                                                               aCCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGc    <  1:229322/148‑1 (MQ=255)
                                                                                                                                                cccGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATTATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTg  <  1:78562/149‑1 (MQ=255)
                                                                                                                                                cccGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCt   <  7:312443/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GTGCCGAAAGTGGCTGGACGCTCACTACCCGAACGTGGAGCGCGTGGCCGTTTCGAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTG  >  NC_002947/1975683‑1975974

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: