Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,816,601 A→T 92.8% 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,6010AT92.8% 998.1 / 50.8 291H281L (CAC→CTC) spoTbifunctional (p)ppGpp synthetase II/ guanosine‑3',5'‑bis pyrophosphate 3'‑pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  ref base A (18/3);  new base T (164/106);  total (182/109)
Fisher's exact test for biased strand distribution p-value = 3.27e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.14e-01

Only 100 of 292 total aligned reads displayed.

GCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCAGA  >  CP009273/3816454‑3816747
                                                                                                                                                   |                                                                                                                                                  
gcgTCAGTGGTCGCGAGAAGCATCTTTATTCGATGTACTGCAAAAGGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGctca                                                                                                                                                  <  2:2491984/150‑1 (MQ=255)
 cgTCAGTGGTCGCGAGATGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAg                                                                                                                                                 >  2:3611582/1‑150 (MQ=255)
 cgTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAg                                                                                                                                                 >  1:2955905/1‑150 (MQ=255)
  gTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGc                                                                                                                                                >  1:577757/1‑150 (MQ=255)
   tCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGcc                                                                                                                                               <  1:1217560/150‑1 (MQ=255)
     aGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTg                                                                                                                                             <  2:4221503/150‑1 (MQ=255)
        ggTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTAc                                                                                                                                          <  1:1367352/150‑1 (MQ=255)
         gTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACa                                                                                                                                         >  2:829737/1‑150 (MQ=255)
           cgcgAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAg                                                                                                                                       >  2:1140173/1‑150 (MQ=255)
             cgAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGcc                                                                                                                                     <  2:3682750/150‑1 (MQ=255)
             cgAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGcc                                                                                                                                     >  2:1450983/1‑150 (MQ=255)
               agaAGCATCTTTATTCGATTTACTGCAAAATTGTGCTCAAAGAGCAGCGATTTCACTCTATCATGGACATCTAAGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCcgc                                                                                                                                   >  2:518756/1‑150 (MQ=255)
                 aaGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGt                                                                                                                                 >  1:1789131/1‑150 (MQ=255)
                   gCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTcc                                                                                                                               <  1:2826714/150‑1 (MQ=255)
                   gCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTcc                                                                                                                               <  1:3220823/150‑1 (MQ=255)
                    cATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCg                                                                                                                              >  2:75885/1‑150 (MQ=255)
                     aTCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCgg                                                                                                                             >  2:594575/1‑150 (MQ=255)
                       cTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGc                                                                                                                           >  1:358441/1‑150 (MQ=255)
                               gATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGaa                                                                                                                   <  2:1924256/150‑1 (MQ=255)
                               gATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGaa                                                                                                                   >  1:4313366/1‑150 (MQ=255)
                                 tttACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAg                                                                                                                 >  2:1814523/1‑150 (MQ=255)
                                   tACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGAc                                                                                                               <  2:1510531/150‑1 (MQ=255)
                                   tACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGAc                                                                                                               <  2:2721545/150‑1 (MQ=255)
                                     cTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACtc                                                                                                             >  2:3039455/1‑149 (MQ=255)
                                      tGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACtat                                                                                                            >  1:2343132/1‑150 (MQ=255)
                                        cAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTatat                                                                                                          <  2:3512369/150‑1 (MQ=255)
                                          aaaTGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCg                                                                                                        <  1:4359745/150‑1 (MQ=255)
                                            aTGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGcc                                                                                                      <  1:174225/150‑1 (MQ=255)
                                              ggTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCAt                                                                                                    >  2:1405642/1‑150 (MQ=255)
                                               gTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCAtt                                                                                                   >  2:2183356/1‑150 (MQ=255)
                                                tGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCGGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTc                                                                                                  >  1:226985/1‑150 (MQ=255)
                                                     aaaGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCaaaa                                                                                             >  1:1950339/1‑150 (MQ=255)
                                                     aaaGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCaaaa                                                                                             >  1:2822811/1‑150 (MQ=255)
                                                     aaaGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCaaaa                                                                                             >  2:2722768/1‑150 (MQ=255)
                                                      aaGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAg                                                                                            <  1:265481/150‑1 (MQ=255)
                                                       agagCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGc                                                                                           >  2:1254877/1‑150 (MQ=255)
                                                        gagCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCg                                                                                          <  2:647548/150‑1 (MQ=255)
                                                         agcagcGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGa                                                                                         >  2:1324848/1‑150 (MQ=255)
                                                         agcagcGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGa                                                                                         <  2:1936539/150‑1 (MQ=255)
                                                         agcagcGTTTTCACTCGATCATGGACATCAACGCTTTCCGCGTGAACGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCAGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGa                                                                                         >  2:1552845/1‑150 (MQ=255)
                                                                ttttCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTa                                                                                  <  1:3955884/150‑1 (MQ=255)
                                                                    cACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAg                                                                              >  2:1321695/1‑150 (MQ=255)
                                                                     aCTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGt                                                                             >  1:3413359/1‑150 (MQ=255)
                                                                         gATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCttt                                                                         >  1:4171557/1‑150 (MQ=255)
                                                                           tCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTCTCAGTCTTTGc                                                                       >  2:2600212/1‑150 (MQ=255)
                                                                            cATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGca                                                                      >  1:1504334/1‑150 (MQ=255)
                                                                              tGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGcaca                                                                    >  2:1069777/1‑150 (MQ=255)
                                                                                 aCATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCt                                                                 >  1:96126/1‑150 (MQ=255)
                                                                                 aCATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCt                                                                 >  2:2246608/1‑150 (MQ=255)
                                                                                 aCATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCt                                                                 >  1:1675307/1‑150 (MQ=255)
                                                                                  cATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTc                                                                <  2:2649612/150‑1 (MQ=255)
                                                                                   aTCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCg                                                               >  1:3925891/1‑150 (MQ=255)
                                                                                   aTCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCg                                                               >  2:3065825/1‑150 (MQ=255)
                                                                                   aTCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCg                                                               <  2:1789131/150‑1 (MQ=255)
                                                                                    tCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCga                                                              <  1:3323200/150‑1 (MQ=255)
                                                                                     cTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCgat                                                             <  1:51639/150‑1 (MQ=255)
                                                                                     cTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCgat                                                             <  1:3720707/150‑1 (MQ=255)
                                                                                      tACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAACGCGAACGGCTATCAGTCTTTGCACACCTCgatg                                                            <  1:2909865/150‑1 (MQ=255)
                                                                                      tACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCgatg                                                            >  2:3149657/1‑150 (MQ=255)
                                                                                      tACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCgatg                                                            >  2:963228/1‑150 (MQ=255)
                                                                                       aCGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCgatga                                                           <  2:1748672/150‑1 (MQ=255)
                                                                                        cGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCgatgat                                                          <  1:1096301/150‑1 (MQ=255)
                                                                                          cTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCg                                                        >  2:3887445/1‑150 (MQ=255)
                                                                                          cTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCg                                                        <  1:772577/150‑1 (MQ=255)
                                                                                          cTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCg                                                        >  1:706583/1‑150 (MQ=255)
                                                                                            ttCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGc                                                      <  2:860112/150‑1 (MQ=255)
                                                                                              ccGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGGTGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGGCTTTGCACACCTCGATGATCGGccc                                                    >  1:4109969/1‑150 (MQ=255)
                                                                                                  gTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCACAAGCGAACGGCTCTCAGTCGTTGCACACCTCGATGATCGGCCCGCAc                                                >  2:1209624/1‑150 (MQ=255)
                                                                                                   tGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACg                                               <  2:407147/150‑1 (MQ=255)
                                                                                                          cAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGcc                                        <  1:3086837/150‑1 (MQ=255)
                                                                                                             tGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGt                                     <  1:3267169/150‑1 (MQ=255)
                                                                                                             tGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGt                                     >  1:480591/1‑150 (MQ=255)
                                                                                                             tGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGt                                     <  1:2068126/150‑1 (MQ=255)
                                                                                                                ttCTGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGa                                  <  2:196706/150‑1 (MQ=255)
                                                                                                                   tGACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGt                               >  1:2307578/1‑150 (MQ=255)
                                                                                                                    gACACCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTc                              >  2:2442182/1‑150 (MQ=255)
                                                                                                                      cacCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCa                            >  2:3268337/1‑150 (MQ=255)
                                                                                                                      cacCTGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCa                            <  2:1867672/150‑1 (MQ=255)
                                                                                                                       acCTGTTATCGCGTGCTGGGCCAGATGATCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAg                           <  2:2979151/150‑1 (MQ=255)
                                                                                                                          tGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATc                        <  1:3617506/150‑1 (MQ=255)
                                                                                                                          tGTTATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATc                        <  1:1424914/150‑1 (MQ=255)
                                                                                                                            ttATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCg                      >  2:3895246/1‑150 (MQ=255)
                                                                                                                             tATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGt                     >  2:912770/1‑150 (MQ=255)
                                                                                                                             tATCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGt                     >  2:1468570/1‑150 (MQ=255)
                                                                                                                             tATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGt                     >  2:1615090/1‑150 (MQ=255)
                                                                                                                             tATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGt                     >  2:3803545/1‑150 (MQ=255)
                                                                                                                              ctCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGACCGGCTATCAGTCTTTGCACACCTCGATGAGCGGCCCGCACGGTGTGCCGGTGGAGGTCCAGATCCGTa                    >  2:2937345/2‑150 (MQ=255)
                                                                                                                              aTCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAACGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTa                    >  2:3418894/1‑150 (MQ=255)
                                                                                                                              aTCGCGTGCTGGGCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTa                    <  1:3132622/150‑1 (MQ=255)
                                                                                                                                cgcgTGCTGGGCCAGATGCTCAGCCTGTCCAAGCCGCGTCCGGGCCGCGTGAAAGACTATAGCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTAcc                  <  1:1020776/150‑1 (MQ=255)
                                                                                                                                       tGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGata           >  2:2301524/1‑150 (MQ=255)
                                                                                                                                          gCCAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATgg        <  1:2722768/150‑1 (MQ=255)
                                                                                                                                            cAGATGCTCAGCCTGTACACGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGAACCGAAGATATGGAc      <  1:101568/150‑1 (MQ=255)
                                                                                                                                            cAGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGAc      <  2:2095787/150‑1 (MQ=255)
                                                                                                                                             aGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGAcc     <  1:2140140/150‑1 (MQ=255)
                                                                                                                                             aGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGAcc     >  1:737323/1‑150 (MQ=255)
                                                                                                                                             aGATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGAcc     >  1:922518/1‑150 (MQ=255)
                                                                                                                                              gATGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCa    <  2:4239558/150‑1 (MQ=255)
                                                                                                                                               aTGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCAg   >  1:1975110/1‑150 (MQ=255)
                                                                                                                                                tGCTCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCAGa  >  2:3119032/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCAGA  >  CP009273/3816454‑3816747

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: