Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,816,601 A→T 93.9% H281L (CAC→CTC)  spoT → bifunctional (p)ppGpp synthetase II/ guanosine‑3',5'‑bis pyrophosphate 3'‑pyrophosphohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,816,6010AT93.9% 1102.0 / 43.0 320H281L (CAC→CTC) spoTbifunctional (p)ppGpp synthetase II/ guanosine‑3',5'‑bis pyrophosphate 3'‑pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  ref base A (12/8);  new base T (147/153);  total (159/161)
Fisher's exact test for biased strand distribution p-value = 3.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 323 total aligned reads displayed.

CGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCAG  >  CP009273/3816453‑3816746
                                                                                                                                                    |                                                                                                                                                 
cgcgTCAGTGGTCGCTAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGTACATCTACGCTTTCCTCGTGATCGTCAATGATTCTGACACCTTTTATCGCGTGCTGGGCCAGATGcac                                                                                                                                                  <  1:775699/150‑1 (MQ=255)
   gTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGc                                                                                                                                               <  1:1006519/150‑1 (MQ=255)
   gTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGc                                                                                                                                               <  1:5082257/150‑1 (MQ=255)
    tCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGcc                                                                                                                                              <  2:2086584/150‑1 (MQ=255)
     cAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCt                                                                                                                                             >  1:2076694/1‑150 (MQ=255)
         ggTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTAc                                                                                                                                         >  2:2657816/1‑150 (MQ=255)
          gTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACa                                                                                                                                        >  1:382612/1‑150 (MQ=255)
           tCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACaa                                                                                                                                       <  2:1876578/150‑1 (MQ=255)
           tCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACaa                                                                                                                                       >  2:4753910/1‑150 (MQ=255)
           tCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGAGCGTCAATGATTCTGACACCTGGTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACac                                                                                                                                       >  2:4243694/1‑149 (MQ=255)
            cgcgAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAg                                                                                                                                      >  1:1363607/1‑150 (MQ=255)
              cgAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGcc                                                                                                                                    <  1:3451301/150‑1 (MQ=255)
              cgAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGcc                                                                                                                                    <  2:4368411/150‑1 (MQ=255)
               gagaAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCcg                                                                                                                                   <  1:1862067/150‑1 (MQ=255)
               gagaAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCcg                                                                                                                                   >  1:4782254/1‑150 (MQ=255)
               gagaAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCcg                                                                                                                                   >  2:1510313/1‑150 (MQ=255)
               gagaAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCcg                                                                                                                                   <  1:4995365/150‑1 (MQ=255)
                agaAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCcgc                                                                                                                                  >  2:2466831/1‑150 (MQ=255)
                agaAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCcgc                                                                                                                                  <  1:4328770/150‑1 (MQ=255)
                agaAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCcgc                                                                                                                                  >  2:1571645/1‑150 (MQ=255)
                   aGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTc                                                                                                                               >  2:406782/1‑150 (MQ=255)
                    gCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTcc                                                                                                                              <  2:2791374/150‑1 (MQ=255)
                    gCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTcc                                                                                                                              <  2:4975378/150‑1 (MQ=255)
                    gCATCTTTATTCGATTGACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTcc                                                                                                                              <  1:5016615/150‑1 (MQ=255)
                     cATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCg                                                                                                                             >  2:490977/1‑150 (MQ=255)
                     cATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCg                                                                                                                             >  2:1926623/1‑150 (MQ=255)
                     cATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCg                                                                                                                             >  1:3873746/1‑150 (MQ=255)
                        cTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGc                                                                                                                          >  2:3573348/1‑150 (MQ=255)
                         tttATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGcc                                                                                                                         >  2:2987631/1‑150 (MQ=255)
                         tttATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGcc                                                                                                                         >  2:4767673/1‑150 (MQ=255)
                         tttATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGcc                                                                                                                         <  1:220110/150‑1 (MQ=255)
                           tATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCcgc                                                                                                                       >  2:3185657/1‑150 (MQ=255)
                           tATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACAACTGTTATCGCGTGATGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCcgc                                                                                                                       >  2:588955/1‑150 (MQ=255)
                            aTTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCcgcg                                                                                                                      >  2:3940087/1‑150 (MQ=255)
                             ttCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGt                                                                                                                     <  1:3965928/150‑1 (MQ=255)
                              tCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTg                                                                                                                    >  2:3153715/1‑150 (MQ=255)
                               cGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGa                                                                                                                   >  1:4055001/1‑150 (MQ=255)
                                   ttACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGa                                                                                                               <  1:479925/150‑1 (MQ=255)
                                    tACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGAc                                                                                                              >  1:3643366/1‑150 (MQ=255)
                                     aCTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACt                                                                                                             <  2:1666810/150‑1 (MQ=255)
                                       tGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACtat                                                                                                           >  1:3708625/1‑150 (MQ=255)
                                       tGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACtat                                                                                                           >  1:4690808/1‑150 (MQ=255)
                                          aaaaTGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATc                                                                                                        >  1:4196885/1‑150 (MQ=255)
                                          aaaaTGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACGATATc                                                                                                        >  1:4243524/1‑150 (MQ=255)
                                            aaTGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGc                                                                                                      <  1:2549994/150‑1 (MQ=255)
                                                gTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCAtt                                                                                                  <  2:4991137/150‑1 (MQ=255)
                                                     cAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCaaa                                                                                             >  2:4660991/1‑150 (MQ=255)
                                                      aaaGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCaaaa                                                                                            >  1:2218747/1‑150 (MQ=255)
                                                         gagCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCg                                                                                         >  1:2398938/1‑150 (MQ=255)
                                                            cagcGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAAc                                                                                      <  2:2585493/150‑1 (MQ=255)
                                                             agcGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACg                                                                                     >  2:1011937/1‑150 (MQ=255)
                                                               cGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGc                                                                                   >  1:1883742/1‑150 (MQ=255)
                                                                 ttttCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTa                                                                                 >  1:2953051/1‑150 (MQ=255)
                                                                        tCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCt                                                                          >  2:2367345/1‑150 (MQ=255)
                                                                               tGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGcaca                                                                   >  1:376106/1‑150 (MQ=255)
                                                                                   cATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTc                                                               <  2:304275/150‑1 (MQ=255)
                                                                                   cATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTc                                                               >  1:2936871/1‑150 (MQ=255)
                                                                                    ttCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCg                                                              <  2:4110304/149‑1 (MQ=255)
                                                                                    aTCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCg                                                              <  1:4446679/150‑1 (MQ=255)
                                                                                        aCGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCgatga                                                          >  2:647367/1‑150 (MQ=255)
                                                                                        aCGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCgatga                                                          <  1:1640939/150‑1 (MQ=255)
                                                                                          gCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATc                                                        <  1:991362/150‑1 (MQ=255)
                                                                                            tttCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCgg                                                      <  2:4612795/150‑1 (MQ=255)
                                                                                             ttCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGc                                                     >  1:1831735/1‑150 (MQ=255)
                                                                                              tCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGcc                                                    >  2:3110408/1‑150 (MQ=255)
                                                                                               ccGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCTATGATCGGccc                                                   >  2:335781/1‑150 (MQ=255)
                                                                                               ccGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGccc                                                   <  1:9020/150‑1 (MQ=255)
                                                                                                cgcgTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTCTATCGCCATTCCAAAAGCGAACGGCTATCAGTATTTGCACACCTCGATGATCGGCCCg                                                  <  1:732168/150‑1 (MQ=255)
                                                                                                cgcgTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCTATGATCGGCCCg                                                  >  1:1686614/1‑150 (MQ=255)
                                                                                                cgcgTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCg                                                  <  2:443446/150‑1 (MQ=255)
                                                                                                 gcgTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGc                                                 >  2:1679540/1‑150 (MQ=255)
                                                                                                   gTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCAc                                               >  2:4389869/1‑150 (MQ=255)
                                                                                                    tGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACg                                              >  2:1362649/1‑150 (MQ=255)
                                                                                                     gATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACgg                                             <  1:490977/150‑1 (MQ=255)
                                                                                                         gTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGtgtg                                         <  1:1473328/150‑1 (MQ=255)
                                                                                                           cAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCTGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGATCGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGcc                                       <  1:4295121/150‑1 (MQ=255)
                                                                                                            aaTGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCTATGATCGGCCCGCACGGTGTGCCg                                      <  2:1184799/150‑1 (MQ=255)
                                                                                                            aaTGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCg                                      <  1:1209882/150‑1 (MQ=255)
                                                                                                            aaTGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCg                                      >  1:3541907/1‑150 (MQ=255)
                                                                                                            aaTGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCg                                      <  1:2250020/150‑1 (MQ=255)
                                                                                                             aTGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCgg                                     >  2:4560026/1‑150 (MQ=255)
                                                                                                                  tCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAg                                <  1:1249045/150‑1 (MQ=255)
                                                                                                                  tCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAg                                <  1:1251969/150‑1 (MQ=255)
                                                                                                                  tCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAg                                <  1:2747053/150‑1 (MQ=255)
                                                                                                                   cTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAgg                               <  2:2433758/150‑1 (MQ=255)
                                                                                                                   cTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAgg                               >  1:2884012/1‑150 (MQ=255)
                                                                                                                    tGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGAGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGt                              <  2:2376037/150‑1 (MQ=255)
                                                                                                                       cacCTGTTATCGCGTGCTGGGCCAGATGCACCGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATCTCGCCAGTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCa                           <  1:458277/150‑1 (MQ=255)
                                                                                                                          cTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGAt                        <  1:2581073/150‑1 (MQ=255)
                                                                                                                            gTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATcc                      <  2:601537/150‑1 (MQ=255)
                                                                                                                              tATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCTATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGt                    <  2:4684911/150‑1 (MQ=255)
                                                                                                                               aTCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTa                   >  2:2619639/1‑150 (MQ=255)
                                                                                                                                 cgcgTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTAcc                 <  2:779090/150‑1 (MQ=255)
                                                                                                                                 cgcgTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTAcc                 <  1:2455906/150‑1 (MQ=255)
                                                                                                                                     tGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAg             <  1:1740314/150‑1 (MQ=255)
                                                                                                                                     tGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAg             <  1:5020842/150‑1 (MQ=255)
                                                                                                                                          ggCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATg        <  2:3873746/150‑1 (MQ=255)
                                                                                                                                           gCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATgg       >  2:3625087/1‑150 (MQ=255)
                                                                                                                                            ccAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGa      <  2:3324042/150‑1 (MQ=255)
                                                                                                                                             cAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGAc     >  1:2312430/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
CGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCAG  >  CP009273/3816453‑3816746

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: