Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,816,601 A→T 90.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,6010AT90.4% 1064.6 / 79.1 333H281L (CAC→CTC) spoTbifunctional (p)ppGpp synthetase II/ guanosine‑3',5'‑bis pyrophosphate 3'‑pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  ref base A (17/15);  new base T (164/137);  total (181/152)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

Only 100 of 333 total aligned reads displayed.

GTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCA  >  CP009273/3816456‑3816745
                                                                                                                                                 |                                                                                                                                                
gTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGc                                                                                                                                              <  2:1713840/150‑1 (MQ=255)
 tCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGcc                                                                                                                                             >  2:4279631/1‑150 (MQ=255)
       gTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACa                                                                                                                                       <  1:876362/150‑1 (MQ=255)
          gcgAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGc                                                                                                                                    >  1:4346153/1‑150 (MQ=255)
           cgAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGcc                                                                                                                                   <  1:1400992/150‑1 (MQ=255)
             agaAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCcgc                                                                                                                                 >  1:4402888/1‑150 (MQ=255)
                aGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTc                                                                                                                              >  2:1273225/1‑150 (MQ=255)
                  cATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCg                                                                                                                            <  1:2665698/150‑1 (MQ=255)
                     cTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGc                                                                                                                         >  1:3950958/1‑150 (MQ=255)
                     cTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGc                                                                                                                         >  1:3432074/1‑150 (MQ=255)
                       ttATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCTACAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCcg                                                                                                                       >  2:4397935/1‑150 (MQ=255)
                         aTTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCcgcg                                                                                                                     <  2:547110/150‑1 (MQ=255)
                          ttCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGt                                                                                                                    >  1:3511203/1‑150 (MQ=255)
                          ttCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGt                                                                                                                    <  2:1908216/150‑1 (MQ=255)
                           tCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTg                                                                                                                   <  2:2084464/150‑1 (MQ=255)
                               tttACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAg                                                                                                               <  1:1226307/150‑1 (MQ=255)
                                ttACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGa                                                                                                              <  2:1078716/150‑1 (MQ=255)
                                 tACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGAc                                                                                                             <  2:2137889/150‑1 (MQ=255)
                                   cTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACta                                                                                                           >  1:3367778/1‑150 (MQ=255)
                                     gCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACtata                                                                                                         <  1:931760/150‑1 (MQ=255)
                                     gCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACtata                                                                                                         >  2:3124825/1‑150 (MQ=255)
                                      cAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTatat                                                                                                        <  2:2975919/150‑1 (MQ=255)
                                      cAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTatat                                                                                                        >  1:2105640/1‑150 (MQ=255)
                                      cAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTatat                                                                                                        <  2:693654/150‑1 (MQ=255)
                                      cAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTatat                                                                                                        >  1:1157358/1‑150 (MQ=255)
                                       aaaaTGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATc                                                                                                       <  1:2363754/150‑1 (MQ=255)
                                       aaaaTGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATc                                                                                                       >  2:3158582/1‑150 (MQ=255)
                                       aaaaTGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATc                                                                                                       >  2:4542803/1‑150 (MQ=255)
                                       aaaaTGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATc                                                                                                       <  1:783370/150‑1 (MQ=255)
                                        aaaTGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCg                                                                                                      >  1:28396/1‑150 (MQ=255)
                                        aaaTGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCg                                                                                                      >  1:3296548/1‑150 (MQ=255)
                                         aaTGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTCCGCTTTCCGCGTGATCGTCAATGATTCTGAGACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGc                                                                                                     >  1:2931966/1‑150 (MQ=255)
                                          aTGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGcc                                                                                                    >  1:2059006/1‑150 (MQ=255)
                                           tGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCa                                                                                                   >  1:1307660/1‑150 (MQ=255)
                                           tGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCa                                                                                                   <  1:4336577/150‑1 (MQ=255)
                                            ggTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCAt                                                                                                  >  2:831997/1‑150 (MQ=255)
                                             gTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCAtt                                                                                                 >  1:1469568/1‑150 (MQ=255)
                                               gCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTcc                                                                                               <  2:3293906/150‑1 (MQ=255)
                                               gCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTcc                                                                                               >  2:741128/1‑150 (MQ=255)
                                                 tCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCaa                                                                                             >  2:2293362/1‑150 (MQ=255)
                                                 tCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCaa                                                                                             >  1:2510168/1‑150 (MQ=255)
                                                   aaaGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCaaaa                                                                                           <  1:3223670/150‑1 (MQ=255)
                                                     agagCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGc                                                                                         >  1:1458521/1‑150 (MQ=255)
                                                     agagCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGc                                                                                         >  1:1049733/1‑150 (MQ=255)
                                                      gagCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCg                                                                                        >  2:2694870/1‑150 (MQ=255)
                                                         cagcGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAAc                                                                                     <  2:2467780/150‑1 (MQ=255)
                                                          agcGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACg                                                                                    <  2:2398598/150‑1 (MQ=255)
                                                          agcGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACg                                                                                    >  1:4508247/1‑150 (MQ=255)
                                                           gcGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCTGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACgg                                                                                   <  1:531091/150‑1 (MQ=255)
                                                            cGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGc                                                                                  <  1:3392622/150‑1 (MQ=255)
                                                             gTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCt                                                                                 >  2:3342608/1‑150 (MQ=255)
                                                             gTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCt                                                                                 >  2:1065645/1‑150 (MQ=255)
                                                             gTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCt                                                                                 >  1:2658655/1‑150 (MQ=255)
                                                                 tCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCa                                                                             >  2:3264336/1‑150 (MQ=255)
                                                                 tCACTCGATCATGGACATCTACGCTATCCGAGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATTCTCAGCCTGTACAAGCCGCGTCCGGGCAGCGAGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCa                                                                             >  2:3655024/1‑150 (MQ=255)
                                                                   aCTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGt                                                                           >  1:2981924/1‑150 (MQ=255)
                                                                     tCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCt                                                                         >  2:3743912/1‑150 (MQ=255)
                                                                     tCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCt                                                                         <  1:1957649/150‑1 (MQ=255)
                                                                       gATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCttt                                                                       >  2:3548238/1‑150 (MQ=255)
                                                                        aTCATGGACCTCGACGCTGTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTAGATCGCCATTGCAAAAGCGAACGGCTATCAGTCTTTg                                                                      <  1:3143393/150‑1 (MQ=255)
                                                                          cATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGca                                                                    >  2:3237796/1‑150 (MQ=255)
                                                                           aTGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGcac                                                                   >  2:3968349/1‑150 (MQ=255)
                                                                            tGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGcaca                                                                  >  1:1389683/1‑150 (MQ=255)
                                                                             ggACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCacac                                                                 >  2:3833766/1‑150 (MQ=255)
                                                                              gACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACAcc                                                                <  2:4359494/150‑1 (MQ=255)
                                                                              gACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACAcc                                                                <  2:3652738/150‑1 (MQ=255)
                                                                               aCATCTTCGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCt                                                               <  2:2262900/150‑1 (MQ=255)
                                                                                cacctccGCTCTCCGCGACATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGCTGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTc                                                              <  2:2354261/144‑1 (MQ=255)
                                                                                cATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTc                                                              >  2:2224344/1‑150 (MQ=255)
                                                                                 aTCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTct                                                             >  1:3744939/1‑149 (MQ=255)
                                                                                 aTCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCg                                                             <  1:2251145/150‑1 (MQ=255)
                                                                                     aCGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCtatga                                                         >  1:4267511/1‑150 (MQ=255)
                                                                                       gCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATc                                                       <  1:1909491/150‑1 (MQ=255)
                                                                                        cTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATAGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCg                                                      >  2:1190298/1‑150 (MQ=255)
                                                                                               cgTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCTATGATCGGCCCGCa                                               <  2:1045667/150‑1 (MQ=255)
                                                                                                gTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCAc                                              >  2:4032242/1‑150 (MQ=255)
                                                                                                 tGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACg                                             >  2:3022796/1‑150 (MQ=255)
                                                                                                         aaTGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCTATGATCGGCCCGCACGGTGTGCCg                                     <  2:1511357/150‑1 (MQ=255)
                                                                                                           tGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGt                                   >  1:3880962/1‑150 (MQ=255)
                                                                                                            gATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCTATGATCGGCCCGCACGGTGTGCCGGtt                                  >  1:1124787/1‑150 (MQ=255)
                                                                                                            gATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCTATGATCGGCCCGCACGGTGTGCCGGtt                                  >  1:687954/1‑150 (MQ=255)
                                                                                                             aTTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCTATGATCGGCCCGCACGGTGTGCCGGTTg                                 >  1:1821628/1‑150 (MQ=255)
                                                                                                             aTTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCTATGATCGGCCCGCACGGTGTGCCGGTTg                                 >  1:202881/1‑150 (MQ=255)
                                                                                                              ttCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGa                                <  2:2931966/150‑1 (MQ=255)
                                                                                                               tCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCTATGATCGGCCCGCACGGTGTGCCGGTTGAg                               <  2:3781749/150‑1 (MQ=255)
                                                                                                                          ttATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCg                    >  2:987033/1‑150 (MQ=255)
                                                                                                                           tATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCTATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGt                   >  1:3657485/1‑150 (MQ=255)
                                                                                                                            aTCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTa                  >  1:2195752/1‑150 (MQ=255)
                                                                                                                             tCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGACCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCTATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTAc                 >  1:2672790/1‑150 (MQ=255)
                                                                                                                              cgcgTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTAcc                <  2:4243022/150‑1 (MQ=255)
                                                                                                                               gcgTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCTATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCg               <  1:637256/150‑1 (MQ=255)
                                                                                                                               gcgTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCg               <  1:1425514/150‑1 (MQ=255)
                                                                                                                                cgTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCTATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGa              <  1:1707587/150‑1 (MQ=255)
                                                                                                                                     tGGGCCAGATGCACAGCCTGTACAAGCCGCGGCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGata         <  2:3452990/150‑1 (MQ=255)
                                                                                                                                      gggCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCTATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGatat        >  1:4233527/1‑150 (MQ=255)
                                                                                                                                       ggCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATg       >  2:3397278/1‑150 (MQ=255)
                                                                                                                                          cAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCTATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGAc    <  1:921219/150‑1 (MQ=255)
                                                                                                                                           aGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGAcc   >  2:1013209/1‑150 (MQ=255)
                                                                                                                                            gATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCa  >  1:167224/1‑150 (MQ=255)
                                                                                                                                            gATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCa  >  1:167103/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCA  >  CP009273/3816456‑3816745

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: