Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 93,176 T→C 12.7% I140V (ATC→GTC)  trpB ← tryptophan synthase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294793,1760TC12.7% 488.7 / 46.3 181I140V (ATC→GTC) trpBtryptophan synthase subunit beta
Reads supporting (aligned to +/- strand):  ref base T (82/76);  new base C (13/10);  total (95/86)
Fisher's exact test for biased strand distribution p-value = 8.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 182 total aligned reads displayed.

GACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTA  >  NC_002947/93027‑93324
                                                                                                                                                     |                                                                                                                                                    
gACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGAt                                                                                                                                                      <  2:2761657/150‑1 (MQ=255)
 aCGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATc                                                                                                                                                     <  2:1610536/150‑1 (MQ=255)
  cGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCa                                                                                                                                                    <  1:3246512/150‑1 (MQ=255)
   gTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCAc                                                                                                                                                   <  1:3724062/150‑1 (MQ=255)
    ttGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACg                                                                                                                                                  <  1:1004914/150‑1 (MQ=255)
     tGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGc                                                                                                                                                 <  2:312541/150‑1 (MQ=255)
       gtgaCCCAATCACGAAGCGCTTCGTTCATGGCGTCTTTGAGCGTGCCATTGCCAGCCGTTACCGGCTTTACTTCGGCACCCAGCAGCTTCATGCGGAAGACGTTTGGCTTCTGACGTTCAACGTCAGTGGCACCCATATAGACCACGCaa                                                                                                                                               <  1:3836165/147‑1 (MQ=255)
       gTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCaa                                                                                                                                               >  1:4572710/1‑150 (MQ=255)
       gTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCaa                                                                                                                                               <  2:3822530/150‑1 (MQ=255)
         aaCCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAagg                                                                                                                                             <  2:1930346/150‑1 (MQ=255)
           cccAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAaggca                                                                                                                                           >  2:4341829/1‑150 (MQ=255)
            ccAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAaggcag                                                                                                                                          <  1:2601599/150‑1 (MQ=255)
             caatCACGAAGCGCTTCGTTCATGGCGTCTTTGAGCGTGCCATTGCCAGCCGTTACCGGCTTTACTTCGGCACCCAGCAGCTTCATGCGGAAGACGTTTGGCTTCTGACGTTCAACGTCAGTGGCACCCATATAGACCACGCAaggcagg                                                                                                                                         >  1:1335683/4‑150 (MQ=255)
             caatCACGAAGCGCTTCGTTCATGGCGTCTTTGAGCGTGCCATTGCCAGCCGTTACCGGCTTTACTTCGGCACCCAGCAGCTTCATGCGGAAGACGTTTGGCTTCTGACGTTCAACGTCAGTGGCACCCATATAGACCACGCAaggcagg                                                                                                                                         >  1:3915333/4‑150 (MQ=255)
              aatCACGAAGCGATTCGTTCATGGCGTCTTTGAGCGTGCCATTGCCAGCCGTTACCGTCTTTACTTCGGCACCCAGCAGCTTCATGCGGAAGACGTTTGGATTCTGACGTTCAACGTCAGTGGCACCCATATAGACCACGCAaggcaggc                                                                                                                                        <  1:1973850/148‑1 (MQ=255)
               gTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGcc                                                                                                                                       <  1:3423204/150‑1 (MQ=255)
                  aCGAAGCGCTTCGTTCATGGCGTCTTTGAGCGTGCCATTGCCAGCCGTTACCGGCTTTACTTCGGCACCCAGCAGCTTCATGCGGAAGACGTTTGGCTTCTGACGTTCAACGTCAGTGGCACCCATATAGACCACGCAAGGCAGGCCGaa                                                                                                                                    >  2:3653837/1‑150 (MQ=255)
                    gCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAgc                                                                                                                                  >  2:296175/1‑150 (MQ=255)
                         gCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGcgcgcg                                                                                                                             >  1:1808243/1‑150 (MQ=255)
                           ttCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCgcggc                                                                                                                           >  2:508078/1‑150 (MQ=255)
                            tCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCgcggcg                                                                                                                          >  2:2246241/1‑150 (MQ=255)
                               ttCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACa                                                                                                                       >  1:3651796/1‑150 (MQ=255)
                               gtCATGGCGTCTTTCAGGGGACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCACGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCTTGTACATCACGCAAGGCAGGCCGTAGCGCGCGGCGACa                                                                                                                       <  1:1026776/149‑1 (MQ=255)
                                  aTGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAgtg                                                                                                                    <  2:269824/150‑1 (MQ=255)
                                  aTGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAgtg                                                                                                                    >  2:3426551/1‑150 (MQ=255)
                                  aTGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAgtg                                                                                                                    >  2:641747/1‑150 (MQ=255)
                                   tGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAgtgg                                                                                                                   >  1:4865055/1‑150 (MQ=255)
                                     gCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAgtggcg                                                                                                                 >  1:4646778/1‑150 (MQ=255)
                                       gTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAgtggcggt                                                                                                               >  2:73447/1‑150 (MQ=255)
                                       gTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAgtggcggt                                                                                                               <  2:536113/150‑1 (MQ=255)
                                        tCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAgtggcggtg                                                                                                              <  1:5085746/150‑1 (MQ=255)
                                          tttCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAgtggcggtggc                                                                                                            >  2:4136700/1‑150 (MQ=255)
                                          tttCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAgtggcggtggc                                                                                                            >  2:2370387/1‑150 (MQ=255)
                                          tttCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAgtggcggtggc                                                                                                            >  1:1506447/1‑150 (MQ=255)
                                          tttCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAgtggcggtggc                                                                                                            >  1:4357228/1‑150 (MQ=255)
                                            tCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGa                                                                                                          <  1:3772555/150‑1 (MQ=255)
                                            tCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGa                                                                                                          >  2:3644992/1‑150 (MQ=255)
                                            tCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGa                                                                                                          >  1:864997/1‑150 (MQ=255)
                                               gggTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGacac                                                                                                       >  1:2380035/1‑150 (MQ=255)
                                               gggTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGacac                                                                                                       <  1:5103121/150‑1 (MQ=255)
                                                  tACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGt                                                                                                    <  2:2581484/150‑1 (MQ=255)
                                                  tACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGt                                                                                                    <  2:2525099/150‑1 (MQ=255)
                                                  tACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGt                                                                                                    <  2:2667805/150‑1 (MQ=255)
                                                              cggTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGgcgc                                                                                        <  2:4912979/150‑1 (MQ=255)
                                                                    cTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGttt                                                                                  >  1:970956/1‑150 (MQ=255)
                                                                       gCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCgg                                                                               <  1:4880423/150‑1 (MQ=255)
                                                                       gCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCgg                                                                               <  2:1214260/150‑1 (MQ=255)
                                                                       gCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCgg                                                                               <  1:2667725/150‑1 (MQ=255)
                                                                         aCGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCg                                                                             >  2:39243/1‑150 (MQ=255)
                                                                          cGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGa                                                                            <  1:177969/150‑1 (MQ=255)
                                                                              ttcGGCACCCAGCAGCTTCATGCGGAAGACGTTTGGCTTCTGACGTTCAACGTCAGTGGCACCCATATAGACCACGCAAGGCAGGCCGAAACGGGCTGCAACGGTTGCTGATGCAACGCCATGCTGACCAGCGCCGGTTTCCGCGAtgat                                                                        >  2:3332227/2‑147 (MQ=255)
                                                                              ttcGGCACCCAGCAGCTTCATGCGGAAGACGTTTGGCTTCTGACGTTCAACGTCAGTGGCACCCATATAGACCACGCAAGGCAGGCCGAAACGGGCTGCAACGGTTGCTGATGCAACGCCATGCTGACCAGCGCCGGTTTCCGCGAtgat                                                                        >  1:2730755/2‑147 (MQ=255)
                                                                               tcGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAgg                                                                       <  2:1342423/150‑1 (MQ=255)
                                                                                  gcacccAGCAGCTTCATGCGGAAGACGTTTGGCTTCTGACGTTCAACGTCAGTGGCACCCATATAGACCACGCAAGGCAGGCCGAAACGGGCTGCAACGGTTGCTGATGCAACGCCATGCTGACCAGCGCCGGTTTCCGCGATGATgcgg                                                                    <  2:506099/147‑2 (MQ=255)
                                                                                  gcacccAGCAGCTTCATGCGGAAGACGTTTGGCTTCTGACGTTCAACGTCAGTGGCACCCATATAGACCACGCAAGGCAGGCCGAAACGGGCTGCAACGGTTGCTGATGCAACGCCATGCTGACCAGCGCCGGTTTCCGCGATGATgcgg                                                                    <  2:4905004/147‑2 (MQ=255)
                                                                                  gcacccAGCAGCTTCATGCGGAAGACGTTTGGCTTCTGACGTTCAACGTCAGTGGCACCCATATAGACCACGCAAGGCAGGCCGAAACGGGCTGCAACGGTTGCTGATGCAACGCCATGCTGACCAGCGCCGGTTTCCGCGATGATgcgg                                                                    >  2:2993748/4‑149 (MQ=255)
                                                                                      ccAGCAGCTTCATGCGGAAGACGTTTGGCTTCTGACGTTCAACGTCAGTGGCACCCATATAGACCACGCAAGGCAGGCCGAAACGGGCTGCAACGGTTGCTGATGCAACGCCATGCTGACCAGCGCCGGTTTCCGCGATGATgcgggttt                                                                >  2:2637736/1‑145 (MQ=255)
                                                                                      ccAGCAGCTTCATGCGGAACATGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCtttt                                                                >  2:1115708/1‑150 (MQ=255)
                                                                                      ccAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCtttt                                                                >  1:1279656/1‑150 (MQ=255)
                                                                                        agcagcTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTg                                                              >  1:692664/1‑150 (MQ=255)
                                                                                             ctccATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCAt                                                         <  2:3774971/147‑1 (MQ=255)
                                                                                             cTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCAt                                                         >  1:457581/1‑150 (MQ=255)
                                                                                              ttCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCg                                                                                            <  2:2498351/114‑1 (MQ=255)
                                                                                              ttCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCg                                                                                            >  1:2498351/1‑114 (MQ=255)
                                                                                              ttCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATg                                                        >  2:90038/1‑150 (MQ=255)
                                                                                               tCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATgc                                                       <  2:1398270/150‑1 (MQ=255)
                                                                                                 aTGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATgcgc                                                     >  2:4790196/1‑150 (MQ=255)
                                                                                                   gCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCtt                                                   <  2:4094128/150‑1 (MQ=255)
                                                                                                    cGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTg                                                  <  1:2097962/150‑1 (MQ=255)
                                                                                                      gAAGACGTTTGGCTTCTGACGTTCAACGTCAGTGGCACCCATATAGACCACGCAAGGCAGGCCGAAACGGGCTGCAACGGTTGCTGATGCAACGCCATGCTGACCAGCGCCGGTTTCCGCGATGATGCGGGTTTTGCCCATGCGCTTGGc                                                <  2:3708381/150‑1 (MQ=255)
                                                                                                      gAACATGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGc                                                >  1:2030292/1‑150 (MQ=255)
                                                                                                        acacGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCca                                              >  1:3976381/1‑150 (MQ=255)
                                                                                                         cacGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCcag                                             <  1:611002/150‑1 (MQ=255)
                                                                                                           cGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCcagca                                           <  2:3157957/150‑1 (MQ=255)
                                                                                                                    gctgGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGc                                  <  2:614322/150‑1 (MQ=255)
                                                                                                                    gctgGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGc                                  <  1:3632405/150‑1 (MQ=255)
                                                                                                                     ctgGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGcc                                 >  2:2119923/1‑150 (MQ=255)
                                                                                                                      tgGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCg                                >  2:3083449/1‑150 (MQ=255)
                                                                                                                        gcgcTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGAt                              <  2:1101857/150‑1 (MQ=255)
                                                                                                                          gcTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGc                            >  2:4218258/1‑150 (MQ=255)
                                                                                                                          gcTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGc                            >  1:2801572/1‑150 (MQ=255)
                                                                                                                            tCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGGCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTTCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAg                          <  1:2804953/150‑1 (MQ=255)
                                                                                                                            tCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAg                          >  1:4248354/1‑150 (MQ=255)
                                                                                                                            tCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAg                          <  2:4272835/150‑1 (MQ=255)
                                                                                                                            tCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAg                          <  2:4273282/150‑1 (MQ=255)
                                                                                                                            tCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAg                          <  1:2805018/150‑1 (MQ=255)
                                                                                                                              gATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAgtt                        >  1:2702052/1‑150 (MQ=255)
                                                                                                                              gATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAgtt                        >  2:2037867/1‑150 (MQ=255)
                                                                                                                              aacgTCAGTGGCACCCATATAGACCACGCAAGGCAGGCCGAAACGGGCTGCAACGGTTGCTGATGCAACGCCATGCTGACCAGCGCCGGTTTCCGCGATGATGCGGGTTTTGCCCATGCGCTTGGCAAGCAGGATCTGGCCGAGGCAatt                        <  1:3332227/147‑1 (MQ=255)
                                                                                                                                tGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAgttgt                      >  1:229906/1‑150 (MQ=255)
                                                                                                                                    ggTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCAc                  <  2:794181/150‑1 (MQ=255)
                                                                                                                                     gTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCAcc                 >  2:3836833/1‑150 (MQ=255)
                                                                                                                                     gTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCAcc                 >  1:2567118/1‑150 (MQ=255)
                                                                                                                                         cACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTtgt             <  2:4903797/150‑1 (MQ=255)
                                                                                                                                            ccATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTgcg          >  2:3411730/1‑150 (MQ=255)
                                                                                                                                               tGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTgcgcgc       <  2:4028166/150‑1 (MQ=255)
                                                                                                                                                  aGACCACGCAAGGCAGGCCGAAACGGGCTGCAACGGTTGCTGATGCAACGCCATGCTGACCAGCGCCGGTTTCCGCGATGATGCGGGTTTTGCCCATGCGCTTGGCAAGCAGGATCTGGCCGagg                             >  2:4722650/1‑123 (MQ=255)
                                                                                                                                                   gATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGt   >  2:672269/1‑150 (MQ=255)
                                                                                                                                                   gATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGt   >  2:1850646/1‑150 (MQ=255)
                                                                                                                                                    aTCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCGCCTTGTGCGCGCCGGTa  >  2:2640750/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTA  >  NC_002947/93027‑93324

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