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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,076,5260CA90.9% 976.0 / 86.3 329D211Y (GAC→TAC) speBagmatinase
Reads supporting (aligned to +/- strand):  ref base C (15/15);  new base A (139/160);  total (154/175)
Fisher's exact test for biased strand distribution p-value = 8.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.87e-01

Only 100 of 329 total aligned reads displayed.

GTACCAGTTTAATAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTC  >  CP009273/3076379‑3076673
                                                                                                                                                   |                                                                                                                                                   
 tACCAGTTTAATAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATc                                                                                                                                                  >  1:2637179/1‑150 (MQ=255)
 tACCAGTTTAATAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATc                                                                                                                                                  >  1:1047187/1‑150 (MQ=255)
 tACCAGTTTAATAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATc                                                                                                                                                  <  2:227573/150‑1 (MQ=255)
   ccAGTTTAATAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCa                                                                                                                                                <  2:3962218/150‑1 (MQ=255)
   ccAGTTTAATAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCa                                                                                                                                                <  2:1523524/150‑1 (MQ=255)
    cAGTTTAATAGCGCGATCGGAGGTCAGGCCGCCACTCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCAc                                                                                                                                               <  2:1863830/150‑1 (MQ=255)
           aTAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATCTCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGa                                                                                                                                        >  2:2497228/1‑150 (MQ=255)
           aTAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGa                                                                                                                                        >  2:618064/1‑150 (MQ=255)
            tAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGAt                                                                                                                                       >  2:304492/1‑150 (MQ=255)
              gcgcGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCg                                                                                                                                     <  1:823030/150‑1 (MQ=255)
                  gATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCa                                                                                                                                 >  2:2178685/1‑150 (MQ=255)
                  gATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCa                                                                                                                                 <  2:1779874/150‑1 (MQ=255)
                    tCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCAcc                                                                                                                               >  2:1420576/1‑150 (MQ=255)
                     cGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCt                                                                                                                              <  1:4349099/150‑1 (MQ=255)
                     cGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCt                                                                                                                              >  1:4220139/1‑150 (MQ=255)
                      ggAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTg                                                                                                                             <  1:3604785/150‑1 (MQ=255)
                      ggAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATACCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTg                                                                                                                             <  1:2888273/150‑1 (MQ=255)
                        aGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGc                                                                                                                           <  1:1524775/150‑1 (MQ=255)
                            cAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGc                                                                                                                       >  2:3526181/1‑150 (MQ=255)
                            cAGGCCGCCAATCACTGGCGTACCGGGGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGc                                                                                                                       >  2:3169875/1‑150 (MQ=255)
                             aGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCg                                                                                                                      >  1:48868/1‑150 (MQ=255)
                             aGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCg                                                                                                                      >  1:436606/1‑150 (MQ=255)
                              ggCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGt                                                                                                                     >  1:3998226/1‑150 (MQ=255)
                               gccgccAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTc                                                                                                                    <  2:380403/150‑1 (MQ=255)
                                ccgccAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTcc                                                                                                                   <  1:3960166/150‑1 (MQ=255)
                                ccgccAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTcc                                                                                                                   <  2:4286928/150‑1 (MQ=255)
                                  gccTATCTCTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAg                                                                                                                 <  1:4337130/150‑1 (MQ=255)
                                  gccAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAg                                                                                                                 >  1:728371/1‑150 (MQ=255)
                                   ccAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGc                                                                                                                >  2:4494954/1‑150 (MQ=255)
                                   ccAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGc                                                                                                                >  2:2091430/1‑150 (MQ=255)
                                        cACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGAGATCAATAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACGTGGGCGTTTAAGTAATCCACGCTGCGATCGTTCACATGGCAGGCGTCCAGCACgga                                                                                                           >  2:2595368/1‑149 (MQ=255)
                                           tGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTaaa                                                                                                        <  1:874808/150‑1 (MQ=255)
                                              cGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGcc                                                                                                     >  1:2091597/1‑150 (MQ=255)
                                              cGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGcc                                                                                                     >  1:3369878/1‑150 (MQ=255)
                                              cGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGcc                                                                                                     <  2:3007089/150‑1 (MQ=255)
                                               gTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCg                                                                                                    <  2:3011981/150‑1 (MQ=255)
                                                tACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGt                                                                                                   >  2:3978813/1‑150 (MQ=255)
                                                   cGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGt                                                                                                <  2:369267/150‑1 (MQ=255)
                                                   cGGAGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGt                                                                                                >  2:2056080/1‑150 (MQ=255)
                                                     gTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCt                                                                                              <  1:3282763/150‑1 (MQ=255)
                                                      tGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCtt                                                                                             >  1:3886660/1‑150 (MQ=255)
                                                      tGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACAGTAAAGCCGTTGTCtt                                                                                             >  2:2893018/1‑150 (MQ=255)
                                                       gCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTaaa                                                                                                        >  1:3468110/1‑138 (MQ=255)
                                                        ccTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTa                                                                                           >  1:4507134/1‑150 (MQ=255)
                                                        ccTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTa                                                                                           <  1:2239323/150‑1 (MQ=255)
                                                          tGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATc                                                                                         <  2:2054793/150‑1 (MQ=255)
                                                          tGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATc                                                                                         <  2:2112756/150‑1 (MQ=255)
                                                           ggTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCa                                                                                        <  1:3272945/150‑1 (MQ=255)
                                                                aaaaGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTACACCGGCATATCAACCACACTCTGTTTCACTTGGGCGATAACTTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACTGTAAAGCCGTTGTCTTTATCAAACTc                                                                                   <  2:2743986/150‑1 (MQ=255)
                                                                aaaaGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTc                                                                                   >  2:2432785/1‑150 (MQ=255)
                                                                aaaaGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTc                                                                                   >  2:1052226/1‑150 (MQ=255)
                                                                aaaaGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTc                                                                                   >  1:4079818/1‑150 (MQ=255)
                                                                aaaaGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTc                                                                                   >  2:3759010/1‑150 (MQ=255)
                                                                  aaGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCgg                                                                                 <  1:2737844/150‑1 (MQ=255)
                                                                     cAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTAc                                                                              <  1:3817215/150‑1 (MQ=255)
                                                                      aGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACg                                                                             >  1:2499978/1‑150 (MQ=255)
                                                                          tCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATa                                                                         <  2:2844102/150‑1 (MQ=255)
                                                                          tCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATa                                                                         <  1:4132557/150‑1 (MQ=255)
                                                                             aGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCa                                                                      <  1:304149/150‑1 (MQ=255)
                                                                             aGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCa                                                                      >  1:3178524/1‑150 (MQ=255)
                                                                                 aGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCt                                                                  <  1:602116/150‑1 (MQ=255)
                                                                                 aGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCt                                                                  >  2:3885081/1‑150 (MQ=255)
                                                                                  gTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTg                                                                 <  1:565866/150‑1 (MQ=255)
                                                                                   tCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGc                                                                >  1:3806037/1‑150 (MQ=255)
                                                                                   tCGATAACAAAACTCAGGTAAACCGCCATATCACCCACAAACTGTATCACTTGGGCGATAACGTAATCCACGCTGCGATCGATCACCTGGCAGGCGTCCAGCACGGTAAAGCCCATGACTTAATCAAACACGGTACGAATACaaaacaga                                                                >  2:3658463/1‑145 (MQ=255)
                                                                                    cGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCa                                                               <  2:4742614/150‑1 (MQ=255)
                                                                                        atCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGa                                                           <  1:1320517/150‑1 (MQ=255)
                                                                                         tCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGaa                                                          >  2:2352064/1‑150 (MQ=255)
                                                                                          cAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAAt                                                         <  1:298079/150‑1 (MQ=255)
                                                                                          cAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAAt                                                         <  2:1311209/150‑1 (MQ=255)
                                                                                           aaaaGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATg                                                        <  2:3465652/150‑1 (MQ=255)
                                                                                                tCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCg                                                   <  1:1972840/150‑1 (MQ=255)
                                                                                                  aGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCggg                                                 >  2:4453077/1‑150 (MQ=255)
                                                                                                  aGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCggg                                                 <  2:1733038/150‑1 (MQ=255)
                                                                                                     tAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCg                                              >  2:1229702/1‑150 (MQ=255)
                                                                                                      aaaCCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGa                                             >  1:619659/1‑150 (MQ=255)
                                                                                                       aaCCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGAt                                            >  1:1149121/1‑150 (MQ=255)
                                                                                                       aaCCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGAt                                            >  2:1591298/1‑150 (MQ=255)
                                                                                                          cGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAg                                         >  1:952552/1‑150 (MQ=255)
                                                                                                           ggCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGa                                        >  1:2344165/1‑150 (MQ=255)
                                                                                                            gCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGAc                                       >  2:3907548/1‑150 (MQ=255)
                                                                                                             cATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGACTGATTCGGGGCGATCAGAcc                                      >  1:162875/1‑150 (MQ=255)
                                                                                                               tatCACCCACACTCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACctt                                    <  1:4281176/150‑1 (MQ=255)
                                                                                                                atCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACcttc                                   <  2:2637179/150‑1 (MQ=255)
                                                                                                                atCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACcttc                                   >  1:4614468/1‑150 (MQ=255)
                                                                                                                atCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACcttc                                   <  1:1240166/150‑1 (MQ=255)
                                                                                                                atCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACcttc                                   >  1:65956/1‑150 (MQ=255)
                                                                                                                atCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACcttc                                   >  2:731012/1‑150 (MQ=255)
                                                                                                                atCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACcttc                                   <  2:1884587/150‑1 (MQ=255)
                                                                                                                      cacaATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCgg                             <  2:48868/150‑1 (MQ=255)
                                                                                                                         aaTCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGgcgc                          <  2:4220139/150‑1 (MQ=255)
                                                                                                                              gTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGtat                     <  1:2832208/150‑1 (MQ=255)
                                                                                                                               tttCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGtata                    >  2:3079167/1‑150 (MQ=255)
                                                                                                                                  cACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGGTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGaa                 <  2:3291633/150‑1 (MQ=255)
                                                                                                                                     ttGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACAt              <  2:436606/150‑1 (MQ=255)
                                                                                                                                      tGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATa             <  2:1837782/150‑1 (MQ=255)
                                                                                                                                          cGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGc         >  2:4013887/1‑150 (MQ=255)
                                                                                                                                           gATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGcc        <  2:4615125/150‑1 (MQ=255)
                                                                                                                                            aTAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCg       >  1:371044/1‑150 (MQ=255)
                                                                                                                                              aaCGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTg     <  1:2442501/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GTACCAGTTTAATAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTC  >  CP009273/3076379‑3076673

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