Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 554,621 C→G 6.0% T186T (ACC→ACG rpsC → 30S ribosomal protein S3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947554,6210CG 6.0% 610.8 / 21.3 184T186T (ACC→ACGrpsC30S ribosomal protein S3
Reads supporting (aligned to +/- strand):  ref base C (68/105);  new base G (4/7);  total (72/112)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.05e-01

Only 100 of 189 total aligned reads displayed.

GTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACCACAAGCACCAGCGCCTCGTAAAAAAGCTGC  >  NC_002947/554472‑554770
                                                                                                                                                     |                                                                                                                                                     
  aaTGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCt                                                                                                                                                     <  2:3491005/150‑1 (MQ=255)
    tGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCATACCCTGc                                                                                                                                                   >  2:1445433/1‑150 (MQ=255)
     gTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCg                                                                                                                                                  >  2:723423/1‑150 (MQ=255)
     gTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCg                                                                                                                                                  >  2:4689201/1‑150 (MQ=255)
     gTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCg                                                                                                                                                  <  2:2470090/150‑1 (MQ=255)
          gTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTCTCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCg                                                                                                                                             >  2:4365424/1‑150 (MQ=255)
           tCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGa                                                                                                                                            <  1:2144446/150‑1 (MQ=255)
                 cATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGa                                                                                                                                      <  2:4979082/150‑1 (MQ=255)
                  aTGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGAc                                                                                                                                     >  2:2326231/1‑150 (MQ=255)
                  aTGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGAc                                                                                                                                     >  1:3589883/1‑150 (MQ=255)
                  aTGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGAc                                                                                                                                     <  2:301408/150‑1 (MQ=255)
                  aTGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGAc                                                                                                                                     <  2:3111208/150‑1 (MQ=255)
                    gAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTa                                                                                                                                   <  2:3033872/150‑1 (MQ=255)
                     aaGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTac                                                                                                                                  >  1:1393556/1‑150 (MQ=255)
                         gTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTacaaca                                                                                                                              <  1:1012360/150‑1 (MQ=255)
                         gTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTacaaca                                                                                                                              <  1:225901/150‑1 (MQ=255)
                             ggTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTa                                                                                                                          >  1:4863801/1‑150 (MQ=255)
                              gTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTAc                                                                                                                         <  1:1603679/150‑1 (MQ=255)
                              gTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTAc                                                                                                                         >  1:2426357/1‑150 (MQ=255)
                                  aGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGt                                                                                                                                                 <  2:37482/122‑1 (MQ=255)
                                  aGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGt                                                                                                                                                 >  1:37482/1‑122 (MQ=255)
                                    aaCGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTCTCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCt                                                                                                                   >  2:3571398/1‑150 (MQ=255)
                                    aaCGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCt                                                                                                                   >  2:3867025/1‑150 (MQ=255)
                                      cGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCTTGCCGATATCGACTACAACACCTACGAAGCTca                                                                                                                 <  2:3336814/150‑1 (MQ=255)
                                      cGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTca                                                                                                                 <  2:4151082/150‑1 (MQ=255)
                                        ccATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTcaca                                                                                                               <  2:3130148/150‑1 (MQ=255)
                                        ccATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTcaca                                                                                                               <  2:4652729/150‑1 (MQ=255)
                                         cATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCacac                                                                                                              <  2:1301472/150‑1 (MQ=255)
                                              gTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACtt                                                                                                         <  2:1659554/150‑1 (MQ=255)
                                                       ccAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGa                                                                                                >  2:2723316/1‑150 (MQ=255)
                                                         aaGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATc                                                                                              <  1:3030023/150‑1 (MQ=255)
                                                          aGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCg                                                                                             <  2:99448/150‑1 (MQ=255)
                                                           gggCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCgg                                                                                            <  1:2726087/150‑1 (MQ=255)
                                                            ggCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGgt                                                                                           <  1:2085076/150‑1 (MQ=255)
                                                            ggCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGgt                                                                                           <  2:4842002/150‑1 (MQ=255)
                                                             gCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCggcg                                                                                          <  2:3194735/150‑3 (MQ=255)
                                                               aTCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGtgtg                                                                                        >  2:1916494/1‑150 (MQ=255)
                                                               aTCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGtgtg                                                                                        >  2:1117442/1‑150 (MQ=255)
                                                                tCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGa                                                                                       <  2:4504061/150‑1 (MQ=255)
                                                                tCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGa                                                                                       <  1:351729/150‑1 (MQ=255)
                                                                  aaGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAg                                                                                     <  2:3944313/150‑1 (MQ=255)
                                                                   aGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAgg                                                                                    >  2:188173/1‑150 (MQ=255)
                                                                      tCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGTTCGTGTGCCTCTGCACACCCTGCGTGCCGTTATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGttt                                                                                 <  2:3678635/150‑1 (MQ=255)
                                                                      tCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGttt                                                                                 <  1:2890724/150‑1 (MQ=255)
                                                                      tCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGttt                                                                                 <  1:3204718/150‑1 (MQ=255)
                                                                      tCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGttt                                                                                 <  2:2316424/150‑1 (MQ=255)
                                                                       ccAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTg                                                                                <  2:3557654/150‑1 (MQ=255)
                                                                         aGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGa                                                                              <  2:193830/150‑1 (MQ=255)
                                                                                cGGTCGTCTTGGCGGCGCTGAAATCGCACGTATGGAATGGTATCGCGAAGGTCGCGTTCCGCTTCACACGCTGCGTGCCGATATCGACTACGGTACGGCAGAAGCGCAGACCGCTTACGGCATCTGCGGCGTGAAGGTTTGGGTATTCaa                                                                       >  2:2993638/1‑150 (MQ=255)
                                                                                cGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCaa                                                                       <  2:1998485/150‑1 (MQ=255)
                                                                                 ggTCGTCTTGGCGGCGCTGAAATCGCACGTATGGAATGGTATCGCGAAGGTCGCGTTCCGCTTCACACGCTGCGTGCCGATATCGACTACGGTACGGCAGAAGCGCAGACCGCTTACGGCATCTGCGGCGTGAAGGTTTGGGTATTCaag                                                                      <  2:2287628/150‑2 (MQ=255)
                                                                                 ggTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCaaa                                                                      <  1:2701612/150‑1 (MQ=255)
                                                                                   tcgtcTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAgg                                                                    >  2:4305277/1‑150 (MQ=255)
                                                                                   tcgtcTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAgg                                                                    >  2:3951318/1‑150 (MQ=255)
                                                                                     gtcTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCg                                                                  <  1:3963661/150‑1 (MQ=255)
                                                                                      tctcGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGa                                                                 <  2:2267106/150‑1 (MQ=255)
                                                                                      tctcGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGa                                                                 <  2:2267249/150‑1 (MQ=255)
                                                                                       ctcGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTCCGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGaa                                                                <  2:836752/150‑1 (MQ=255)
                                                                                       ctcGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGaa                                                                <  2:3166057/150‑1 (MQ=255)
                                                                                        tcGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTa                                                                                                                                   >  2:3901725/1‑82 (MQ=255)
                                                                                           gcggcgCTGAAATCGCACGTATGGAATGGTATCGCGAAGGTCGCGTTCCGCTTCACACGCTGCGTGCCGATATCGACTACGGTACGGCAGAAGCGCAGACCGCTTACGGCATCTGCGGCGTGAAGGTTTGGGTATTCAAGGGCGagatcc                                                            <  2:24053/150‑6 (MQ=255)
                                                                                            cggCGCTGAAATCGCACGTATGGAATGGTATCGCGAAGGTCGCGTTCCGCTTCACACGCTGCGTGCCGATATCGACTACGGTACGGCAGAAGCGCAGACCGCTTACGGCATCTGCGGCGTGAAGGTTTGGGTATTCAAGGGCGagatcct                                                           <  2:4466956/150‑7 (MQ=255)
                                                                                             ggTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTAtt                                                          <  2:1722938/150‑1 (MQ=255)
                                                                                              gTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATtg                                                         <  1:2630505/150‑1 (MQ=255)
                                                                                              gTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCGACGAAGCTCACACCACTTCCGGTGTGATCGGTGTGAAGGTTTGGATCTTCAACGGCGAAGTTATtg                                                         <  2:564618/150‑1 (MQ=255)
                                                                                                  tGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGAACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATtggtgg                                                     >  2:789685/1‑150 (MQ=255)
                                                                                                   gagaTTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTggtggt                                                    >  1:593076/1‑150 (MQ=255)
                                                                                                      aTTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGc                                                 <  1:2864161/150‑1 (MQ=255)
                                                                                                           tCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCaaga                                            >  2:2431395/1‑150 (MQ=255)
                                                                                                                ccGAGTGGTCGCGCGAAGGTCGTGTGCCTCTGCACCCCCTGCGTGCCGAGATCGACTACAACACCTACGAAGCTCACACCACTTACGGGGTGATCGGTGTGAAGGGTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCACGAAGAAc                                       >  2:4576269/1‑150 (MQ=255)
                                                                                                                   aGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGa                                    <  1:1161200/150‑1 (MQ=255)
                                                                                                                      ggTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAAc                                 <  1:100444/150‑1 (MQ=255)
                                                                                                                          tCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACcaca                             >  1:2789411/1‑150 (MQ=255)
                                                                                                                            gcgAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACCACAAg                           >  2:2082618/1‑150 (MQ=255)
                                                                                                                            gcgAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACCACAAg                           >  1:4212197/1‑150 (MQ=255)
                                                                                                                              gAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACCACAAGCa                         >  1:3799763/1‑150 (MQ=255)
                                                                                                                              gAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTCACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACCACAAGCa                         >  1:1216877/1‑150 (MQ=255)
                                                                                                                                 ggTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACCACAAGCACCa                      >  1:5063600/1‑150 (MQ=255)
                                                                                                                                   tCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGgt                                                                                           <  1:17212/79‑1 (MQ=255)
                                                                                                                                   tCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGgt                                                                                           >  2:17212/1‑79 (MQ=255)
                                                                                                                                     gtgtGCCACTGCACACCATGAGTGCCTATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACCACAAGCACCAgcgc                  >  2:122519/1‑150 (MQ=255)
                                                                                                                                      tgtgCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACCACAAGCACCAGCGcc                 >  1:3021150/1‑150 (MQ=255)
                                                                                                                                      tgtgCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACCACAAGCACCAGCGcc                 >  2:2541856/1‑150 (MQ=255)
                                                                                                                                       gtgCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACCACAAGCACCAGCGCCt                <  1:3446098/150‑1 (MQ=255)
                                                                                                                                        tgCCTCTGCACACCCTGCGTGCCGATATCGACTACAACAACTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACCACAAGCACCAGCGCCTc               <  1:4911358/150‑1 (MQ=255)
                                                                                                                                          ccTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACCACAAGCACCAGCGCCTCGt             <  1:3951139/150‑1 (MQ=255)
                                                                                                                                           ctctGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAAcc                                <  2:932734/130‑1 (MQ=255)
                                                                                                                                           ctctGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAAcc                                >  1:932734/1‑130 (MQ=255)
                                                                                                                                           ctctGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACCACAAGCACCAGCGCCTCGTa            >  2:5031451/1‑150 (MQ=255)
                                                                                                                                           ctctGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACCACAAGCACCAGCGCCTCGTa            >  1:548960/1‑150 (MQ=255)
                                                                                                                                           ctctGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACCACAAGCACCAGCGCCTCGTa            >  1:4450452/1‑150 (MQ=255)
                                                                                                                                           ctctGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACCACAAGCACCAGCGCCTCGTa            >  2:5031231/1‑150 (MQ=255)
                                                                                                                                            tctGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACCACAAGCACCAGCGCCTCGTaa           <  2:1803902/150‑1 (MQ=255)
                                                                                                                                               gCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACCACAAGCACCAGCGCCTCGTaaaaa        >  2:653568/1‑150 (MQ=255)
                                                                                                                                               gCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACCACAAGCACCAGCGCCTCGTaaaaa        <  1:727764/150‑1 (MQ=255)
                                                                                                                                                cacaCCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACCACAAGCACCAGCGCCTCGTaaaaaa       <  1:1309998/150‑1 (MQ=255)
                                                                                                                                                   acCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACCACAAGCACCAGCGCCTCGTAAAAAAgct    <  1:1928389/150‑1 (MQ=255)
                                                                                                                                                   acCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACCACAAGCACCAGCGCCTCGTAAAAAAgct    <  2:557208/150‑1 (MQ=255)
                                                                                                                                                    cccTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACCACAAGCACCAGCGCCTCGTAAAAAAgctg   <  2:1706100/150‑1 (MQ=255)
                                                                                                                                                     ccTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACCACAAGCACCAGCGCCTCGTAAAAAAgctgc  <  1:4585211/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACCACAAGCACCAGCGCCTCGTAAAAAAGCTGC  >  NC_002947/554472‑554770

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: