Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,816,601 A→T 81.4% 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,6010AT81.4% 855.1 / 182.7 305H281L (CAC→CTC) spoTbifunctional (p)ppGpp synthetase II/ guanosine‑3',5'‑bis pyrophosphate 3'‑pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  ref base A (32/25);  new base T (154/94);  total (186/119)
Fisher's exact test for biased strand distribution p-value = 4.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 312 total aligned reads displayed.

CGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCAGAT  >  CP009273/3816455‑3816748
                                                                                                                                                  |                                                                                                                                                   
cgTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGctca                                                                                                                                                   <  2:2596886/149‑1 (MQ=255)
   cAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGcc                                                                                                                                                <  1:2601471/149‑1 (MQ=255)
        gTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTAc                                                                                                                                           <  2:1638540/149‑1 (MQ=255)
         tCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACa                                                                                                                                          >  2:4441696/1‑149 (MQ=255)
           gcgAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAg                                                                                                                                        >  2:877702/1‑149 (MQ=255)
              agaAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCcg                                                                                                                                     >  2:1248777/1‑149 (MQ=255)
              agaAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCcg                                                                                                                                     >  2:2527118/1‑149 (MQ=255)
                aaGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCcgcg                                                                                                                                   >  1:1492390/1‑149 (MQ=255)
                  gCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTc                                                                                                                                 >  1:4659006/1‑149 (MQ=255)
                  gCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTc                                                                                                                                 >  1:3824918/1‑149 (MQ=255)
                   cATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTcc                                                                                                                                <  1:780540/149‑1 (MQ=255)
                    aTCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCg                                                                                                                               <  2:3967458/149‑1 (MQ=255)
                    aTCTTTATTCGATTTACCGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCg                                                                                                                               >  1:3610629/1‑149 (MQ=255)
                     tCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCgg                                                                                                                              <  1:3596596/149‑1 (MQ=255)
                       tttATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTTCGGGc                                                                                                                            >  2:875836/1‑149 (MQ=255)
                       tttATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTTCGGGc                                                                                                                            >  2:2186133/1‑149 (MQ=255)
                       tttATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGc                                                                                                                            <  1:286811/149‑1 (MQ=255)
                        ttATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGcc                                                                                                                           <  2:240035/149‑1 (MQ=255)
                         tATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCcg                                                                                                                          >  1:993322/1‑149 (MQ=255)
                         tATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCcg                                                                                                                          >  1:2749088/1‑149 (MQ=255)
                          aTTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCcgc                                                                                                                         >  2:4403055/1‑149 (MQ=255)
                             cGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTg                                                                                                                      <  1:866636/149‑1 (MQ=255)
                                tttACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGaaa                                                                                                                   >  1:732620/1‑149 (MQ=255)
                                tttACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGaaa                                                                                                                   <  1:1778856/149‑1 (MQ=255)
                                tttACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGaaa                                                                                                                   <  2:4197002/149‑1 (MQ=255)
                                tttACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGAGCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGaaa                                                                                                                   >  2:755513/1‑149 (MQ=255)
                                   aCTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGAc                                                                                                                <  1:3100669/149‑1 (MQ=255)
                                    cTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGCTCATGGACATCTACGCTTTCCTCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACt                                                                                                               <  2:4960638/149‑1 (MQ=255)
                                    cTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACt                                                                                                               <  2:3716365/149‑1 (MQ=255)
                                    cTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACt                                                                                                               >  1:3663361/1‑149 (MQ=255)
                                     tGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACta                                                                                                              >  2:830104/1‑149 (MQ=255)
                                       cAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACtata                                                                                                            <  1:2370293/149‑1 (MQ=255)
                                       cAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACtata                                                                                                            <  2:4849196/149‑1 (MQ=255)
                                         aaaTGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATc                                                                                                          <  2:4287243/149‑1 (MQ=255)
                                         aaaTGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATc                                                                                                          >  2:5050213/1‑149 (MQ=255)
                                         aaaTGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATc                                                                                                          <  2:4701894/149‑1 (MQ=255)
                                           aTGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGc                                                                                                        >  2:3665070/1‑149 (MQ=255)
                                           aTGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGc                                                                                                        >  1:4732061/1‑149 (MQ=255)
                                              gTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCAt                                                                                                     <  2:3250427/149‑1 (MQ=255)
                                                gCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTc                                                                                                   >  1:3968602/1‑149 (MQ=255)
                                                gCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTc                                                                                                   <  1:3440017/149‑1 (MQ=255)
                                                 cTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTcc                                                                                                  <  2:3103133/149‑1 (MQ=255)
                                                   cAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCaa                                                                                                >  2:3491469/1‑149 (MQ=255)
                                                      agagCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAg                                                                                             >  1:4946349/1‑149 (MQ=255)
                                                        agcagcGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCg                                                                                           >  2:2753251/1‑149 (MQ=255)
                                                           agcGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAAc                                                                                        >  2:2079966/1‑149 (MQ=255)
                                                            gcGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACg                                                                                       >  1:401363/1‑149 (MQ=255)
                                                            gcGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACg                                                                                       <  2:4671150/149‑1 (MQ=255)
                                                            gcGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACg                                                                                       >  2:364724/1‑149 (MQ=255)
                                                                   cACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCa                                                                                <  1:2398152/149‑1 (MQ=255)
                                                                       cGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCt                                                                            >  1:3294185/1‑149 (MQ=255)
                                                                         aTCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCttt                                                                          >  1:2020732/1‑149 (MQ=255)
                                                                          tCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTg                                                                         >  1:1388863/1‑149 (MQ=255)
                                                                          tCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTg                                                                         >  1:1130336/1‑149 (MQ=255)
                                                                          tCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTg                                                                         >  1:4793521/1‑149 (MQ=255)
                                                                              ggACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCGGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGcaca                                                                     <  2:1925732/149‑1 (MQ=255)
                                                                                aCATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACAcc                                                                   >  1:1749255/1‑149 (MQ=255)
                                                                                aCATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACAcc                                                                   >  1:170784/1‑149 (MQ=255)
                                                                                aCATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACAcc                                                                   >  2:1932582/1‑149 (MQ=255)
                                                                                  aTCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTc                                                                 >  2:1017123/1‑149 (MQ=255)
                                                                                  aTCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTc                                                                 <  1:4560192/149‑1 (MQ=255)
                                                                                       cGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCgatga                                                            <  2:2405502/149‑1 (MQ=255)
                                                                                       cGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCgatga                                                            <  2:2667782/149‑1 (MQ=255)
                                                                                         cTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTCTCAGTCTTTGCACACCTCGATGATc                                                          >  1:1003691/1‑149 (MQ=255)
                                                                                         cTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATc                                                          >  2:383282/1‑149 (MQ=255)
                                                                                          tttCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCg                                                         >  2:5203503/1‑149 (MQ=255)
                                                                                            tCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGc                                                       >  2:1525321/1‑149 (MQ=255)
                                                                                            tCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGc                                                       >  2:866933/1‑149 (MQ=255)
                                                                                            tCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGc                                                       <  1:4102792/149‑1 (MQ=255)
                                                                                             ccGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATCATCGGcc                                                      <  1:2999758/149‑1 (MQ=255)
                                                                                               gcgTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCg                                                    >  2:1877064/1‑149 (MQ=255)
                                                                                               gcgTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCg                                                    >  2:46386/1‑149 (MQ=255)
                                                                                                 gTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCa                                                  >  2:4675574/1‑149 (MQ=255)
                                                                                                   gATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACg                                                >  1:1310102/1‑149 (MQ=255)
                                                                                                   gATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACg                                                >  2:3000360/1‑149 (MQ=255)
                                                                                                        tCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGtgtg                                           >  1:2036119/1‑149 (MQ=255)
                                                                                                        tCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGtgtg                                           >  1:519878/1‑149 (MQ=255)
                                                                                                         cAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGc                                          <  1:2146050/149‑1 (MQ=255)
                                                                                                         cAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGc                                          <  1:431042/149‑1 (MQ=255)
                                                                                                         cAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGc                                          >  2:247099/1‑149 (MQ=255)
                                                                                                          aaTGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGcc                                         >  2:3413943/1‑149 (MQ=255)
                                                                                                                  tGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAgg                                 <  2:2723863/149‑1 (MQ=255)
                                                                                                                   gACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGt                                >  2:1746914/1‑149 (MQ=255)
                                                                                                                      acCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCa                             >  1:4857724/1‑149 (MQ=255)
                                                                                                                       ccTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAg                            >  1:4430430/1‑149 (MQ=255)
                                                                                                                        cTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGa                           <  2:4055525/149‑1 (MQ=255)
                                                                                                                            tATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCg                       >  1:910372/1‑149 (MQ=255)
                                                                                                                             aTCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGt                      >  2:3350426/1‑149 (MQ=255)
                                                                                                                                   tGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGaa                <  1:3769760/149‑1 (MQ=255)
                                                                                                                                      tGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGat             <  2:375448/149‑1 (MQ=255)
                                                                                                                                       gggCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGata            <  1:2274675/149‑1 (MQ=255)
                                                                                                                                        ggCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGatat           <  2:1022336/149‑1 (MQ=255)
                                                                                                                                         gCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATg          >  1:2241247/1‑149 (MQ=255)
                                                                                                                                             gATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGAcc      <  2:4474328/149‑1 (MQ=255)
                                                                                                                                             gATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTATGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGACGATATGGAcc      >  1:2063601/1‑149 (MQ=255)
                                                                                                                                             gATGCTCAGCCTGTACAAGCCGCGTCCGCGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACCGCTATCAGTCTATGCACACCTCGATGATCGGCCCGCACGGAGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGAcc      >  2:4511807/1‑149 (MQ=255)
                                                                                                                                             gATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGAcc      >  1:3733276/1‑149 (MQ=255)
                                                                                                                                             gATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGAcc      <  2:328132/149‑1 (MQ=255)
                                                                                                                                                 cacaGCCTGTACAAGCCGCGTCCGGGCCGCGTTAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCAGAt  <  1:1227100/149‑1 (MQ=255)
                                                                                                                                                 cacaGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCAGAt  <  1:247099/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCAGAT  >  CP009273/3816455‑3816748

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: