Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,913,996 T→G 8.1% Y744D (TAC→GAC)  yebT → MCE domain protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,913,9960TG 8.1% 660.6 / 47.6 234Y744D (TAC→GAC) yebTMCE domain protein
Reads supporting (aligned to +/- strand):  ref base T (106/109);  new base G (7/12);  total (113/121)
Fisher's exact test for biased strand distribution p-value = 3.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

Only 100 of 238 total aligned reads displayed.

CGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTACCTGGATGGCTTAAGCATTATTGTTGAAGCGCCGGAAGCCGGTTCGTTAGGCATCGGTACGCCTGTGCTGTTCCGTGGTCTGGAAGTCGGTACGGTTACAGGAATGACGCTGGGGACATTGTCAGATCGCGTGATGATTGCGATGC  >  CP009273/1913847‑1914143
                                                                                                                                                     |                                                                                                                                                   
cgcgCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGtt                                                                                                                                                     <  1:776267/150‑1 (MQ=255)
 gcgcTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTa                                                                                                                                                    <  1:3434071/150‑1 (MQ=255)
 gcgcTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTa                                                                                                                                                    <  1:3433776/150‑1 (MQ=255)
   gcTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTAcc                                                                                                                                                  >  2:3357410/1‑150 (MQ=255)
    cTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTACCt                                                                                                                                                 >  2:3720240/1‑150 (MQ=255)
     ttCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAAGTACAAGAGGCCACCATTACTGATTCGCGTTACCTg                                                                                                                                                >  2:804232/1‑150 (MQ=255)
      tctcAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTACCTgg                                                                                                                                               <  2:3497217/150‑1 (MQ=255)
      tctcAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTACCTgg                                                                                                                                               <  2:1223620/150‑1 (MQ=255)
                 acacCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTACCTGGATGGCTTAAGc                                                                                                                                    <  2:2453585/150‑1 (MQ=255)
                    ccGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTACCTGGATGGCTTAAGCatt                                                                                                                                 <  1:178686/150‑1 (MQ=255)
                    ccGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTACCTGGATGGCTTAAGCatt                                                                                                                                 <  1:178873/150‑1 (MQ=255)
                       cAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTACCTGGATGGCTTAAGCattatt                                                                                                                              <  1:558434/150‑1 (MQ=255)
                       cAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTACCTGGATGGCTTAAGCattatt                                                                                                                              >  2:2034130/1‑150 (MQ=255)
                       cAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTACCTGGATGGCTTAAGCattatt                                                                                                                              >  2:2035670/1‑150 (MQ=255)
                       cAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTACCTGGATGGCTTAAGCattatt                                                                                                                              >  2:3772960/1‑150 (MQ=255)
                         aaTTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTACCTGGATGGCTTAAGCATTAttgt                                                                                                                            >  2:473083/1‑150 (MQ=255)
                             tCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTACCTGGATGGCTTAAGCATTATTGTTGaa                                                                                                                        <  1:488568/150‑1 (MQ=255)
                             tCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTACCTGGATGGCTTAAGCATTATTGTTGaa                                                                                                                        <  1:488528/150‑1 (MQ=255)
                             tCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTACCTGGATGGCTTAAGCATTATTGTTGaa                                                                                                                        <  1:488412/150‑1 (MQ=255)
                               ggCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTACCTGGATGGCTTAAGCATTATTGTTGAAgc                                                                                                                      <  1:3044976/150‑1 (MQ=255)
                               ggCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTGACCTGGATGGCTTAAGCATTATTGTTGAAgc                                                                                                                      <  2:1670258/150‑1 (MQ=255)
                                 cAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTACCTGGATGGCTTAAGCATTATTGTTGAAgcgc                                                                                                                    >  1:2487414/1‑150 (MQ=255)
                                 cAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTACCTGGATGGCTTAAGCATTATTGTTGAAgcgc                                                                                                                    >  1:188213/1‑150 (MQ=255)
                                   gCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTACCTGGATGGCTTAAGCATTATTGTTGAAGCGCCg                                                                                                                  <  2:396158/150‑1 (MQ=255)
                                     tGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTGACCTGGATGGCTTAAGCATTATTGTTGAAGCGCCGGa                                                                                                                <  1:2464966/150‑1 (MQ=255)
                                     tGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTGACCTGGATGGCTTAAGCATTATTGTTGAAGCGCCGGa                                                                                                                <  1:2464959/150‑1 (MQ=255)
                                         gTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTACCTGGATGGCTTAAGCATTATTGTTGAAGCGCCGGAAGcc                                                                                                            <  1:1948493/150‑1 (MQ=255)
                                          ttGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTACCTGGATGGCTTAAGCATTATTGTTGAAGCGCCGGAAGCCg                                                                                                           >  2:3756110/1‑150 (MQ=255)
                                            gAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTACCTGGATGGCTTAAGCATTATTGTTGAAGCGCCGGAAGCCGGt                                                                                                         <  2:1943763/150‑1 (MQ=255)
                                             aGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTACCTGGATGGCTTAAGCATTATTGTTGAAGCGCCGGAAGCCGGtt                                                                                                        <  1:2904085/150‑1 (MQ=255)
                                             aGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTACCTGGATGGCTTAAGCATTATTGTTGAAGCGCCGGAAGCCGGtt                                                                                                        <  2:1696716/150‑1 (MQ=255)
                                               cATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTACCTGGATGGCTTAAGCATTATTGTTGAAGCGCCGGAAGCCGGTTCg                                                                                                      <  2:2674455/150‑1 (MQ=255)
                                               cATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTACCTGGATGGCTTAAGCATTATTGTTGAAGCGCCGGAAGCCGGTTCg                                                                                                      <  2:10658/150‑1 (MQ=255)
                                               cATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTGACCTGGATGGCTTAAGCATTATTGTTGAAGCGCCGGAAGCCGGTTCg                                                                                                      <  1:3693869/150‑1 (MQ=255)
                                                aTCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTACCTGGATGGCTTAAGCATTATTGTTGAAGCGCCGGAAGCCGGTTCGt                                                                                                     >  2:3903644/1‑150 (MQ=255)
                                                 tCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTACCTGGATGGCTTAAGCATTATTGTTGAAGCGCCGGAAGCCGGTTCGtt                                                                                                    <  2:505488/150‑1 (MQ=255)
                                                 tCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTACCTGGATGGCTTAAGCATTATTGTTGAAGCGCCGGAAGCCGGTTCGtt                                                                                                    <  2:4299055/150‑1 (MQ=255)
                                                      aTACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTACCTGGATGGCTTAAGCATTATTGTTGAAGCGCCGGAAGCCGGTTCGTTAGGCa                                                                                               <  1:110006/150‑1 (MQ=255)
                                                      aTACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTACCTGGATGGCTTAAGCATTATTGTTGAAGCGCCGGAAGCCGGTTCGTTAGGCa                                                                                               <  2:3579217/150‑1 (MQ=255)
                                                      aTACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTACCTGGATGGCTTAAGCATTATTGTTGAAGCGCCGGAAGCCGGTTCGTTAGGCa                                                                                               <  2:2559468/150‑1 (MQ=255)
                                                      aTACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTACCTGGATGGCTTAAGCATTATTGTTGAAGCGCCGGAAGCCGGTTCGTTAGGCa                                                                                               <  1:2414393/150‑1 (MQ=255)
                                                       tACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTACCTGGATGGCTTAAGCATTATTGTTGAAGCGCCGGAAGCCGGTTCGTTAGGCAt                                                                                              <  2:2121885/150‑1 (MQ=255)
                                                        aCTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTACCTGGATGGCTTAAGCATTATTGTTGAAGCGCCGGAAGCCGGTTCGTTAGGCATc                                                                                             >  1:2611715/1‑150 (MQ=255)
                                                          tATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTACCTGGATGGCTTAAGCATTATTGTTGAAGCGCCGGAAGCCGGTTCGTTAGGCATCgg                                                                                           >  1:556888/1‑150 (MQ=255)
                                                          tATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTACCTGGATGGCTTAAGCATTATTGTTGAAGCGCCGGAAGCCGGTTCGTTAGGCATCgg                                                                                           >  1:2841157/1‑150 (MQ=255)
                                                            tcctccAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGAGTCGCGTTACCTGGATGGCTTAAGCATTATTGTTGAAGCGCCGGAAGCCGGTTCGTTAGGCATCGGTa                                                                                         <  2:3996698/150‑1 (MQ=255)
                                                             cctccAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTACCTGGATGGCTTAAGCATTATTGTTGAAGCGCCGGAAGCCGGTTCGTTAGGCATCGGTAc                                                                                        <  2:1459946/150‑1 (MQ=255)
                                                             cctccAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTACCTGGATGGCTTAAGCATTATTGTTGAAGCGCCGGAAGCCGGTTCGTTAGGCATCGGTAc                                                                                        >  1:2060139/1‑150 (MQ=255)
                                                               tccAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTACCTGGATGGCTTAAGCATTATTGTTGAAGCGCCGGAAGCCGGTTCGTTAGGCATCGGTACGc                                                                                      <  1:1346551/150‑1 (MQ=255)
                                                                  aGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTACCTGGATGGCTTAAGCATTATTGTTGAAGCGCCGGAAGCCGGTTCGTTAGGCATCGGTACGCCtg                                                                                   >  2:762003/1‑150 (MQ=255)
                                                                  aGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTACCTGGATGGCTTAAGCAATATTGTTGAAGCGCCGGAAGCCGGTTCGTTAGGCATCGGTACGCCtg                                                                                   <  1:3315173/150‑1 (MQ=255)
                                                                      gTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTACCTGGATGGCTTAAGCATTATTGTTGAAGCGCCGGAAGCCGGTTCGTTAGGCATCGGTACGCCTGTGCt                                                                               >  1:3818589/1‑150 (MQ=255)
                                                                      gTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTACCTGGATGGCTTAAGCATTATTGTTGAAGCGCCGGAAGCCGGTTCGTTAGGCATCGGTACGCCTGTGCt                                                                               <  1:3859168/150‑1 (MQ=255)
                                                                       tataTCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTACCTGGATGGCTTAAGCATTATTGTTGAAGCGCCGGAAGCCGGTTCGTTAGGCATCGGTACGCCTGTGCTg                                                                              <  1:1811090/150‑1 (MQ=255)
                                                                       tataTCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTACCTGGATGGCTTAAGCATTATTGTTGAAGCGCCGGAAGCCGGTTCGTTAGGCATCGGTACGCCTGTGCTg                                                                              >  1:48039/1‑150 (MQ=255)
                                                                       tataTCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTACCTGGATGGCTTAAGCATTATTGTTGAAGCGCCGGAAGCCGGTTCGTTAGGCATCGGTACGCCTGTGCTg                                                                              >  1:2861524/1‑150 (MQ=255)
                                                                           tCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTACCTGGATGGCTTAAGCATTATTGTTGAAGCGCCGGAAGCCGGTTCGTTAGGCATCGGTACGCCTGTGCTGTTcc                                                                          <  2:3715478/150‑1 (MQ=255)
                                                                            cAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTACCTGGATGGCTTAAGCATTATTGTTGAAGCGCCGGAAGCCGGTTCGTTAGGCATCGGTACGCCTGTGCTGTTCCg                                                                         >  1:1991618/1‑150 (MQ=255)
                                                                            cAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTACCTGGATGGCTTAAGCATTATTGTTGAAGCGCCGGAAGCCGGTTCGTTAGGCATCGGTACGCCTGTGCTGTTCCg                                                                         >  1:1991466/1‑150 (MQ=255)
                                                                              aCGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTACCTGGATGGCTTAAGCATTATTGTTGAAGCGCCGGAAGCCGGTTCGTTAGGCATCGGTACGCCTGTGCTGTTCCGTg                                                                       >  2:3784994/1‑150 (MQ=255)
                                                                               cGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTACCTGGATGGCTTAAGCATTATTGTTGAAGCGCCGGAAGCCGGTTCGTTAGGCATCGGTACGCCTGTGCTGTTCCGTgg                                                                      >  2:416935/1‑150 (MQ=255)
                                                                                    aaCCAGGCCGGGGCCATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTACCTGGATGGCTTAAGCATGATTGTTGAAGCGCCGGAAGCCGGTTCGTTAGGCATCGGTACGCCTGTGCTGTTCCGTGGTCTgg                                                                 >  2:1367083/1‑150 (MQ=255)
                                                                                    aaCCAGGCCGGGGCAATCCTCGCCGCGACTTTGTATTACAAGAGGCCACCATTACTGATTCGCGTTACCTGGATGGCTTAAGCATTATTGTTGAAGCGCCGGAAGCCGGTTCGTTAGGCATCGGTACGCCTGTGCTGTTCCGTGGTCTgg                                                                 <  1:3971758/150‑1 (MQ=255)
                                                                                    aaCCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTACCTGGATGGCTTAAGCATTATTGTTGAAGCGCCGGAAGCCGGTTCGTTAGGCATCGGTACGCCTGTGCTGTTCCGTGGTCTgg                                                                 <  1:3791100/150‑1 (MQ=255)
                                                                                       cAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTACCTGGATGGCTTAAGCATTATTGTTGAAGCGCCGGAAGCCGGTTCGTTAGGCATCGGTACGCCTGTGCTGTTCCGTGGTCTGGAAg                                                              >  1:1442565/1‑150 (MQ=255)
                                                                                        aGGCCGTGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTACCTGGATGGCTTAAGCATTATTGTTGAAGCGCCGGAAGCCGGTTCGTTAGGCATCGGTACGCCTGTGCTGTTCCGTGGTCTGGAAGt                                                             >  1:2094183/1‑150 (MQ=255)
                                                                                        aGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTACCTGGATGGCTTAAGCATTATTGTTGAAGCGCCGGAAGCCGGTTCGTTAGGCATCGGTACGCCTGTGCTGTTCCGTGGTCTGGAAGt                                                             >  2:3288977/1‑150 (MQ=255)
                                                                                            cGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTACCTGGATGGCTTAAGCATTATTGTTGAAGCGCCGGAAGCCGGTTCGTTAGGCATCGGTACGCCTGTGCTGTTCCGTGGTCTGGAAGTCGGt                                                         <  2:3248789/150‑1 (MQ=255)
                                                                                            cGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTGACCTGGATGGCTTAAGCATTATTGTTGAAGCGCCGGAAGCCGGTTCGTTAGGCATCGGTACGCCTGTGCTGTTCCGTGGTCTGGAAGTCGGt                                                         <  1:391807/150‑1 (MQ=255)
                                                                                                gCAATCCTCGCCGCGACTTTGAATTACAATAGGCCACCATTACTGATTCGCGTTACCTGGATGGCTTAAGCATTATTGTTGAAGCGCCGGAAGCCGGTTCGTTAGGCATCGGTACGCCTGTGCTGTTCCGTGGTCTGGAAGTCGGTACgg                                                     >  1:3354332/1‑150 (MQ=255)
                                                                                                gCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTACCTGGATGGCTTAAGCATTATTGTTGAAGCGCCGGAAGCCGGTTCGTTAGGCATCGGTACGCCTGTGCTGTTCCGTGGTCTGGAAGTCGGTACgg                                                     >  1:1073226/1‑150 (MQ=255)
                                                                                                   aTCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTACCTGGATGGCTTAAGCATTATTGTTGAAGCGCCGGAAGCCGGTTCGTTAGGCATCGGTACGCCTGTGCTGTTCCGTGGTCTGGAAGTCGGTACGGTTa                                                  >  2:160907/1‑150 (MQ=255)
                                                                                                        cgccgcGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTACCTGGATGGCTTAAGCATTATTGTTGAAGCGCCGGAAGCCGGTTCGTTAGGCATCGGTACGCCTGTGCTGTTCCGTGGTCTGGAAGTCGGTACGGTTACAGGa                                             <  2:101112/150‑1 (MQ=255)
                                                                                                         gccgcGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTACCTGGATGGCTTAAGCATTATTGTTGAAGCGCCGGAAGCCGGTTCGTTAGGCATCGGTACGCCTGTGCTGTTCCGTGGTCTGGAAGTCGGTACGGTTACAGGaa                                            >  2:2092671/1‑150 (MQ=255)
                                                                                                         gccgcGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTACCTGGATGGCTTAAGCATTATTGTTGAAGCGCCGGAAGCCGGTTCGTTAGGCATCGGTACGCCTGTGCTGTTCCGTGGTCTGGAAGTCGGTACGGTTACAGGaa                                            >  2:4004654/1‑150 (MQ=255)
                                                                                                         gccgcGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTACCTGGATGGCTTAAGCATTATTGTTGAAGCGCCGGAAGCCGGTTCGTTAGGCATCGGTACGCCTGTGCTGTTCCGTGGTCTGGAAGTCGGTACGGTTACAGGaa                                            >  2:1750813/1‑150 (MQ=255)
                                                                                                         gccgcGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTACCTGGATGGCTTAAGCATTATTGTTGAAGCGCCGGAAGCCGGTTCGTTAGGCATCGGTACGCCTGTGCTGTTCCGGGGTCTGGAAGTCGGTACGGTTACAGGaa                                            >  1:3703969/1‑150 (MQ=255)
                                                                                                             cgACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTGACCTGGATGGCTTAAGCATTATTGTTGAAGCGCCGGAAGCCGGTTCGTTAGGCATCGGTACGCCTGTGCTGTTCCGTGGTCTGGAAGTCGGTACGGTTACAGGAATGAc                                        >  1:1163194/1‑150 (MQ=255)
                                                                                                                 tttGAATTACAAGAGGCCACCATTACTGATTCGCGTTACCTGGATGGCTTAAGCATTATTGTTGAAGCGCCGGAAGCCGGTTCGTTAGGCCTCGGTACGCCTGTGCTGTTCCGTGGTCTGGAAGTCGGTACGGTTACAGGAATGACGCTg                                    >  2:59493/1‑150 (MQ=255)
                                                                                                                   tGAATTACAAGAGGCCACCATTACTGATTCGCGTTACCTGGATGGCTTAAGCATTATTGTTGAAGCGCCGGAAGCCGGTTCGTTAGGCATCGGTACGCCTGTGCTGTTCCGTGGTCTGGAAGTCGGTACGGTTACAGGAATGACGCTggg                                  <  1:4402868/150‑1 (MQ=255)
                                                                                                                      aTTACAAGAGGCCACCATTACTGATTCGCGTTACCTGGATGGCTTAAGCATTATTGTTGAAGCGCCGGAAGCCGGTTCGTTAGGCATCGGTACGCCTGTGCTGTTCCGTGGTCTGGAAGTCGGTACGGTTACAGGAATGACGCTGGGGAc                               >  2:3262960/1‑150 (MQ=255)
                                                                                                                      aTTACAAGAGGCCACCATTACTGATTCGCGTTACCTGGATGGCTTAAGCATTATTGTTGAAGCGCCGGAAGCCGGTTCGTTAGGCATCGGTACGCCTGTGCTGTTCCGTGGTCTGGAAGTCGGTACGGTTACAGGAATGACGCTGGGGAc                               >  2:3262478/1‑150 (MQ=255)
                                                                                                                         aCAAGAGGCCACCATTACTGATTCGCGTTACCTGGATGGCTTAAGCATTATTGTTGAAGCGCCGGAAGCCGGTTCGTTAGGCATCGGTACGCCTGTGCTGTTCCGTGGTCTGGAAGTCGGTACGGTTACAGGAATGACGCTGGGGACAtt                            >  1:3898439/1‑150 (MQ=255)
                                                                                                                          cAAGAGGCCACCATTACTGATTCGCGTTACCTGGATGGCTTAAGCATTATTGTTGAAGCGCCGGAAGCCGGTTCGTTAGGCATCGGTACGCCTGTGCTGTTCCGTGGTCTGGAAGTCGGTACGGTTACAGGAATGACGCTGGGGACATTg                           >  2:52641/1‑150 (MQ=255)
                                                                                                                            agagGCCACCATTACTGATTCGCGTTACCTGGATGGCTTAAGCATTATTGTTGAAGCGCCGGAAGCCCGTTCGTTAGGCATCGGTACGCCTGTGCTGTTCCGTGGTCTGGAAGTCGGTACGGTTACAGGAATGACGCTGGGGACATTGTc                         >  2:4371846/1‑150 (MQ=255)
                                                                                                                                 ccaccaTTACTGATTCGCGTTACCTGGATGGCTTAAGCATTATTGTTGAAGCGCCGGAAGCCGGTTCGTTAGGCATCGGTACGCCTGTGCTGTTCCGTGGTCTGGAAGTCGGTACGGTTACAGGAATGACGCTGGGGACATTGTCAGATc                    >  1:2656286/1‑150 (MQ=255)
                                                                                                                                  caccaTTACTGATTCGCGTTACCTGGATGGCTTAAGCATTATTGTTGAAGCGCCGGAAGCCGGTTCGTTAGGCATCGGTACGCCTGTGCTGTTCCGTGGTCTGGAAGTCGGTACGGTTACAGGAATGACGCTGGGGACATTGTCAGATcg                   <  1:1743917/150‑1 (MQ=255)
                                                                                                                                  caccaTTACTGATTCGCGTTACCTGGATGGCTTAAGCATTATTGTTGAAGCGCCGGAAGCCGGTTCGTTAGGCATCGGTACGCCTGTGCTGTTCCGTGGTCTGGAAGTCGGTACGGTTACAGGAATGACGCTGGGGACATTGTCAGATcg                   >  1:2248099/1‑150 (MQ=255)
                                                                                                                                     caTTACTGAGTCGCGTTACCTGGATGGCTTAAGCATTATTGTTGAAGCGCCGGAAGCCGGTTCGTTAGGCATCGGTACGCCTGTGCTGTTCCGTGGTCTGGAAGTCGGTACGGTTACAGGAATGACGCTGGGGACATTGTCAGATCGCGt                <  1:2982631/150‑1 (MQ=255)
                                                                                                                                      aTTACTGATTCGCGTTACCTGGATGGCTTAAGCATTATTGTTGAAGCGCCGGAAGCCGGTTCGTTAGGCATCGGTACGCCTGTGCTGTTCCGTGGTCTGGAAGTCGGTACGGTTACAGGAATGACGCTGGGGACATTGTCAGATCGCGtg               >  1:1536176/1‑150 (MQ=255)
                                                                                                                                          cTGATTCGCGTGACCTGGATGGCTTAAGCATTATTGTTGAAGCGCCGGAAGCCGGTTCGTTAGGCATCGGTACGCCTGTGCTGTTCCGTGGTCTGGAAGTCGGTACGGTTACAGGAATGACGCTGGGGACATTGTCAGATCGCGtgatga           >  2:4479273/1‑150 (MQ=255)
                                                                                                                                           tGATTCGCGTTACCTGGATGGCTTAAGCATTATTGTTGAAGCGCCGGAAGCCGGTTCGTTAGGCATCGGTACGCCTGTGCTGTTCCGTGGTCTGGAAGTCGGTACGGTTACAGGAATGACGCTGGGGACATTGTCAGATCGCGTgatgat          <  1:674710/150‑1 (MQ=255)
                                                                                                                                           tGATTCGCGTTACCTGGATGGCTTAAGCATTATTGTTGAAGCGCCGGAAGCCGGTTCGTTAGGCATCGGTACGCCTGTGCTGTTCCGTGGTCTGGAAGTCGGTACGGTTACAGGAATGACGCTGGGGACATTGTCAGATCGCGTgatgat          <  2:3845182/150‑1 (MQ=255)
                                                                                                                                            gATTCGCGTTACCTGGATGGCTTAAGCATTATTGTTGAAGCGCCGGAAGCCGGTTCGTTAGGCATCGGTACGCCTGTGCTGTTCCGTGGTCTGGAAGTCGGTACGGTTACAGGAATGACGCTGGGGACATTGTCAGATCGCGTGATGAtt         <  2:4439103/150‑1 (MQ=255)
                                                                                                                                             aTTCGCGTGACCTGGATGGCTTAAGCATTATTGTTGAAGCGCCGGAAGCCGGTTCGTTAGGCATCGGTACGCCTGTGCTGTTCCGTGGTCTGGAAGTCGGTACGGTTACAGGAATGACGCTGGGGACATTGTCAGATCGCGTGATGATTg        <  2:4134649/150‑1 (MQ=255)
                                                                                                                                             aTTCGCGTGACCTGGATGGCTTAAGCATTATTGTTGAAGCGCCGGAAGCCGGGTCGTTAGGCATCGGTACGCCTGTGCTGTTCCGTGGTCTGGAAGTCGGTACGGTTACAGGAATGACGCTGGGGACATTGTCAGATCGCGTGATGATTg        <  2:2066638/150‑1 (MQ=255)
                                                                                                                                             aTTCGAGTGACCTGGATGGCTTAAGCATTATTGTTGTAGCGCCGGAAGCCGGTTCGTTAGGCATCGGTACGCCTGTGCTGTTCCGTGGTCTGGAAGTCGGTACGGTTACAGGAATGACGCTGGGGACATTGTCAGATCGCGTGATGATTg        <  2:4134621/150‑1 (MQ=255)
                                                                                                                                                cgcgTTACCTGGATGGCTTAAGCATTATTGTTGAAGCGCCGGAAGCCGGTTCGTTAGGCATCGGTACGCCTGTGCTGTTCCGTGGTCTGGAAGTCGGTACGGTTACAGGAATGACGCTGGGGACATTGTCAGATCGCGTGATGATTGCGa     <  2:3578400/150‑1 (MQ=255)
                                                                                                                                                 gcgTTACCTGGATGGCTTAAGCATTATTGTTGAAGCGCCGGAAGCCGGTTCGTTAGGCATCGGTACGCCTGTGCTGTTCCGTGGTCTGGAAGTCGGTACGGTTACAGGAATGACGCTGGGGACATTGTCAGATCGCGTGATGATTGCGAt    <  1:473041/150‑1 (MQ=255)
                                                                                                                                                   gTTACCTGGATGGCTTAAGCATTATTGTTGAAGCGCCGGAAGCCGGTTCGTTAGGCATCGGTACGCCTGTGCTGTTCCGTGGTCTGGAAGTCGGTACGGTTACAGGAATGACGCTGGGGACATTGTCAGATCGCGTGATGATTGCGATgc  >  1:2218836/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTACCTGGATGGCTTAAGCATTATTGTTGAAGCGCCGGAAGCCGGTTCGTTAGGCATCGGTACGCCTGTGCTGTTCCGTGGTCTGGAAGTCGGTACGGTTACAGGAATGACGCTGGGGACATTGTCAGATCGCGTGATGATTGCGATGC  >  CP009273/1913847‑1914143

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