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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,076,5260CA50.6% 47.1 / inf 306D211Y (GAC→TAC) speBagmatinase
Reads supporting (aligned to +/- strand):  ref base C (87/64);  new base A (69/86);  total (156/150)
Fisher's exact test for biased strand distribution p-value = 2.28e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01

Only 100 of 308 total aligned reads displayed.

CGTACCAGTTTAATAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCAA  >  CP009273/3076378‑3076675
                                                                                                                                                    |                                                                                                                                                     
 gTACCAGTTTAATAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCAt                                                                                                                                                     <  1:3232569/150‑1 (MQ=255)
    ccAGTTTAATAGCGCGTTCGGAGGTCAGGCCGCCTATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCTGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCa                                                                                                                                                  <  1:3289562/150‑1 (MQ=255)
    ccAGTTTAATAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCa                                                                                                                                                  >  1:10658/1‑150 (MQ=255)
      aGTTTACTAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACg                                                                                                                                                <  1:1446653/150‑1 (MQ=255)
      aGTTTAATAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACg                                                                                                                                                >  1:2643418/1‑150 (MQ=255)
      aGTTTAATAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACg                                                                                                                                                <  1:5203251/150‑1 (MQ=255)
       gTTTAATAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGc                                                                                                                                               >  2:3692389/1‑150 (MQ=255)
        tttAATAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACTATCTGTGTCACTTGGGCGATAACGTCATCCACGCt                                                                                                                                              >  2:3159922/1‑150 (MQ=255)
        tttAATAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCt                                                                                                                                              <  2:5056106/150‑1 (MQ=255)
        tttAATAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCt                                                                                                                                              >  2:4293739/1‑150 (MQ=255)
          tAATAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGc                                                                                                                                            <  1:340019/150‑1 (MQ=255)
              atcgcgATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATc                                                                                                                                        <  2:2987547/148‑1 (MQ=255)
                cgcgATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCCGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATCACGTAATCCACGCTGCGATCGt                                                                                                                                      >  1:4559288/1‑150 (MQ=255)
                   gATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCa                                                                                                                                   >  1:3091970/1‑150 (MQ=255)
                   gATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGGCAGGTACACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCa                                                                                                                                   >  1:2203020/1‑150 (MQ=255)
                      cGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCt                                                                                                                                >  2:3482036/1‑150 (MQ=255)
                        gAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTgg                                                                                                                              <  2:4239432/150‑1 (MQ=255)
                          ggTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCa                                                                                                                            <  2:4748022/150‑1 (MQ=255)
                           gTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAg                                                                                                                           <  1:1245349/150‑1 (MQ=255)
                             cAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGc                                                                                                                         <  1:3111235/150‑1 (MQ=255)
                                 ccgccAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTcc                                                                                                                     <  1:1960256/150‑1 (MQ=255)
                                 ccgccAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTcc                                                                                                                     <  1:3565757/150‑1 (MQ=255)
                                  cgccAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCa                                                                                                                    <  1:716774/150‑1 (MQ=255)
                                   gccAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAg                                                                                                                   <  1:775013/150‑1 (MQ=255)
                                   gccAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACATGGCAGGCGTCCAg                                                                                                                   >  2:380171/1‑150 (MQ=255)
                                      aaTCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCAc                                                                                                                >  2:1594295/1‑150 (MQ=255)
                                      aaTCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCAc                                                                                                                >  2:1594345/1‑150 (MQ=255)
                                            tGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTaaa                                                                                                          >  1:4916900/1‑150 (MQ=255)
                                             ggCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAg                                                                                                         >  1:2440812/1‑150 (MQ=255)
                                                 tACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGt                                                                                                     <  1:33746/150‑1 (MQ=255)
                                                 tACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGt                                                                                                     <  1:34042/150‑1 (MQ=255)
                                                    cGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGt                                                                                                  <  1:934254/150‑1 (MQ=255)
                                                     ggTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTc                                                                                                 <  2:5010570/150‑1 (MQ=255)
                                                     ggTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTc                                                                                                 <  2:846813/150‑1 (MQ=255)
                                                     ggTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTc                                                                                                 >  2:2179787/1‑150 (MQ=255)
                                                      gTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCt                                                                                                >  2:5202577/1‑150 (MQ=255)
                                                      gTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCt                                                                                                >  1:1314302/1‑150 (MQ=255)
                                                      gTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCt                                                                                                <  2:2705605/150‑1 (MQ=255)
                                                      gTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCt                                                                                                <  2:2705701/150‑1 (MQ=255)
                                                           tGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATc                                                                                           <  2:4439029/150‑1 (MQ=255)
                                                             gTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCaa                                                                                         >  1:1550659/1‑150 (MQ=255)
                                                                  aaaGCAGGATCCAGTCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCAAGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCg                                                                                    <  2:1939702/150‑1 (MQ=255)
                                                                  aaaGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACCATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCg                                                                                    <  2:2331066/150‑1 (MQ=255)
                                                                    aGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGt                                                                                  >  1:1116669/1‑150 (MQ=255)
                                                                     gCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTa                                                                                 <  2:759279/150‑1 (MQ=255)
                                                                      cAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTAc                                                                                <  1:2828301/150‑1 (MQ=255)
                                                                      cAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTAc                                                                                <  1:4129932/150‑1 (MQ=255)
                                                                       aGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACg                                                                               <  2:2817525/150‑1 (MQ=255)
                                                                       aGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACg                                                                               >  1:4992363/1‑150 (MQ=255)
                                                                       aGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACg                                                                               <  2:2552098/150‑1 (MQ=255)
                                                                          aTCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAAt                                                                            <  1:1649855/150‑1 (MQ=255)
                                                                           tCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATa                                                                           <  2:3067575/150‑1 (MQ=255)
                                                                           tCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATa                                                                           <  2:2260058/150‑1 (MQ=255)
                                                                              aGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCAGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCa                                                                        >  1:4460390/1‑150 (MQ=255)
                                                                               ggCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCaa                                                                       <  2:3322832/150‑1 (MQ=255)
                                                                                gCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAAt                                                                      <  2:1333894/150‑1 (MQ=255)
                                                                                gCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAAt                                                                      <  1:1403074/150‑1 (MQ=255)
                                                                                      gATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCAc                                                                <  1:2381758/150‑1 (MQ=255)
                                                                                       atatCAAAAGTCAGGTAAGCCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACg                                                               >  1:1196404/1‑150 (MQ=255)
                                                                                        tatCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGTACgg                                                              >  2:2017626/1‑150 (MQ=255)
                                                                                           cAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAAt                                                           <  1:3826929/150‑1 (MQ=255)
                                                                                           cAAAAGTCAGGTAAACCGGCATATCACCCACAACCTGATTCACTGGGCCGCTAACGTAATACACGCTGCGAACGTTCACCTGGCAGGCGACCAGCACGGTAAAGCCGTTGTATTTATCAATCTCGGTACGACTACCCATCTGCACGGacc                                                           >  2:4828456/1‑148 (MQ=255)
                                                                                            aaaaGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATg                                                          >  1:4343909/1‑150 (MQ=255)
                                                                                                gTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTc                                                      <  2:4503425/150‑1 (MQ=255)
                                                                                                  cAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCgg                                                    <  1:636539/150‑1 (MQ=255)
                                                                                                  cAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCgg                                                    >  2:806190/1‑150 (MQ=255)
                                                                                                  cAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCgg                                                    >  2:829408/1‑150 (MQ=255)
                                                                                                   aGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCggg                                                   <  2:572356/150‑1 (MQ=255)
                                                                                                     gaaaaCCCTCATATCACACACACTCAGTTTCACTAGCGCGATAACGTAATACAAGCTGCGATCGATCACCTGGCAGGCGTCCAGCACGGTAAAACCGTTGACTTTATAAAACTCGGACCCAATACCAATCTGCACGGAAAGAATCGGGTc                                                 >  2:4043452/3‑150 (MQ=255)
                                                                                                     gTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTc                                                 <  1:4815722/150‑1 (MQ=255)
                                                                                                      tAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCg                                                <  1:2718829/150‑1 (MQ=255)
                                                                                                             gCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGAc                                         <  1:3481970/150‑1 (MQ=255)
                                                                                                               atatCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACct                                       <  2:3036086/150‑1 (MQ=255)
                                                                                                                tatCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACctt                                      >  1:985733/1‑150 (MQ=255)
                                                                                                                  tCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACcttct                                    >  1:4094512/1‑150 (MQ=255)
                                                                                                                  tCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACcttct                                    >  1:597788/1‑150 (MQ=255)
                                                                                                                     cccACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTc                                 >  2:3999027/1‑150 (MQ=255)
                                                                                                                       cacaATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCgg                               >  1:3932997/1‑150 (MQ=255)
                                                                                                                       cacaATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCgg                               >  2:1719790/1‑150 (MQ=255)
                                                                                                                        acaATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGgc                              <  1:2752274/150‑1 (MQ=255)
                                                                                                                         caATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTAGgcg                             >  1:952014/1‑150 (MQ=255)
                                                                                                                          aaTCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGgcgc                            >  1:2563896/1‑150 (MQ=255)
                                                                                                                           aTCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGcgcg                           >  1:5135598/1‑150 (MQ=255)
                                                                                                                             cTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGt                         >  2:1200782/1‑150 (MQ=255)
                                                                                                                             cTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGt                         >  2:3108732/1‑150 (MQ=255)
                                                                                                                               gTTTCACTTTCGCGAAAACGTAATCCACGATGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGACGTTGTCTTTATAAAACACGGTACGAATAAAAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCTGCGCGGaat                       >  2:1108049/1‑150 (MQ=255)
                                                                                                                               gTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGtat                       <  2:1591688/150‑1 (MQ=255)
                                                                                                                                tttCACTTGGGAGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGtata                      <  2:726724/150‑1 (MQ=255)
                                                                                                                                 ttCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAg                     >  2:3585473/1‑150 (MQ=255)
                                                                                                                                 ttCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAg                     <  2:1435576/150‑1 (MQ=255)
                                                                                                                                   cACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGaa                   <  1:4000739/150‑1 (MQ=255)
                                                                                                                                    aCTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAAc                  >  2:4695538/1‑150 (MQ=255)
                                                                                                                                      ttGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACAt                <  2:4155615/150‑1 (MQ=255)
                                                                                                                                        gggCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAg              >  2:3857042/1‑150 (MQ=255)
                                                                                                                                           cGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGc           <  2:4872117/150‑1 (MQ=255)
                                                                                                                                           cGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGGCCAGCACGGTAAAGCCGGTGTCTTTATCACACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGc           <  2:2275834/150‑1 (MQ=255)
                                                                                                                                           cGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGc           >  2:4625121/1‑150 (MQ=255)
                                                                                                                                               aaCGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTg       <  1:1897274/150‑1 (MQ=255)
                                                                                                                                                aCGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTgg      <  2:5134028/150‑1 (MQ=255)
                                                                                                                                                aCGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTgg      <  1:1740696/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CGTACCAGTTTAATAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCAA  >  CP009273/3076378‑3076675

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑