Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 93,210 G→T 7.5% A128A (GCC→GCA trpB ← tryptophan synthase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294793,2100GT 7.5% 551.9 / 24.1 169A128A (GCC→GCAtrpBtryptophan synthase subunit beta
Reads supporting (aligned to +/- strand):  ref base G (86/70);  new base T (6/7);  total (92/77)
Fisher's exact test for biased strand distribution p-value = 5.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 195 total aligned reads displayed.

TGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCG  >  NC_002947/93062‑93358
                                                                                                                                                    |                                                                                                                                                    
tGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAgtgg                                                                                                                                                     >  1:3303982/1‑150 (MQ=255)
 cgCGTCACTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCCAGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAgtggc                                                                                                                                                    <  2:3257484/149‑1 (MQ=255)
  gCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGACGCGCGCGGCGACAgtggcg                                                                                                                                                   >  2:489149/1‑150 (MQ=255)
       tttCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGATCACGTTCGCCTGCTGGCGCTCGCTGTCGGTGGCACCCATGTAGATCTCGCAAGGCAGGCCGAAGCGCGCGGCGACAgtggcggtggc                                                                                                                                              >  1:610411/1‑150 (MQ=255)
       tttCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAgtggcggtggc                                                                                                                                              >  1:1020967/1‑150 (MQ=255)
       tttCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAgtggcggtggc                                                                                                                                              <  2:1570/150‑1 (MQ=255)
       tttCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAgtggcggtggc                                                                                                                                              >  1:3492401/1‑150 (MQ=255)
        ttCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAggggcggtggcg                                                                                                                                             >  2:4461640/1‑150 (MQ=255)
           aGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGaca                                                                                                                                          >  2:1202865/1‑150 (MQ=255)
            gggTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGacac                                                                                                                                         >  1:3901461/1‑150 (MQ=255)
               tACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGt                                                                                                                                      <  1:2247505/150‑1 (MQ=255)
                   ggTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCAt                                                                                                                                  <  1:5318542/150‑1 (MQ=255)
                   ggTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCAt                                                                                                                                  <  1:5318425/150‑1 (MQ=255)
                       ccGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCg                                                                                                                              <  2:1674948/150‑1 (MQ=255)
                         ggcggTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGgc                                                                                                                            >  2:1778446/1‑150 (MQ=255)
                            ggTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGcc                                                                                                                         <  2:4082154/150‑1 (MQ=255)
                                  tGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTc                                                                                                                   >  1:960933/1‑150 (MQ=255)
                                     cACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGc                                                                                                                <  1:1601233/150‑1 (MQ=255)
                                      tttacttcGGCACCCAGCAGCTTCATGCGGAAGACGTTTGGCTTCTGACGTTCAACGTCAGTGGCACCCATATAGACCACGCAAGGCAGGCCGAAACGGGCTGCAACGGTTGCTGATGCAACGCCATGCTGACCAGCGCCGGTTTCcgcg                                                                                                               >  2:946012/7‑150 (MQ=255)
                                      aCGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGTGCCGGTTTCGGCg                                                                                                               <  2:1834480/150‑1 (MQ=255)
                                      aCGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCg                                                                                                               >  1:1399276/1‑150 (MQ=255)
                                       cGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGa                                                                                                              <  2:2254889/150‑1 (MQ=255)
                                        gATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCAGTTTCGGCGAt                                                                                                             >  2:2559112/1‑150 (MQ=255)
                                         aTCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGGCGGTGGCACCCATGTAGATCACGCAAGGCCGGCCGACGCGCGCGGCGACCGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATc                                                                                                            >  2:298255/1‑150 (MQ=255)
                                          tctcGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCa                                                                                                           <  2:5380206/150‑1 (MQ=255)
                                           ttcGGCACCCAGCAGCTTCATGCGGAAGACGTTTGGCTTCTGACGTTCAACGTCAGTGGCACCCATATAGACCACGCAAGGCAGGCCGAAACGGGCTGCAACGGTTGCTGATGCAACGCCATGCTGACCAGCGCCGGTTTCCGCGAtgat                                                                                                          >  1:4534855/2‑147 (MQ=255)
                                           ctcGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAg                                                                                                          <  1:1120538/150‑1 (MQ=255)
                                                cgcCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCt                                                                                                     >  1:2571897/1‑150 (MQ=255)
                                                    cagcagCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCttttt                                                                                                 <  2:127624/150‑1 (MQ=255)
                                                    cagcagCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCttttt                                                                                                 >  1:2648231/1‑150 (MQ=255)
                                                    cagcagCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCttttt                                                                                                 <  2:4409882/150‑1 (MQ=255)
                                                    cagcagCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCttttt                                                                                                 <  1:3999072/150‑1 (MQ=255)
                                                     agcagcTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTg                                                                                                <  1:2143037/150‑1 (MQ=255)
                                                          cTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCAt                                                                                           <  1:1094864/150‑1 (MQ=255)
                                                           ttCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATg                                                                                          >  1:3198752/1‑150 (MQ=255)
                                                              aTGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATgcgc                                                                                       >  1:3825468/1‑150 (MQ=255)
                                                               tGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCa                                                                                                           >  2:688068/1‑129 (MQ=255)
                                                               tGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCa                                                                                                           <  1:688068/129‑1 (MQ=255)
                                                               tGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCt                                                                                      >  1:2354733/1‑150 (MQ=255)
                                                                gCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCtt                                                                                     <  2:3538905/150‑1 (MQ=255)
                                                                gCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCtt                                                                                     <  1:774676/150‑1 (MQ=255)
                                                                      cacGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCcag                                                                               <  2:5222851/150‑1 (MQ=255)
                                                                        cGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCcagca                                                                             <  1:194442/150‑1 (MQ=255)
                                                                         gTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTg                                                                                                >  2:137408/1‑130 (MQ=255)
                                                                         gTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCcagcag                                                                            >  2:2912915/1‑150 (MQ=255)
                                                                            cGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCAc                                                                         >  2:4014297/1‑150 (MQ=255)
                                                                            cGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCAc                                                                         <  1:4184727/150‑1 (MQ=255)
                                                                                  ctgGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGcc                                                                   <  1:3383445/150‑1 (MQ=255)
                                                                                   tgGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCg                                                                  >  1:548976/1‑150 (MQ=255)
                                                                                     gcgcTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGAt                                                                >  2:1079777/1‑150 (MQ=255)
                                                                                     gcgcTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGAt                                                                >  2:4247756/1‑150 (MQ=255)
                                                                                     gcgcTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGAt                                                                >  2:3820204/1‑150 (MQ=255)
                                                                                       gcTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGc                                                              >  2:4829034/1‑150 (MQ=255)
                                                                                       gcTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGc                                                              >  2:1791305/1‑150 (MQ=255)
                                                                                        cTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCa                                                             >  2:760183/1‑150 (MQ=255)
                                                                                        cTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCa                                                             <  1:2427726/150‑1 (MQ=255)
                                                                                        cTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCa                                                             >  2:3913644/1‑150 (MQ=255)
                                                                                           gATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAgtt                                                          >  1:4682546/1‑150 (MQ=255)
                                                                                           aacgTCAGTGGCACCCATATAGACCACGCAAGGCAGGCCGAAACGGGCTGCAACGGTTGCTGATGCAACGCCATGCTGACCAGCGCCGGTTTCCGCGATGATGCGGGTTTTGCCCATGCGCTTGGCAAGCAGGATCTGGCCGAGGCAatt                                                          <  1:1296038/147‑1 (MQ=255)
                                                                                               tcagTGGCACCCATATAGACCACGCAAGGCAGGCCGAAACGGGCTGCAACGGTTGCTGATGCAACGCCATGCTGACCAGCGCCGGTTTCCGCGATGATGCGGGTTTTGCCCATGCGCTTGGCAAGCAGGATCTGGCCGAGGCAattgtta                                                      <  2:3024314/147‑2 (MQ=255)
                                                                                                  gTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCAcc                                                   <  2:1686963/150‑1 (MQ=255)
                                                                                                  gTGGCACCCATATAGACCACGCAAGGCAGGCCGAAACGGGCTGCAACGGTTGCTGATGCAACGCCATGCTGACCAGCGCCGGTTTCCGCGATGATGCGGGTTTTGCCCATGCGCTTGGCAAGCAGGATCTGGCCGAGGCAattgttaatc                                                   >  1:3459595/1‑146 (MQ=255)
                                                                                                   tGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCt                                                  >  1:1312396/1‑150 (MQ=255)
                                                                                                      cACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTtgt                                               <  1:1206719/150‑1 (MQ=255)
                                                                                                        cccATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTgc                                             <  1:1712426/150‑1 (MQ=255)
                                                                                                         ccATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCAAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTgcg                                            >  2:4885444/1‑150 (MQ=255)
                                                                                                          cATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTgcgc                                           >  1:2767397/1‑150 (MQ=255)
                                                                                                           aTGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTgcgcg                                          >  2:1735634/1‑150 (MQ=255)
                                                                                                           aTGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTgcgcg                                          >  2:3321481/1‑150 (MQ=255)
                                                                                                              tAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCg                                       <  2:4448358/150‑1 (MQ=255)
                                                                                                                gATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGt                                     <  2:3651322/150‑1 (MQ=255)
                                                                                                                gATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGt                                     <  2:3577549/150‑1 (MQ=255)
                                                                                                                gATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGt                                     >  2:577965/1‑150 (MQ=255)
                                                                                                                   cACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAgttg                                                         >  1:3342989/1‑127 (MQ=255)
                                                                                                                   cACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAgttg                                                         <  2:3342989/127‑1 (MQ=255)
                                                                                                                   cACGCAAGGCAGGCCGAAACGGGCTGCAACGGTTGCTGATGCAACGCCATGCTGACCAGCGCCGGTTTCCGCGATGATGCGGGTTTTGCCCATGCGCTTGGCAAGCAGGATCTGGCCGAGGCAATTGTTAATCTTGTGGCTGCCGGTATg                                  <  1:4405528/150‑1 (MQ=255)
                                                                                                                     cGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGt                                >  2:2192097/1‑150 (MQ=255)
                                                                                                                      gCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGtt                               >  1:4276336/1‑150 (MQ=255)
                                                                                                                      gCAAGGCAGGCCGAAACGGGCTGCAACGGTTGCTGATGCAACGCCATGCTGACCAGCGCCGGTTTCCGCGATGATGCGGGTTTTGCCCATGCGCTTGGCAAGCAGGATCTGGCCGAGGCAATTGTTAATCTTGTGGCTGCCGGTATGGtt                               >  1:2948899/1‑150 (MQ=255)
                                                                                                                        aaGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGa                             >  2:1520583/1‑150 (MQ=255)
                                                                                                                        aaGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTCGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGa                             <  1:507703/150‑1 (MQ=255)
                                                                                                                            caggcCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCTGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGctc                         >  2:2363318/1‑150 (MQ=255)
                                                                                                                            caggcCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGctc                         >  2:3167655/1‑150 (MQ=255)
                                                                                                                             aggcCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGCTGAGctct                        >  1:1372309/1‑150 (MQ=255)
                                                                                                                               gcCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCgg                                      <  2:1262762/134‑1 (MQ=255)
                                                                                                                                ccGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTg                                                                                    >  2:880294/1‑87 (MQ=255)
                                                                                                                                ccGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTg                                                                                    <  1:880294/87‑1 (MQ=255)
                                                                                                                                ccGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCa                     <  2:668709/150‑1 (MQ=255)
                                                                                                                                   aaGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGc                  <  1:2707144/150‑1 (MQ=255)
                                                                                                                                     gcgcgcGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCtt                <  2:2672009/150‑1 (MQ=255)
                                                                                                                                     gcgcgcGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCtt                >  2:1655942/1‑150 (MQ=255)
                                                                                                                                       gcgcgGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTga              <  1:3313713/150‑1 (MQ=255)
                                                                                                                                        cgcgGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTgaa             <  2:3906596/150‑1 (MQ=255)
                                                                                                                                           ggcgACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCTGTATGGTTGAGCTCTTCACGCTTgaagaa          <  1:1165575/150‑1 (MQ=255)
                                                                                                                                              gACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAAt       >  2:2538564/1‑150 (MQ=255)
                                                                                                                                                 aGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCtt    >  1:4525807/1‑150 (MQ=255)
                                                                                                                                                 aGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCtt    >  1:4525732/1‑150 (MQ=255)
                                                                                                                                                  gtggcggtggcgACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTc   <  1:3957467/150‑1 (MQ=255)
                                                                                                                                                   tggcggtggcgACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTg                                                                                    >  2:2161446/1‑68 (MQ=255)
                                                                                                                                                   tggcggtggcgACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCg  >  2:5346256/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCG  >  NC_002947/93062‑93358

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: