Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 554,621 C→G 5.5% 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 5.5% 1092.7 / 35.5 324T186T (ACC→ACG) rpsC30S ribosomal protein S3
Reads supporting (aligned to +/- strand):  ref base C (123/182);  new base G (8/10);  total (132/192)
Fisher's exact test for biased strand distribution p-value = 8.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 332 total aligned reads displayed.

TAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACCACAAGCACCAGCGCCTCGTAAAAAAGCTGC  >  NC_002947/554473‑554770
                                                                                                                                                    |                                                                                                                                                     
tAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGACAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACAccc                                                                                                                                                      <  2:6228354/150‑1 (MQ=255)
   tGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGc                                                                                                                                                   <  1:5151385/150‑1 (MQ=255)
   cattCCGTCGCGCCATGAAGCGTGCCGTTCAGTCGGCAATGCGTCTCGGTGCTGAAGGTATCCGTATCAACTGCTCCGGTCGTCTTGGCGGCGCTGAAATCGCACGTATGGCATGGTATCGCGAAGGTCGCGTTCCGCTTCACACGCTGc                                                                                                                                                   >  2:6453142/3‑150 (MQ=255)
   cattCCGTCGCGCCATGAAGCGTGCCGTTCAGTCGGCAATGCGTCTCGGTGCTGAAGGTATCCGTATCAACTGCTCCGGTCGTCTTGGCGGCGCTGAAATCGCACGTATGGAATGGTATCGCGAAGGTCGCGTTCCGCTTCACACGCTGc                                                                                                                                                   <  1:654689/148‑1 (MQ=255)
   cattCCGTCGCGCCATGAAGCGTGCCGTTCAGTCGGCAATGCGTCTCGGTGCTGAAGGTATCCGTATCAACTGCTCCGGTCGTCTTGGCGGCGCTGAAATCGCACGTATGGAATGGTATCGCGAAGGTCGCGTTCCGCTTCACACGCTGc                                                                                                                                                   <  1:5827995/148‑1 (MQ=255)
   cattCCGTCGCGCCATGAAGCGTGCCGTTCAGTCGGCAATGCGTCTCGGTGCTGAAGGTATCCGTATCAACTGCTCCGGTCGTCTTGGCGGCGCTGAAATCGCACGTATGGAATGGTATCGCGAAGGTCGCGTTCCGCTTCACACGCTGc                                                                                                                                                   <  2:7610194/148‑1 (MQ=255)
    gTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCg                                                                                                                                                  >  2:4231342/1‑150 (MQ=255)
       ccGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGc                                                                                                                                               <  2:5262654/150‑1 (MQ=255)
          tCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGa                                                                                                                                            <  1:4300693/150‑1 (MQ=255)
            gcgcCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGata                                                                                                                                          >  1:1328603/1‑150 (MQ=255)
              gcCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATc                                                                                                                                        >  1:9105871/1‑150 (MQ=255)
                cATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGa                                                                                                                                      <  1:2550704/150‑1 (MQ=255)
                cATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGa                                                                                                                                      >  1:1102822/1‑150 (MQ=255)
                  tGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACt                                                                                                                                    <  1:8746069/150‑1 (MQ=255)
                   gAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTa                                                                                                                                   <  1:1415964/150‑1 (MQ=255)
                    aaGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTac                                                                                                                                  >  1:7494066/1‑150 (MQ=255)
                     aGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTaca                                                                                                                                 <  2:4486912/150‑1 (MQ=255)
                       cGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTacaac                                                                                                                               <  2:4951818/150‑1 (MQ=255)
                          gCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACAcc                                                                                                                            <  2:2779292/150‑1 (MQ=255)
                            ggTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTa                                                                                                                          <  2:7727520/150‑1 (MQ=255)
                             gTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTAc                                                                                                                         <  2:6925085/150‑1 (MQ=255)
                               aCAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAGGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGa                                                                                                                       >  2:7335066/1‑150 (MQ=255)
                                       ccATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTcaca                                                                                                               <  1:2696661/150‑1 (MQ=255)
                                       ccATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTcaca                                                                                                               <  1:40521/150‑1 (MQ=255)
                                          tGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCAcacca                                                                                                            <  1:4971198/150‑1 (MQ=255)
                                          tGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCAcacca                                                                                                            <  2:4810277/150‑1 (MQ=255)
                                          tGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCAcacca                                                                                                            <  2:5700534/150‑1 (MQ=255)
                                             gTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACtt                                                                                                         >  1:4328776/1‑150 (MQ=255)
                                              tATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTa                                                                                                        >  2:1598174/1‑150 (MQ=255)
                                                ttGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCAcaccgtaagtg                                                                                                      <  2:5230005/150‑8 (MQ=255)
                                                ttGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACg                                                                                                      <  1:3759094/150‑1 (MQ=255)
                                                   gTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGgtg                                                                                                   <  1:7350100/150‑1 (MQ=255)
                                                          gggCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCgg                                                                                            <  2:4781050/150‑1 (MQ=255)
                                                           ggCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGgt                                                                                           <  2:3132413/150‑1 (MQ=255)
                                                           ggCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGgt                                                                                           <  1:8302562/150‑1 (MQ=255)
                                                           ggCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGAGTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGgt                                                                                           <  2:3291531/150‑1 (MQ=255)
                                                             cATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGgtgt                                                                                         <  2:2982202/150‑1 (MQ=255)
                                                              aTCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGtgtg                                                                                        >  1:7166270/1‑150 (MQ=255)
                                                                cAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGaa                                                                                      >  2:4354062/1‑150 (MQ=255)
                                                                   gATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGt                                                                                   >  1:2267773/1‑150 (MQ=255)
                                                                     tCCAGGTGAGCGGTCGTCTCGGCGGTGCTTAGATTGCTAGTACCGAGTGGTATCGCGAATGTAGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACTAAGCTCACACCACTTACTGTGTGATCGGTGTGAAGGttt                                                                                 <  1:5636007/150‑1 (MQ=255)
                                                                     tCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGttt                                                                                 <  1:5194309/150‑1 (MQ=255)
                                                                     tCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGttt                                                                                 <  1:6882994/150‑1 (MQ=255)
                                                                     tCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGttt                                                                                 <  1:5996812/150‑1 (MQ=255)
                                                                     tCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGttt                                                                                 <  1:7597949/150‑1 (MQ=255)
                                                                     tCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGttt                                                                                 <  2:1828765/150‑1 (MQ=255)
                                                                      ccAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTATAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTg                                                                                <  1:4251888/150‑1 (MQ=255)
                                                                      ccAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTg                                                                                >  1:6189610/1‑150 (MQ=255)
                                                                         ggTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGAt                                                                             <  1:2646213/150‑1 (MQ=255)
                                                                             aGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTc                                                                         <  1:6774021/150‑1 (MQ=255)
                                                                              ccGGTCGTCTTGGCGGCGCTGAAATCGCACGTATGGAATGGTATCGCGAAGGTCGCGTTCCGCTTCACACGCTGCGTGCCGATATCGACTACGGTACGGCAGAAGCGCAGACCGCTTACGGCATCTGCGGCGTGAAGGTTTGGGTATTCa                                                                        >  2:1071972/2‑150 (MQ=255)
                                                                                ggTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCaan                                                                      <  1:3500781/150‑2 (MQ=255)
                                                                                ggTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCaaa                                                                      <  2:8802395/150‑1 (MQ=255)
                                                                                ggTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCaaa                                                                      <  2:1496441/150‑1 (MQ=255)
                                                                                ggTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCaaa                                                                      <  2:3758630/150‑1 (MQ=255)
                                                                                 gtcgtcTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAg                                                                     <  2:5888810/150‑1 (MQ=255)
                                                                                  tcgtcTTGGCGGCGCTGAAATCGCACGTATGGAATGGTATCGCGAAGGTCGCGTTCCGCTTCACACGCTGCGTGCCGATATCGACTACGGTACGGCAGAAGCGCAGACCGCTTACGGCATCTGCGGCGTGAAGGTTTGGGTATTCAAggg                                                                    >  1:1864144/1‑150 (MQ=255)
                                                                                  tcgtcTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAgg                                                                    <  2:741386/150‑1 (MQ=255)
                                                                                     tctcGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGa                                                                 >  1:5537669/1‑150 (MQ=255)
                                                                                      cttggcggcGCTGAAATCGCACGTATGGAATGGTATCGCGAAGGTCGCGTTCCGCTTCACACGCTGCGTGCCGATATCGACTACGGTACGGCAGAAGCGCAGACCGCTTACGGCATCTGCGGCGTGAAGGTTTGGGTATTCAAGGGCGag                                                                <  2:2012829/147‑2 (MQ=255)
                                                                                      ctcGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGaa                                                                <  1:5370618/150‑1 (MQ=255)
                                                                                       tcGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAg                                                               >  1:1477188/1‑150 (MQ=255)
                                                                                         ggcggTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGtt                                                             >  1:7906680/1‑150 (MQ=255)
                                                                                          gcggTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTa                                                            <  2:9086213/150‑1 (MQ=255)
                                                                                            ggTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTAtt                                                          <  2:5211138/150‑1 (MQ=255)
                                                                                            ggTGCTGAGATTGCTCGTACCGAGTGGTATCGAGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTAtt                                                          <  1:7059019/150‑1 (MQ=255)
                                                                                             gTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATtg                                                         <  1:383306/150‑1 (MQ=255)
                                                                                             gTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATtg                                                         <  2:5582899/150‑1 (MQ=255)
                                                                                              tGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATtgg                                                        <  1:4512581/150‑1 (MQ=255)
                                                                                              tGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATtgg                                                        <  1:5240185/150‑1 (MQ=255)
                                                                                                 tGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATtggtgg                                                     <  2:2838710/150‑1 (MQ=255)
                                                                                                 tGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATtggtgg                                                     >  2:2212661/1‑150 (MQ=255)
                                                                                                  gagaTTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTggtggt                                                    <  2:8993774/150‑1 (MQ=255)
                                                                                                     aTTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGc                                                 >  1:6534846/1‑150 (MQ=255)
                                                                                                        gCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCaa                                              <  2:8149972/150‑1 (MQ=255)
                                                                                                          tCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCaaga                                            <  1:6799596/150‑1 (MQ=255)
                                                                                                              aCCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAgaagaa                                        <  1:7685022/150‑1 (MQ=255)
                                                                                                               ccGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAAc                                       <  1:4460200/150‑1 (MQ=255)
                                                                                                               ccGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAAc                                       <  1:3117918/150‑1 (MQ=255)
                                                                                                                cGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACt                                      >  1:2911467/1‑150 (MQ=255)
                                                                                                                    tGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGaaa                                  <  1:6687265/150‑1 (MQ=255)
                                                                                                                    tGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGaaa                                  <  1:8156513/150‑1 (MQ=255)
                                                                                                                      gTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAAcc                                >  1:8434841/1‑150 (MQ=255)
                                                                                                                          cgcgAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACCACaa                            <  2:7155757/150‑1 (MQ=255)
                                                                                                                          cgcgAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACCACaa                            <  2:6654929/150‑1 (MQ=255)
                                                                                                                           gcgAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACCACAAg                           >  1:3044101/1‑150 (MQ=255)
                                                                                                                            cgAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACCACAAGc                          >  1:958999/1‑150 (MQ=255)
                                                                                                                            cgAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACCACAAGc                          >  2:7060444/1‑150 (MQ=255)
                                                                                                                              aaGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACCACAAGCAc                        <  1:236365/150‑1 (MQ=255)
                                                                                                                                 gTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACCACAAGCACCAg                     <  2:1328603/150‑1 (MQ=255)
                                                                                                                                   cGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACCACAAGCACCAgcg                   >  1:5394636/1‑150 (MQ=255)
                                                                                                                                      gtgCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACCACAAGCACCAGCGCCt                <  2:606445/150‑1 (MQ=255)
                                                                                                                                       tgCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACCACAAGCACCAGCGCCTc               <  2:6570543/150‑1 (MQ=255)
                                                                                                                                       gatatctGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACCACAAGCACCAGCGCCTc               >  1:6157570/5‑150 (MQ=255)
                                                                                                                                          ctctGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACCACAAGCACCAGCGCCTCGTa            >  1:2480609/1‑150 (MQ=255)
                                                                                                                                          ctctGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACCACAAGCACCAGCGCCTCGTa            >  1:1754579/1‑150 (MQ=255)
                                                                                                                                               cacaCCCTGCGTGCCGATATCGACTACAACACCTACTAAGCTCACACCACTTACGGTGTGATCGGTGTTAATGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACCACAAGCACCAGCGCCTCGTaaaaaa       <  2:5544129/150‑1 (MQ=255)
                                                                                                                                                acacCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGaaa                                  >  2:2482650/1‑122 (MQ=255)
                                                                                                                                                acacCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACCACAAGCACCAGCGCCTCGTAAAAAAg      >  2:1143498/1‑150 (MQ=255)
                                                                                                                                                  acCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACCACAAGCACCAGCGCCTCGTAAAAAAgct    <  2:6793038/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACCACAAGCACCAGCGCCTCGTAAAAAAGCTGC  >  NC_002947/554473‑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: ‑