Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,851,606 C→T 12.0% H111H (CAC→CAT ilvD → dihydroxy‑acid dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,851,6060CT12.0% 652.4 / 60.2 244H111H (CAC→CATilvDdihydroxy‑acid dehydratase
Reads supporting (aligned to +/- strand):  ref base C (111/102);  new base T (18/11);  total (130/114)
Fisher's exact test for biased strand distribution p-value = 3.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 245 total aligned reads displayed.

GAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGACCAAGCTGGCCAGCCACGGCCTGGA  >  NC_002947/5851460‑5851755
                                                                                                                                                  |                                                                                                                                                     
gAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCTGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTgc                                                                                                                                                    >  2:2917708/1‑150 (MQ=255)
  aaTCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTgcgc                                                                                                                                                  >  2:896158/1‑150 (MQ=255)
   aTCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGcc                                                                                                                                                 <  2:1386717/150‑1 (MQ=255)
    tCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGACTACATGGTCAATGCGCACTGCGCCg                                                                                                                                                >  2:4727274/1‑150 (MQ=255)
      gAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGAc                                                                                                                                              <  1:5353919/150‑1 (MQ=255)
      gAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGAc                                                                                                                                              <  2:2712821/150‑1 (MQ=255)
       aaCGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACg                                                                                                                                             >  1:4291689/1‑150 (MQ=255)
             ccGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCg                                                                                                                                       <  1:1021278/150‑1 (MQ=255)
             ccGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCg                                                                                                                                       >  2:3428386/1‑150 (MQ=255)
              cGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGAAGCCATCGt                                                                                                                                      <  1:2290805/150‑1 (MQ=255)
                  ggcgtggcCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGt                                                                                                                                  <  1:501804/150‑1 (MQ=255)
                  ggcgtggcCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGt                                                                                                                                  <  1:2501513/150‑1 (MQ=255)
                     gttgcCAAGGAATTCAACACCATTGCGGTCGATGACGGTATCGCGATGGGCCATGACGGCATGCTCTATTCGCTGCCATCGCGCGAGCTGATTTCCGACTCGGTTGAATATATGGTCAATGCCCATTGCGCCGACGCGATGGTCTGCAtt                                                                                                                               <  2:3466051/147‑2 (MQ=255)
                      tggcCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCAGTgcgcattgac                                                                                                                                            >  2:3186631/1‑130 (MQ=255)
                      tggcCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTAtct                                                                                                                              <  1:3108021/150‑1 (MQ=255)
                       ggcCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTAtctc                                                                                                                             <  1:3939924/150‑1 (MQ=255)
                        gcCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTcc                                                                                                                            >  1:1561621/1‑150 (MQ=255)
                         ccAAGGAATTCAACACCATTGCGGTCGATGACGGTATCGCGATGGGCCATGACGGCATGCTCTATTCGCTGCCATCGCGCGAGCTGATTTCCGACTCGGTTGAATATATGGTCAATGCCCATTGCGCCGACGCGATGGTCTGCATTTCCa                                                                                                                           >  1:4895703/1‑150 (MQ=255)
                         ccAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCa                                                                                                                           >  1:2371923/1‑150 (MQ=255)
                          cAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCaa                                                                                                                          >  1:1735020/1‑150 (MQ=255)
                          cAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCaa                                                                                                                          <  2:3775550/150‑1 (MQ=255)
                           aaGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAACTCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCGGTCTCTCCagc                                                                                                                         >  1:4252971/1‑148 (MQ=255)
                           aaGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAAc                                                                                                                         >  1:3989257/1‑150 (MQ=255)
                           aaGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAAc                                                                                                                         <  1:3584211/150‑1 (MQ=255)
                              gAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGc                                                                                                                      <  1:2778654/150‑1 (MQ=255)
                                 ttCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGATGCCGTGGAATACAGGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGAc                                                                                                                   >  2:3115610/1‑150 (MQ=255)
                                    aaCACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAg                                                                                                                >  2:3416449/1‑150 (MQ=255)
                                    aaCACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAg                                                                                                                <  1:4407549/150‑1 (MQ=255)
                                       ccccccgcgGCCCAGGACGCCAGCGCCACGGGCCACGCCGGCAGGCTGGACTCGCGGCCAAGCCGCGACCTCCTTGCCGACGCCGTGGAATACATGGTCACTGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATc                                                                                                             <  2:3249215/145‑1 (MQ=255)
                                         cATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACCCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCAc                                                                                                           <  1:1859869/150‑1 (MQ=255)
                                          aTCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCAcc                                                                                                          >  2:1715489/1‑150 (MQ=255)
                                             gcgGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCAccccc                                                                                                       >  1:4333412/1‑150 (MQ=255)
                                               ggTCGATGACGGTATCGCGATGGGCCATGACGGCATGCTCTATTCGCTGCCATCGCGCGAGCTGATTTCCGACTCGGTTGAATATATGGTCAATGCCCATTGCGCCGACGCGATGGTCTGCATTTCCAACTGCGACAAGATCACTCCCgg                                                                                                     >  1:119/1‑150 (MQ=255)
                                                  cGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCAt                                                                                                  <  2:5171544/150‑1 (MQ=255)
                                                  cGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCAt                                                                                                  <  2:4549498/150‑1 (MQ=255)
                                                    aTGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGc                                                                                                >  2:5363595/1‑150 (MQ=255)
                                                         ggCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGt                                                                                                                                      <  1:3030964/107‑1 (MQ=255)
                                                            aTCGCGATGGGCCATGACGGCATGCTCTATTCGCTGCCATCGCGCGAGCTGATTTCCGACTCGGTTGAATATATGGTCAATGCCCATTGCGCCGACGCGATGGTCTGCATTTCCAACTGCGACAAGATCACTCCCGGCATGCTGATGgca                                                                                        <  2:4052093/150‑2 (MQ=255)
                                                            aTCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGgcc                                                                                        >  2:3159401/1‑150 (MQ=255)
                                                                cgaTGGGCCATGACGGCATGCTCTATTCGCTGCCATCGCGCGAGCTGATTTCCGACTCGGTTGAATATATGGTCAATGCCCATTGCGCCGACGCGATGGTCTGCATTTCCAACTGCGACAAGATCACTCCCGGCATGCTGATGgcagcac                                                                                    >  2:4302630/3‑148 (MQ=255)
                                                                ccATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGccc                                                                                    >  1:1701168/1‑150 (MQ=255)
                                                                  aTGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGa                                                                                                               <  1:584539/121‑1 (MQ=255)
                                                                     ggCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCATCTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTgcgc                                                                               >  2:457998/1‑150 (MQ=255)
                                                                     ggCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTgcgc                                                                               >  2:5276845/1‑150 (MQ=255)
                                                                       ccACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCt                                                                             <  1:3807758/150‑1 (MQ=255)
                                                                        cACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTg                                                                            >  1:1145783/1‑150 (MQ=255)
                                                                           gacgGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAAc                                                                         <  1:4968924/150‑1 (MQ=255)
                                                                           gacgGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAAc                                                                         >  1:64272/1‑150 (MQ=255)
                                                                                cATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATccc                                                                    >  2:3706396/1‑150 (MQ=255)
                                                                                 aTGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCg                                                                   <  1:2108601/150‑1 (MQ=255)
                                                                                 aTGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCg                                                                   <  2:2820247/150‑1 (MQ=255)
                                                                                 aTGCTCTATTCGCTGCCATCGCGCGAGCTGATTTCCGACTCGGTTGAATATATGGTCAATGCCCATTGCGCCGACGCGATGGTCTGCATTTCCAACTGCGACAAGATCACTCCCGGCATGCTGATGGCAGCACTCCGCCTCAATATCCCg                                                                   <  2:5105939/150‑1 (MQ=255)
                                                                                   gCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGt                                                                 <  2:1961143/150‑1 (MQ=255)
                                                                                   gCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGt                                                                 >  1:1487970/1‑150 (MQ=255)
                                                                                     tGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGa                                                               >  1:462901/1‑150 (MQ=255)
                                                                                     tGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGa                                                               <  2:4366195/150‑1 (MQ=255)
                                                                                         cTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCtt                                                           >  2:3946881/1‑150 (MQ=255)
                                                                                          tCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTc                                                          <  1:4522177/150‑1 (MQ=255)
                                                                                          tCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTc                                                          <  2:247858/150‑1 (MQ=255)
                                                                                           cGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCg                                                         <  1:645670/150‑1 (MQ=255)
                                                                                           cGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCg                                                         <  2:1633461/150‑1 (MQ=255)
                                                                                                ccAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTcc                                                    <  2:1865496/150‑1 (MQ=255)
                                                                                                ccAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTcc                                                    >  1:4121323/1‑150 (MQ=255)
                                                                                                 cAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCcg                                                   >  1:718563/1‑150 (MQ=255)
                                                                                                  aaGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCcgg                                                  >  1:3250791/1‑150 (MQ=255)
                                                                                                   aGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCcggc                                                 <  2:2464588/150‑1 (MQ=255)
                                                                                                       gcgAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGcc                                                                                                                                                 <  1:4791501/50‑1 (MQ=255)
                                                                                                       gcgAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGcc                                                                                                                                                 >  2:4791501/1‑50 (MQ=255)
                                                                                                        tgataTCTTTGGCGAAGCCGTGGTATTCATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGccc                                            <  2:4927375/146‑1 (MQ=255)
                                                                                                          aaaTCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGa                                          <  1:4156069/150‑1 (MQ=255)
                                                                                                           aaTCAGTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGAt                                         <  1:3688929/150‑1 (MQ=255)
                                                                                                             tCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATgg                                       <  1:3488891/150‑1 (MQ=255)
                                                                                                               aTTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGaa                                     <  2:2307805/150‑1 (MQ=255)
                                                                                                                 tGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGc                                   <  1:2917708/150‑1 (MQ=255)
                                                                                                                  gCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCAGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGcc                                  >  2:391192/1‑150 (MQ=255)
                                                                                                                   ccGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCg                                 >  1:3186178/1‑150 (MQ=255)
                                                                                                                   ccGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCg                                 <  2:2371923/150‑1 (MQ=255)
                                                                                                                    cGACTCGGTTGAATATATGGTCAATGCCCATTGCGCCGACGCGATGGTCTGCATTTCCAACTGCGACAAGATCACTCCCGGCATGCTGATGGCAGCACTCCGCCTCAATATCCCGGTCGTGTTCGTTTCGGGCGGGCCGATGGAAGCAgg                                <  1:3652076/150‑1 (MQ=255)
                                                                                                                       cGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCaa                             <  1:1423691/150‑1 (MQ=255)
                                                                                                                        tcggttgAATATATGGTCAATGCCCATTGCGCCGACGCGATGGTCTGCATTTCCAACTGCGACAAGATCACTCCCGGCATGCTGATGGCAGCACTCCGCCTCAATATCCCGGTCGTGTTCGTTTCGGGCGGGCCGATGGACGCAGGCAAg                            >  2:4150810/7‑150 (MQ=255)
                                                                                                                        tcggttgAATATATGGTCAATGCCCATTGCGCCGACGCGATGGTCTGCATTTCCAACTGCGACAAGATCACTCCCGGCATGCTGATGGCAGCACTCCGCCTCAATATCCCGGTCGTGTTCGTTTCGGGCGGGCCGATGGAAGCAGGCAAg                            >  2:4151117/7‑150 (MQ=255)
                                                                                                                         ccGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGa                           <  1:1547191/150‑1 (MQ=255)
                                                                                                                           gTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGAcc                         >  2:385826/1‑150 (MQ=255)
                                                                                                                             ggAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGACCaa                       <  1:2250593/150‑1 (MQ=255)
                                                                                                                             ggAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGACCaa                       >  1:1984610/1‑150 (MQ=255)
                                                                                                                              gAATATATGGTCAATGCCCATTGCGCCGACGCGATGGTCTGCATTTCCAACTGCGACAAGATCACTCCCGGCATGCTGATGGCAGCACTCCGCCTCAATATCCCGGTCGTGTTCGTTTCgggcgg                                               <  2:1279463/125‑1 (MQ=255)
                                                                                                                              gAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGGCAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGACCAAg                      >  1:2133114/1‑150 (MQ=255)
                                                                                                                                 tACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGaa                                                                          >  2:3592334/1‑95 (MQ=255)
                                                                                                                                 tACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGACCAAGCTg                   <  1:755326/150‑1 (MQ=255)
                                                                                                                                    aTGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGACCAAGCTGgcc                >  1:2358662/1‑150 (MQ=255)
                                                                                                                                    aTGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGACCAAGCTGgcc                <  1:1715594/150‑1 (MQ=255)
                                                                                                                                         cAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGACCAAGCTGgccagcca           >  1:5141088/1‑150 (MQ=255)
                                                                                                                                         cAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGACCAAGCTGgccagcca           >  2:3661201/1‑150 (MQ=255)
                                                                                                                                           aTGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGACCAAGCTGGCCAGCCACg         >  2:4496075/1‑150 (MQ=255)
                                                                                                                                             gcgcACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGGCAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGACCAAGCTGGCCAGCCACGGc       <  2:2133114/150‑1 (MQ=255)
                                                                                                                                             gcgcACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGACCAAGCTGGCCAGCCACGGc       >  1:4562676/1‑150 (MQ=255)
                                                                                                                                               gcACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGACCAAGCTGGCCAGCCACGGCCt     >  1:4697685/1‑150 (MQ=255)
                                                                                                                                               gcACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGACCAAGCTGGCCAGCCACGGCCt     >  1:4875473/1‑150 (MQ=255)
                                                                                                                                                cACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGACCAAGCTGGCCAGCCACGGCCTg    >  2:4429569/1‑150 (MQ=255)
                                                                                                                                                 aCTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGACCAAGCTGGCCAGCCACGGCCTgg   <  2:839135/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGACCAAGCTGGCCAGCCACGGCCTGGA  >  NC_002947/5851460‑5851755

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: