Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 8,600 T→C 66.8% D121D (GAT→GAC)  talB → transaldolase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092738,6000TC66.8% 207.1 / 184.9 184D121D (GAT→GAC) talBtransaldolase B
Reads supporting (aligned to +/- strand):  ref base T (29/32);  new base C (56/67);  total (85/99)
Fisher's exact test for biased strand distribution p-value = 8.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.60e-24

Only 100 of 184 total aligned reads displayed.

CGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCG  >  CP009273/8453‑8747
                                                                                                                                                   |                                                                                                                                                   
   gaccgacAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTgg                                                                                                                                                >  1:322707/1‑150 (MQ=255)
        acAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCAAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTAt                                                                                                                                             >  2:382544/1‑148 (MQ=255)
             cTGGCGGTGAACATCGGTCTCGAAATTCTGAAACTGATCCCGGGCCGTATCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAAc                                                                                                                                      >  1:15351/1‑150 (MQ=255)
                  aGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCg                                                                                                                                 >  2:481254/1‑150 (MQ=255)
                    tAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTa                                                                                                                               >  1:819232/1‑150 (MQ=255)
                      aacatcggTCTCGAAATTCTGAAACTGATCCCGGGCCGTATCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCAtt                                                                                                                             <  2:925892/144‑1 (MQ=255)
                      aacatcggTCTCGAAATTCTGAAACTGATCCCGGGCCGTATCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCAt                                                                                                                              <  2:319030/143‑1 (MQ=255)
                       atatTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTc                                                                                                                            <  1:519728/150‑1 (MQ=255)
                        catcggTCTCGAAATTCTGAAACTGATCCCGGGCCGTATCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCt                                                                                                                           >  1:101974/5‑150 (MQ=255)
                         atcggTCTCGAAATTCTGAAACTGATCCCGGGCCGTATCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCt                                                                                                                           >  2:572725/4‑149 (MQ=255)
                         atTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCt                                                                                                                           >  1:189461/1‑149 (MQ=255)
                          tcggTCTCGAAATTCTGAAACTGATCCCGGGCCGTATCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGa                                                                                                                         <  2:709733/148‑1 (MQ=255)
                             gTCTCGAAATTCTGAAACTGATCCCGGGCCGTATCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCa                                                                                                                      <  2:154375/150‑1 (MQ=255)
                                 cgAAATTCTGAAACTGATCCCGGGCCGTATCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACt                                                                                                                  >  1:15218/2‑150 (MQ=255)
                                  gAAATTCTGAAACTGATCCCGGGCCGTATCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACt                                                                                                                  <  1:895760/149‑1 (MQ=255)
                                  gAAATTCTGAAACTGATCCCGGGCCATATCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACt                                                                                                                  <  2:241912/149‑1 (MQ=255)
                                  gAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACt                                                                                                                  <  2:870582/149‑1 (MQ=255)
                                   aaaTCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTgg                                                                                                                <  1:427874/150‑1 (MQ=255)
                                    aaTTCTGAAACTGATCCCGGGCCGTATCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGc                                                                                                               >  1:395826/1‑150 (MQ=255)
                                     attcTGAAACTGATCCCGGGCCGTATCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGcc                                                                                                              >  1:282456/4‑149 (MQ=255)
                                     attcTGAAACTGATCCCGGGCCGTATCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGcc                                                                                                              >  1:549645/4‑149 (MQ=255)
                                      ttcTGAAACTGATCCCGGGCCGTATCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGcct                                                                                                             <  1:938442/148‑3 (MQ=255)
                                         tGAAACTGATCCCGGGCCGTATCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTcca                                                                                                          >  1:945407/1‑148 (MQ=255)
                                             aCTGATCCCGGGCCGTATCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTg                                                                                                      <  2:714272/150‑1 (MQ=255)
                                                 atcccGGGCCGTATCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAg                                                                                                  <  1:141450/147‑1 (MQ=255)
                                                       ggCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATc                                                                                            >  1:599510/1‑150 (MQ=255)
                                                           gTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGtg                                                                                        >  1:12935/1‑150 (MQ=255)
                                                            tATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGtgc                                                                                       >  1:553146/1‑150 (MQ=255)
                                                            tATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGtgc                                                                                       <  1:796210/150‑1 (MQ=255)
                                                             aTCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCgcgcc                                                                                      >  2:353990/1‑149 (MQ=255)
                                                              tctcCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCgcgcc                                                                                      <  2:986450/149‑2 (MQ=255)
                                                                tcAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGtgctg                                                                                     >  1:560642/1‑148 (MQ=255)
                                                                 ccaccgAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGcgccgct                                                                                   >  2:274696/6‑148 (MQ=255)
                                                                  caccgAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGcgccgctg                                                                                  <  1:209908/145‑3 (MQ=255)
                                                                  caccgAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGcgccgctg                                                                                  <  1:682371/145‑3 (MQ=255)
                                                                    ccgAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGaa                                                                                >  1:170670/3‑149 (MQ=255)
                                                                    cTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAAc                                                                               <  1:978900/150‑1 (MQ=255)
                                                                    cTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAAc                                                                               <  2:424896/150‑1 (MQ=255)
                                                                         gTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTg                                                                          <  1:715074/150‑1 (MQ=255)
                                                                             acgcacgTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGaa                                                                       >  2:802604/6‑149 (MQ=255)
                                                                                cgcgTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGaaa                                                                      >  1:663844/1‑147 (MQ=255)
                                                                                 gcgTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAaaga                                                                  <  2:579709/150‑1 (MQ=255)
                                                                                 gcgTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAaaga                                                                  <  1:707099/150‑1 (MQ=255)
                                                                                 acgTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAaag                                                                   >  2:441445/2‑149 (MQ=255)
                                                                                 acgTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAaag                                                                   >  1:93407/2‑149 (MQ=255)
                                                                                       ctCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCAt                                                            >  2:93557/2‑150 (MQ=255)
                                                                                        tCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCa                                                             <  1:639691/148‑1 (MQ=255)
                                                                                        tCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCa                                                             <  1:831789/148‑1 (MQ=255)
                                                                                        tCCTATGACACCGAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCa                                                             <  1:414239/148‑1 (MQ=255)
                                                                                          cTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCAtt                                                           <  2:642492/148‑2 (MQ=255)
                                                                                          cTATGACACCGAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCaa                                                         >  1:401619/1‑150 (MQ=255)
                                                                                          cTAAGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCAtta                                                          <  2:395826/149‑3 (MQ=255)
                                                                                            aTGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACt                                                       >  2:109455/1‑150 (MQ=255)
                                                                                            aTGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACt                                                       >  1:692024/1‑150 (MQ=255)
                                                                                             tGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACt                                                       <  1:904880/149‑1 (MQ=255)
                                                                                             tGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTg                                                      <  2:294147/150‑1 (MQ=255)
                                                                                              gACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTg                                                      >  2:726019/1‑149 (MQ=255)
                                                                                              gACACCGAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTg                                                      >  1:756412/1‑149 (MQ=255)
                                                                                               acacCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTa                                                    <  1:460794/150‑1 (MQ=255)
                                                                                               acacCGAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTa                                                    <  2:553146/150‑1 (MQ=255)
                                                                                                cacCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTaa                                                   >  2:141530/1‑150 (MQ=255)
                                                                                                  ccGAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAAcc                                                 <  2:647144/150‑1 (MQ=255)
                                                                                                  ccGAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAAcc                                                 <  2:12935/150‑1 (MQ=255)
                                                                                                   cGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGCACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCt                                                >  1:789607/1‑150 (MQ=255)
                                                                                                    gAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTg                                               <  2:46284/150‑1 (MQ=255)
                                                                                                      aGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGAc                                             <  2:223235/150‑1 (MQ=255)
                                                                                                        cGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGAcgc                                           >  1:300901/1‑148 (MQ=255)
                                                                                                        cGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGAcgc                                           <  1:101803/150‑3 (MQ=255)
                                                                                                         gTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGct                                          <  1:349242/150‑1 (MQ=255)
                                                                                                         gTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGct                                          <  1:112324/150‑1 (MQ=255)
                                                                                                         gTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGct                                          >  1:405/1‑150 (MQ=255)
                                                                                                          tctatcgCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGctg                                         <  2:821156/144‑1 (MQ=255)
                                                                                                              tcgCGAAAGCCAAGCGCATCATCAAGCTGTACAACGCCGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGctgctg                                      <  2:15218/147‑1 (MQ=255)
                                                                                                              tcgCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGctgctg                                      <  2:521670/147‑1 (MQ=255)
                                                                                                              tcgCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGctgctg                                      <  2:666305/147‑1 (MQ=255)
                                                                                                              tcgCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGt                                     <  1:380878/148‑1 (MQ=255)
                                                                                                               cgCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGtt                                    >  2:231087/2‑150 (MQ=255)
                                                                                                                gCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGtt                                    <  1:565898/149‑1 (MQ=255)
                                                                                                                gCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTtc                                   <  1:572725/150‑1 (MQ=255)
                                                                                                                  gAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTtctc                                 <  2:945407/150‑1 (MQ=255)
                                                                                                                  gAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTtctc                                 <  1:984350/150‑1 (MQ=255)
                                                                                                                  gAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTtctc                                 >  2:549671/1‑150 (MQ=255)
                                                                                                                    aaGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCt                               >  2:556175/1‑150 (MQ=255)
                                                                                                                     aGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCtt                              <  2:415496/150‑1 (MQ=255)
                                                                                                                        aaaaCGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCCCCTTCGc                           >  2:167985/1‑150 (MQ=255)
                                                                                                                           gcgcatcaTCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCt                          >  2:674542/8‑148 (MQ=255)
                                                                                                                            cgcatcaTCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAg                       >  1:13626/7‑150 (MQ=255)
                                                                                                                            cgcatcaTCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAg                       <  2:884985/144‑1 (MQ=255)
                                                                                                                            cGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAg                       >  2:758561/1‑150 (MQ=255)
                                                                                                                             gcatcaTCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAgg                      <  1:536294/145‑1 (MQ=255)
                                                                                                                             gcatcaTCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAgg                      <  1:763688/145‑1 (MQ=255)
                                                                                                                              catcaTCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGc                     >  2:193805/5‑150 (MQ=255)
                                                                                                                               atcaTCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCt                    >  1:218176/4‑150 (MQ=255)
                                                                                                                                      aaCTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCtt             >  1:437692/1‑150 (MQ=255)
                                                                                                                                       gctgtACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTtg            >  2:548585/5‑150 (MQ=255)
                                                                                                                                        ctgtACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTtgc           >  2:683559/4‑149 (MQ=255)
                                                                                                                                           tACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTgcgcc        <  2:746036/150‑2 (MQ=255)
                                                                                                                                           tACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTgcgcc        >  2:964113/1‑149 (MQ=255)
                                                                                                                                             cAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTgcgccg       >  2:774050/1‑147 (MQ=255)
                                                                                                                                                cGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGc   >  1:26426/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCG  >  CP009273/8453‑8747

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