Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,106,752 T→G 7.0% Q183P (CAG→CCG)  PP_1878 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,106,7520TG 7.0% 1043.6 / 11.5 404Q183P (CAG→CCG) PP_1878hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (178/179);  new base G (3/34);  total (189/215)
Fisher's exact test for biased strand distribution p-value = 3.40e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 407 total aligned reads displayed.

CTCGTCTTCTGCCGGGGCCATGCGGTCCAGTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCAAGAAGCCCTCATGGGCAGCGCGGGTCACCTGA  >  NC_002947/2106604‑2106901
                                                                                                                                                    |                                                                                                                                                     
 tCGTCTTCTGCCGGGGCCATGCGGTCCAGTACCTGGGCATTGGCCCGGGCCAACCGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCCGGCGCGCCTgg                                                                                                                                                     >  2:3422842/1‑150 (MQ=255)
 tCGTCTTCTGCCGGGGCCATGCGGTCCAGTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTgg                                                                                                                                                     >  2:4385205/1‑150 (MQ=255)
 tCGTCTTCTGCCGGGGCCATGCGGTCCAGTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTgg                                                                                                                                                     >  2:4602370/1‑150 (MQ=255)
  cGTCTTCTGCCGGGGCCATGCGGTCCAGTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGc                                                                                                                                                    >  2:5218942/1‑150 (MQ=255)
   gTCTTCTGCCGGGGCCATGCGGTCCAGTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGcc                                                                                                                                                   <  1:2933484/150‑1 (MQ=255)
     cttctGCCGTGGCCATGCGGTCCAGTACCTGGGCATTGGCCCGGGCAAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACAACTTCGGCTCCGGCTGATTGCGCAAATACTGCACGCCGCGCTCTAAGGCCAGGCGCGCCTGGCCgg                                                                                                                                                 <  1:6908132/150‑1 (MQ=255)
     cttctGCCGGGGCCATGCGGTCCAGTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCgg                                                                                                                                                 <  1:6087516/150‑1 (MQ=255)
      ttctGCCGGGGCCATGCGGTCCAGTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCggg                                                                                                                                                <  2:2550498/150‑1 (MQ=255)
      ttctGCCGGGGCCAAGCGGTCCAGTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCTCGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCggg                                                                                                                                                <  2:7212999/150‑1 (MQ=255)
       tctGCCGGGGCCATGCGGTCCAGTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGc                                                                                                                                               <  1:2644800/150‑1 (MQ=255)
           ccGGGGCCATGCGGTCCAGTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTAATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCCGGCGCGCCAGGCCgggcgggc                                                                                                                                           >  1:6391766/1‑150 (MQ=255)
             ggggCCATGCGGTCCAGTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCt                                                                                                                                         >  2:579458/1‑150 (MQ=255)
               ggcccggCGGTCCAGGCCCGGGGCATTGGCCCGGGCCACCAGTGTCTCACGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCttt                                                                                                                                       <  1:6639843/144‑1 (MQ=255)
                gcccggCGGTCCAGTACCTGGGCATTGGCCCGGGCCAACCGTTTCGCACGATACGTCAGCTGCTCTTCATACACCGTCGGCTCCGGCTGCTTGCGCAACTACTGCACGCCGCGCTCGACGGCCCGGCGCGCCGGGCCGGGCTGGCCtttt                                                                                                                                      <  1:7156417/145‑1 (MQ=255)
                gCCATGCGGTCCAGTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCtttt                                                                                                                                      <  1:2327241/150‑1 (MQ=255)
                gCCATGCGGTCCAGTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCtttt                                                                                                                                      <  1:6765773/150‑1 (MQ=255)
                gCCATGCGGTCCAGTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCtttt                                                                                                                                      <  1:1036221/150‑1 (MQ=255)
                 ccATGCGGTCCAGTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCttttt                                                                                                                                     <  1:3722426/150‑1 (MQ=255)
                   aTGCGGTCCAGTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGc                                                                                                                                   >  1:3025974/1‑150 (MQ=255)
                   aTGCGGTCCAGTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGc                                                                                                                                   <  2:1071890/150‑1 (MQ=255)
                    tGCGGTCCAGTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCa                                                                                                                                  <  2:7273427/150‑1 (MQ=255)
                        gTCCAGTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGc                                                                                                                              >  2:1138500/1‑150 (MQ=255)
                        gTCCAGTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGc                                                                                                                              >  2:4264045/1‑150 (MQ=255)
                        gTCCAGTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCATCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCAGGCCGGGCTGGCCTTTTTGAAGGGc                                                                                                                              >  1:1297873/1‑150 (MQ=255)
                          ccAGTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCt                                                                                                                            <  1:4862595/150‑1 (MQ=255)
                            aGTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGt                                                                                                                          <  2:2789761/150‑1 (MQ=255)
                             gTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGtt                                                                                                                         >  2:6418066/1‑150 (MQ=255)
                             gTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTGCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCGGGCCGGGCTGGCCTTTTTGCAGGGCCTGtt                                                                                                                         <  1:334492/150‑1 (MQ=255)
                                  tGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCg                                                                                                                    >  1:4623115/1‑150 (MQ=255)
                                   gggCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGa                                                                                                                   <  1:1197921/150‑1 (MQ=255)
                                   gggCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCACAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGa                                                                                                                   <  2:6590872/150‑1 (MQ=255)
                                     gCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAgg                                                                                                                 <  2:2822639/150‑1 (MQ=255)
                                     gCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAgg                                                                                                                 <  1:5614942/150‑1 (MQ=255)
                                       atgggCCCGGGCAAACAGTTTCTCAAGATACGTCAGCTGCTCTTCCGCCCCCTTCGGCTCCGGCTGCTTGCGCAAATACGGCACGCCGCGCTCGAAGGCCCGGCGCGCCGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGgtt                                                                                                               <  1:2477851/147‑1 (MQ=255)
                                         tGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGATCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGgttgt                                                                                                             <  2:2611540/150‑1 (MQ=255)
                                          ggCCCGGGCCAACCGTTTCTCACGATCCGTCAGCTGCTCTTCCTCCACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCCCGCCGCGCTCGAAGGCCAGGCGCGCCGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGgttgtt                                                                                                            <  1:4398734/150‑1 (MQ=255)
                                          ggCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGgttgtt                                                                                                            <  2:490511/150‑1 (MQ=255)
                                          ggCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGgttgtt                                                                                                            <  1:3535939/150‑1 (MQ=255)
                                           gCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGgttgtt                                                                                                            >  1:4072859/1‑149 (MQ=255)
                                           gCCCGGGCAAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACCCCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGttgttg                                                                                                           <  1:5139441/150‑1 (MQ=255)
                                             ccGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGACGGCCCGGCGCGCCGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTgaa                                                                                                         <  1:2724323/150‑1 (MQ=255)
                                             ccGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTgaa                                                                                                         <  1:3123687/150‑1 (MQ=255)
                                             ccGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTgaa                                                                                                         <  2:6887561/150‑1 (MQ=255)
                                             ccGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTACATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTgaa                                                                                                         <  1:5322654/150‑1 (MQ=255)
                                              cGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGgttgttttag                                                                                                        >  1:3447230/1‑146 (MQ=255)
                                                ggCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTgaagaa                                                                                                      >  1:2618447/1‑150 (MQ=255)
                                                     ccAGTTGCTCAAGATACGTCAGCTGCTCTTCATACCCCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCta                                                                                                 <  2:7129561/149‑1 (MQ=255)
                                                     aCAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCta                                                                                                 <  2:4707232/150‑1 (MQ=255)
                                                     aCAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGGTGTTGAAGAACTCta                                                                                                 >  1:1403743/1‑150 (MQ=255)
                                                         tttCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATg                                                                                             >  2:2088913/1‑150 (MQ=255)
                                                            ctcAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATgcag                                                                                          <  2:2388542/150‑1 (MQ=255)
                                                            cccAAGATACGTCTGCTGCTCTTCATACCCCTTCGGCTCGGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATgcag                                                                                          <  1:6986395/148‑1 (MQ=255)
                                                             tcAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTg                                                                                                                        >  2:2671179/1‑119 (MQ=255)
                                                             tcAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATgcagc                                                                                         >  2:3525317/1‑150 (MQ=255)
                                                              cAAGATCCTTCAGCTGCTCTTCATACACCTTCGGCTCCGTCTGCTTGCTCAAATACTGCACGCCGCGCTCGATGGCCCGGCGCGCCGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATgcagca                                                                                        <  1:3247395/150‑1 (MQ=255)
                                                              cAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTgaagaa                                                                                                      <  1:5552237/136‑1 (MQ=255)
                                                               aaGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATgcagcag                                                                                       >  2:6017125/1‑150 (MQ=255)
                                                                aGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATgcagcagc                                                                                      >  1:6499315/1‑150 (MQ=255)
                                                                 gATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATgcagcagca                                                                                     >  2:3874360/1‑150 (MQ=255)
                                                                      gTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAgg                                                                                >  1:4312770/1‑150 (MQ=255)
                                                                       tCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCttttt                                                                                                                                     <  2:1400394/96‑1 (MQ=255)
                                                                                    aTACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATTGCGGATTTc                                                                  >  1:1893478/1‑150 (MQ=255)
                                                                                    aTACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTc                                                                  >  2:5031560/1‑150 (MQ=255)
                                                                                     tACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTcc                                                                 <  2:3988633/150‑1 (MQ=255)
                                                                                     tACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTcc                                                                 <  2:6343691/150‑1 (MQ=255)
                                                                                     tACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTcc                                                                 <  1:5218942/150‑1 (MQ=255)
                                                                                      acacCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCt                                                                >  2:4051699/1‑150 (MQ=255)
                                                                                       cacCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCCGGATATGGCGGATTTCCtt                                                               >  1:6850450/1‑150 (MQ=255)
                                                                                       cacCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCtt                                                               <  2:2735753/150‑1 (MQ=255)
                                                                                        acCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGAAATGGCGGATTTCCTTg                                                              >  1:7349013/1‑150 (MQ=255)
                                                                                         ccTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGa                                                             <  2:841826/150‑1 (MQ=255)
                                                                                         ccTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGa                                                             <  2:6389938/150‑1 (MQ=255)
                                                                                          cTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGAc                                                            >  1:3698830/1‑150 (MQ=255)
                                                                                          cTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGAc                                                            >  2:535878/1‑150 (MQ=255)
                                                                                              ggCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCTGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCaa                                                        <  1:4212851/150‑1 (MQ=255)
                                                                                                   cGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCtt                                                   <  1:2088913/150‑1 (MQ=255)
                                                                                                      ctgctTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGc                                                <  2:5412610/150‑1 (MQ=255)
                                                                                                         ctTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTc                                             <  2:1788921/150‑1 (MQ=255)
                                                                                                          ttGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCg                                            >  2:6557220/1‑150 (MQ=255)
                                                                                                               cAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTTACCCAATGCTTTGCCTCGACAGc                                       >  1:451267/1‑150 (MQ=255)
                                                                                                                aaaTACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCg                                      >  2:2490315/1‑150 (MQ=255)
                                                                                                                  aTACTGCACGCCGCGCTCGAAGGCCCGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCggg                                    <  1:5829861/150‑1 (MQ=255)
                                                                                                                         aCGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCaagaag                             >  1:5500340/1‑150 (MQ=255)
                                                                                                                          cgccgcGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCAAGAAGc                            <  2:5551323/150‑1 (MQ=255)
                                                                                                                           gccgcGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCAAGAAGcc                           >  2:4241322/1‑150 (MQ=255)
                                                                                                                              gcgcTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATTGTGGATTTACTTGACCCAATGCTTGGCCTCGACAGCGGGCAAGAAGCCCTc                        >  1:3775823/1‑150 (MQ=255)
                                                                                                                              gcgcTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCAAGAAGCCCTc                        >  2:3345850/1‑150 (MQ=255)
                                                                                                                              gcgcTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCAAGAAGCCCTc                        <  2:284961/150‑1 (MQ=255)
                                                                                                                                gcTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCAAGAAGCCCTCAt                      <  1:1609181/150‑1 (MQ=255)
                                                                                                                                 cagcccggCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCAAGAAGCCCTCATg                     >  1:3635457/7‑150 (MQ=255)
                                                                                                                                  tCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCAAGAAGCCCTCATgg                    <  2:1639/150‑1 (MQ=255)
                                                                                                                                  tCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCAAGAAGCCCTCATgg                    >  2:5091331/1‑150 (MQ=255)
                                                                                                                                  gcGACGGCCCGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCAAGAAGCCCTCATgg                    <  1:890589/149‑1 (MQ=255)
                                                                                                                                   cGAAGGCCAGGCGCGCCGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCAAGAAGCCCTCATggg                   <  2:3025974/150‑1 (MQ=255)
                                                                                                                                       ggCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCAAGAAGCCCTCATGGGCAgc               <  1:4730413/150‑1 (MQ=255)
                                                                                                                                         ccAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCAAGAAGCCCTCATGGGCAgcgc             <  2:1271030/150‑1 (MQ=255)
                                                                                                                                          cAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCAAGAAGCCCTCATGGGCAGcgcg            >  2:5915196/1‑150 (MQ=255)
                                                                                                                                                 gcCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCAAGAAGCCCTCATGGGCAGCGCGGGTCAcc     <  1:4987243/150‑1 (MQ=255)
                                                                                                                                                    tGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCAAGAAGCCCTCATGGGCAGCGCGGGTCACCTGa  <  1:5918951/150‑1 (MQ=255)
                                                                                                                                                    tGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCAAGAAGCCCTCATGGGCAGCGCGGGTCACCTGa  <  1:6723227/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CTCGTCTTCTGCCGGGGCCATGCGGTCCAGTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCAAGAAGCCCTCATGGGCAGCGCGGGTCACCTGA  >  NC_002947/2106604‑2106901

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