Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,076,526 C→A 82.9% D211Y (GAC→TAC)  speB ← agmatinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,076,5260CA82.9% 591.4 / 134.2 246D211Y (GAC→TAC) speBagmatinase
Reads supporting (aligned to +/- strand):  ref base C (20/22);  new base A (92/112);  total (112/134)
Fisher's exact test for biased strand distribution p-value = 8.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

Only 100 of 246 total aligned reads displayed.

GTACCAGTTTAATAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGT  >  CP009273/3076379‑3076672
                                                                                                                                                   |                                                                                                                                                  
gTACCAGTTTAATAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAAt                                                                                                                                                  >  1:2770835/1‑150 (MQ=255)
 tACCAGTTTAATAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATc                                                                                                                                                 <  1:1207703/150‑1 (MQ=255)
  aCCAGTTTAATAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATcc                                                                                                                                                >  2:2404214/1‑150 (MQ=255)
     aGTTTAATAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACg                                                                                                                                             >  2:2218280/1‑150 (MQ=255)
       tttAATAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCt                                                                                                                                           <  1:807730/150‑1 (MQ=255)
        ttAATAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTg                                                                                                                                          >  2:877240/1‑150 (MQ=255)
        ttAATAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTg                                                                                                                                          >  2:877116/1‑150 (MQ=255)
          aaTAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCg                                                                                                                                        <  1:368392/150‑1 (MQ=255)
                 cgATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTc                                                                                                                                 >  1:532636/1‑150 (MQ=255)
                  gATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCa                                                                                                                                >  2:3173438/1‑150 (MQ=255)
                   aTCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCAc                                                                                                                               <  2:3350380/150‑1 (MQ=255)
                     cGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCt                                                                                                                             <  1:304349/150‑1 (MQ=255)
                         ggTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATGTCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCa                                                                                                                         <  1:1836260/150‑1 (MQ=255)
                         ggTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCa                                                                                                                         <  1:4368523/150‑1 (MQ=255)
                           tCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGTTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAgg                                                                                                                       <  1:2708077/150‑1 (MQ=255)
                                  gccAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAg                                                                                                                <  2:534354/150‑1 (MQ=255)
                                  gccAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAg                                                                                                                <  2:1160623/150‑1 (MQ=255)
                                  gccAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAg                                                                                                                >  2:1877183/1‑150 (MQ=255)
                                   ccAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGc                                                                                                               >  2:388566/1‑150 (MQ=255)
                                   ccAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGc                                                                                                               >  2:1193433/1‑150 (MQ=255)
                                       tCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACgg                                                                                                           >  1:3213445/1‑150 (MQ=255)
                                        cACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGt                                                                                                          >  1:4517455/1‑150 (MQ=255)
                                        cACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGt                                                                                                          <  1:811763/150‑1 (MQ=255)
                                        cACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGt                                                                                                          >  2:4029442/1‑150 (MQ=255)
                                        cACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGt                                                                                                          <  1:4438389/150‑1 (MQ=255)
                                          cTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTaa                                                                                                        >  2:1516304/1‑150 (MQ=255)
                                            ggCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAg                                                                                                      <  1:769457/150‑1 (MQ=255)
                                             gCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGc                                                                                                     <  2:31261/150‑1 (MQ=255)
                                             gCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGc                                                                                                     <  1:2528084/150‑1 (MQ=255)
                                             gCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGc                                                                                                     <  2:1854371/150‑1 (MQ=255)
                                              cGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGcc                                                                                                    <  2:96651/150‑1 (MQ=255)
                                               gTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCg                                                                                                   <  1:1461084/150‑1 (MQ=255)
                                                   cGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGt                                                                                               <  2:2770835/150‑1 (MQ=255)
                                                   cGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGt                                                                                               >  1:2030807/1‑150 (MQ=255)
                                                    ggTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTc                                                                                              >  2:3667077/1‑150 (MQ=255)
                                                        ccTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTa                                                                                          >  1:2935167/1‑150 (MQ=255)
                                                        ccTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTa                                                                                          <  2:2500776/150‑1 (MQ=255)
                                                        ccTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTa                                                                                          <  1:1036307/150‑1 (MQ=255)
                                                           ggTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCa                                                                                       <  2:7203/150‑1 (MQ=255)
                                                               cAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACt                                                                                   >  2:2155914/1‑150 (MQ=255)
                                                               cAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACt                                                                                   <  2:4462419/150‑1 (MQ=255)
                                                                  aaGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCgg                                                                                <  1:1097275/150‑1 (MQ=255)
                                                                  aaGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCgg                                                                                >  1:2036070/1‑150 (MQ=255)
                                                                   aGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGt                                                                               >  1:2389615/1‑150 (MQ=255)
                                                                     cAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTAc                                                                             >  1:1884341/1‑150 (MQ=255)
                                                                     cAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTAc                                                                             >  2:2329324/1‑150 (MQ=255)
                                                                        gATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGaa                                                                          <  2:2658909/150‑1 (MQ=255)
                                                                         aTCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAAt                                                                         <  1:3354204/150‑1 (MQ=255)
                                                                           ccAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATAc                                                                       <  1:3418436/150‑1 (MQ=255)
                                                                            cAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATAcc                                                                      <  2:3110592/150‑1 (MQ=255)
                                                                              ggCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCaa                                                                    <  2:3467037/150‑1 (MQ=255)
                                                                               gCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAAt                                                                   <  1:976575/150‑1 (MQ=255)
                                                                                cAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATc                                                                  >  2:3310213/1‑150 (MQ=255)
                                                                                 aGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCt                                                                 >  2:3824711/1‑150 (MQ=255)
                                                                                 aGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCt                                                                 <  2:2590604/150‑1 (MQ=255)
                                                                                  gTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTg                                                                <  1:3580645/150‑1 (MQ=255)
                                                                                    cGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCa                                                              >  2:1151665/1‑150 (MQ=255)
                                                                                    cGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCa                                                              >  2:4329551/1‑150 (MQ=255)
                                                                                     tatatCAAAAGTCAGGTAATCCGGCATATCACCCACAATCTGTTTCACTTGGGAGATAACGTAATCCACGCTGCGATCGTTCACCTGTCTGTAGTCCAGCACGGTATAGCCGTTGTCTTTATCAAACTAGGTACGAATACCAATCTGCAc                                                             <  2:203951/149‑1 (MQ=255)
                                                                                     gATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCAc                                                             <  1:694791/150‑1 (MQ=255)
                                                                                     gATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCAc                                                             >  1:1909992/1‑150 (MQ=255)
                                                                                     gATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCAc                                                             >  2:657787/1‑150 (MQ=255)
                                                                                      atatCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACg                                                            >  1:1985517/1‑150 (MQ=255)
                                                                                      atatCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACg                                                            <  2:1301124/150‑1 (MQ=255)
                                                                                       tatCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACgg                                                           >  2:3024326/1‑150 (MQ=255)
                                                                                          cAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAAt                                                        <  2:1601177/150‑1 (MQ=255)
                                                                                            aaaGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGa                                                      >  1:2809589/1‑150 (MQ=255)
                                                                                             aaGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGAt                                                     >  1:3988340/1‑150 (MQ=255)
                                                                                              aGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGAtt                                                    <  1:2241594/150‑1 (MQ=255)
                                                                                                     tAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGACTGATTCGGGTCg                                             >  2:2102996/1‑150 (MQ=255)
                                                                                                      aaaCCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGa                                            >  2:4089998/1‑150 (MQ=255)
                                                                                                        aCCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATc                                          <  2:2935811/150‑1 (MQ=255)
                                                                                                          cGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAg                                        >  1:4243387/1‑150 (MQ=255)
                                                                                                           ggCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGa                                       >  2:2804864/1‑150 (MQ=255)
                                                                                                           ggCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGa                                       >  1:2310790/1‑150 (MQ=255)
                                                                                                              atatCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACct                                    <  2:1270830/150‑1 (MQ=255)
                                                                                                                  cACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACcttctt                                <  2:4159643/150‑1 (MQ=255)
                                                                                                                   aCCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCttt                               >  1:373752/1‑150 (MQ=255)
                                                                                                                   aCCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCttt                               <  2:458131/150‑1 (MQ=255)
                                                                                                                     ccACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCg                             >  2:2602999/1‑150 (MQ=255)
                                                                                                                     ccACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCg                             >  2:337827/1‑150 (MQ=255)
                                                                                                                       acaATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTACTTTCGgc                           >  1:635085/1‑150 (MQ=255)
                                                                                                                          aTCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGcgcg                        >  1:2216640/1‑150 (MQ=255)
                                                                                                                            cTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGt                      <  1:962993/150‑1 (MQ=255)
                                                                                                                            cTGTTTCACTTGGGCGATAACGTAATCCACACTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGt                      <  1:4040762/150‑1 (MQ=255)
                                                                                                                             tGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGta                     <  2:2000513/150‑1 (MQ=255)
                                                                                                                               tttCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGtata                   >  2:3489204/1‑150 (MQ=255)
                                                                                                                               tttCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGtata                   <  1:1513013/150‑1 (MQ=255)
                                                                                                                                 tCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGa                 <  1:2571823/150‑1 (MQ=255)
                                                                                                                                  cACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGaa                >  1:3070348/1‑150 (MQ=255)
                                                                                                                                  cACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGaa                >  1:3780482/1‑150 (MQ=255)
                                                                                                                                   aCTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAAc               >  1:4145446/1‑150 (MQ=255)
                                                                                                                                    cTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACa              <  1:1446379/150‑1 (MQ=255)
                                                                                                                                     ttGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACAt             <  1:1131060/150‑1 (MQ=255)
                                                                                                                                        ggCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGt          <  1:4353607/150‑1 (MQ=255)
                                                                                                                                        ggCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCGTTATCAAACTCGGTACGAATACCAATCTGCACGGACTGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTCTAGAACATAGt          >  2:4238712/1‑150 (MQ=255)
                                                                                                                                             tAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGt     >  1:2886815/1‑150 (MQ=255)
                                                                                                                                              aaCGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTg    >  1:2990963/1‑150 (MQ=255)
                                                                                                                                              aaCGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTg    <  2:3773050/150‑1 (MQ=255)
                                                                                                                                                cGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGt  <  2:956783/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GTACCAGTTTAATAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGT  >  CP009273/3076379‑3076672

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: ‑