Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 3,722,437 G→A 6.8% D150D (GAC→GAT U712_19190 ← PTS system sucrose‑specific EIIBC component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068813,722,4370GA 6.8% 229.8 / 4.5 73D150D (GAC→GATU712_19190PTS system sucrose‑specific EIIBC component
Reads supporting (aligned to +/- strand):  ref base G (28/40);  new base A (2/3);  total (30/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.63e-01

CATGGATTCAAAAGCGTCGGATGAATCATGATACCTCCGATGACCGCTCCCAAGTAAGGGTTGCTTCCAAACTCTTTTGAAGCACTGACACCGATTAAAATCGGCAAGAAAATAAATGCCGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCATTGATCATACCCAATAATCCCATTAATAGGCCGCTGGCCACAATAGCCGGGATGATTGGAACAAAAATATTAGAAAGCGTTTTCGCAAAT  >  CP006881/3722303‑3722575
                                                                                                                                      |                                                                                                                                          
cATTGATTCAAAAGCGTCGGATGAATCATGATACCTCCGATGACCGCTCCCAAGTAAGGGTTGCTTCCAAACTCTTTTGAAGCACTGACACCGATTAAAATCGGCAAGAAAATAAATGCCGCACTTGAGAACATGTCGa                                                                                                                                        <  2:230152/139‑1 (MQ=255)
cATGGATTCAAAAGCGTCGGATGAATCATGATACCTCCGATGACCGCTCCCAAGTAAGGGTTGCTTCCAAACTCTTTTGAAGCACTGACACCGATTAAAATCGGCAAGAAAATAAATGCCGCACTTGAGAACATGTCGa                                                                                                                                        <  6:302202/139‑1 (MQ=255)
cATGGATTCAAAAGCGTCGGATGAATCATGATACCTCCGATGACCGCTCCCAAGTAAGGGTTGCTTACAAACTCTTTTGAAGCACTGACACCGATTAAAATCGGCAAGAAAATAAATGCCGCACTTGAGAACATGTCGa                                                                                                                                        <  2:76207/139‑1 (MQ=255)
   ggATTCAAAAGCGTCGGATGAATCATGATACCTCCGATGACCGCTCCCAAGTAAGGGTTGCTTCCAAACTCTTTTGAAGCACTGACACCGATTAAAATCGGCAAGAAAATAAATGCCGCACTTGAGAACATGTCGagc                                                                                                                                      <  2:339037/138‑1 (MQ=255)
        cAAAAGCGTCGGATGAATCATGATACCTCCGATGACCGCTCCCAAGTAAGGGTTGCTTCCAAACTCTTTTGAAGCCCTGACACCGATTAAAATCGGCAAGAAAATAAATGCCGCACTTGAGAACATGTCGAGCAGCTgc                                                                                                                                >  8:80499/1‑139 (MQ=255)
            aGCGGCGGATGAATCATGATACCTCCGATGACCGCTCCCAAGTAAGGGTTGCTTCCAAACTCTTTTGAAGCACTGACACCGATTAAAATCGGCAAGAAAATAAATGCCGCACTTGAGAATATGTCGAGCAGCTgcagca                                                                                                                            <  8:346902/139‑1 (MQ=255)
              cGTCGGATGAATCATGATACCTCCGATGACCGCTCCCAAGTAAGGGTTGCTTCCAAACTCTTTTGAAGCACTGACACCGATTAAAATCGGCAAGAAAATAAATGCCGCACTTGAGAACATGTCGAGCAGCTGCagcagc                                                                                                                          >  1:21777/1‑139 (MQ=255)
                  ggATGAATCATGATACCTCCGATGACCGCTCCCAAGTAAGGGTTGCTTCCAAACTCTTTTGAAGCACTGACACCGATTAAAATCGGCAAGAAAATAAATGCCGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCa                                                                                                                       <  6:148927/138‑1 (MQ=255)
                  ggAAGAATCATGATTCCCCCGATGACCGCTCCCAAGTATGGATTGCTGCCAAACTCTTTTGAAGCGCTGACCCCGATTAGAATCGGCAAGAAAATGAATGCTGCACTTGAAAACATATCAAGCAGCTGCAAAAGCGCAg                                                                                                                      <  4:228730/139‑1 (MQ=255)
                     tGAATCATGATACCTCCGATGACCGCTCCCAAGTAAGGGTTGCTTCCAAACTCTTTTGAAGCACTGACACCGATTAAAATCGGCAAGAAAATAAATGCCGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAAt                                                                                                                   <  2:120336/139‑1 (MQ=255)
                           aTGATACCTCCGATGACCGCTCCCAAGTAAGGGTTGCTTCCAAACTCTTTTGAAGCACTGACACCGATTAAAATCGGCAAGAAAATAAATGCCGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGc                                                                                                             <  8:38257/139‑1 (MQ=255)
                                 ccTCCGATGACCGCTCCCAAGTAAGGGTTGCTTCCAAACTCTTTTGAAGCACTGACACCGATTAAAATCGGCAAGAAAATAAATGCCGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATc                                                                                                        <  3:154343/138‑1 (MQ=255)
                                  cTCCGCTGACCGCTCCCAAGTAAGGGTTGCTTCCAAACTCTGTTGACGCACTGACACCGATTAAAATCGGCAAGAAAATAAATGACGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATcc                                                                                                       >  2:258283/1‑138 (MQ=255)
                                    ccGATGACCGCTCCCAAGTAAGGGTTGCTTCCAAACTCTTTTGAAGCACTGACACCGATTAAAATCGGCAAGAAAATAAATGCCGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAAt                                                                                                    >  7:4465/1‑139 (MQ=255)
                                       aTGACCGCTCCCAAGTATGGATTGCTGCCAAACTCTTTTGAAGCGCTGACCCCGATTAGAATCGGCAAGAAAATGAATGCTGCACTTGAAAACATATCAAGCAGCTGCAAAAGCGCAGAATCCTTGCTCATCCAATGaa                                                                                                 <  4:140674/139‑1 (MQ=255)
                                       aTGACCGCTCCCAAGTAAGGGTTGCTTCCAAACTCTTTTGAAGCACTGACACCGATTAAAATCGGCAAGAAAATAAATGCCGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGa                                                                                                  >  2:242875/1‑138 (MQ=255)
                                       aTGACAGCGCCCCAGTAAGGGTTGCCTCCAAACTCTTTTGACGCACGGACACCGATTAACATCGGCAAGAAAATAAAAGCCGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTAGCTCATCCAATGaa                                                                                                 <  6:23463/139‑1 (MQ=255)
                                            cGCTCCCAAGTAAGGGTTGCTTCCAAACTCTTTTGAAGCACTGACACCGATTAAAATCGGCAAGAAAATAAATGCCGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCa                                                                                            >  7:72115/1‑139 (MQ=255)
                                            cGCTCCCAAGTAAGGGTTGCTTCCAAACTCTTTTGAAGCACTGACACCGATTAAAATCGGCAAGAAAATAAATGCCGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCa                                                                                            <  2:294128/139‑1 (MQ=255)
                                            cGCTCCCAAGTAAGGGTTGCTTCCAAACTCTTTTGAAGCACTGACACCGATTAAAATCGGCAAGAAAATAAATGCCGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCa                                                                                            <  3:294479/139‑1 (MQ=255)
                                              cTCCCAAGTATGGATTGCTGCCAAACTCTTTTGAAGCGCTGACCCCGATTAGAATCGGCAAGAAAATGAATGCTGCACTTGAAAACATATCAAGCAGCTGCAAAAGCGCAGAATCCTTGCTCATCCAATGAAACGCAtt                                                                                          >  3:150166/1‑139 (MQ=255)
                                              cTCCCAAGTAAGGGTTGCTTCCAAACTCTTTTGAAGCACTGACACCGATTAAAATCGGCAAGAAAATAACTGCCGCACTTGAGCACATGTCGAGCAGCTGCAGCAGCGCAGCATCCTTGCTCATCCAATGAAACGCAtt                                                                                          >  5:71663/1‑139 (MQ=255)
                                              cTCCCAAGTAAGGGTTGCTTCCAAACTCTTTTGAAGCACTGACACCGATTAAAATCGGCAAGAAAATAAATGCCGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCa                                                                                            <  5:293440/137‑1 (MQ=255)
                                              cTCCCAAGTAAGGGTTGCTTCCAAACTCTTTTGAAGCACTGACACCGATTAAAATCGGCAAGAAAATAAATGCCGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCa                                                                                            >  6:293440/1‑137 (MQ=255)
                                              cTCCCAAGTAAGGGTTGCTTCCAAACTCTTTTGAAGCACTGACACCGATTAAAATCGGCAAGAAAATAAATGCCGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCAtt                                                                                          >  1:242854/1‑139 (MQ=255)
                                              cTCCCAAGTAAGGGTTGCTTCCAAACTCTTTTGAAGCACTGACACCGATTAAAATCGGCAAGAAAATAAATGCCGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCAAATGAAACGCAtt                                                                                          >  8:314891/1‑139 (MQ=255)
                                               tCCCAAGTAAGGGTTGCTTCCAAACTCTTTTGAAGCACTGACACCGATTAAAATCGGCAAGAAAATAAATGCCGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCAtt                                                                                          >  4:247917/1‑138 (MQ=255)
                                                cccAAGTAAGGGGTGCTTCCTAACTCTTTTGAAGCACTGACACCAATTAAAATAGGCAAGAAAATAAATGCCGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCATTGa                                                                                        <  8:295219/139‑1 (MQ=255)
                                                 ccAAGTAAGGGTTGCTTCCAAACTCTTTTGAAGCACTGACACCGATTAAAATCGGCAAGAAAATAAATGCCGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCATTGAt                                                                                       <  2:21777/139‑1 (MQ=255)
                                                     gTAAGGGTTGCTTCCAAACTCTTTTGAAGCACTGACACCTATTAAAATCGGCAAGAAAATAAAAGCCGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCATTGATCATa                                                                                   <  1:53714/139‑1 (MQ=255)
                                                     gTAAGGGTTGCTTCCAAACTCTTTTGAAGCACTGACACCGTTTAAAATCGGCAAGAAAATAAATGCCGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCATTGATCATa                                                                                   <  3:6457/139‑1 (MQ=255)
                                                     gTAAGGGTTGCTTCCAAACTCTTTTGAAGCACTGACACCGATTAAAATCGGCAAGAAAATAAATGCCGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCATTGATCATa                                                                                   <  2:228706/139‑1 (MQ=255)
                                                     gTAAGGGTTGCTTCCAAACTCTTTTGAAGCACTGACACCGATTAAAATCGGCAAGAAAATAAATGCCGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCATTGATCATa                                                                                   <  8:4465/139‑1 (MQ=255)
                                                     gTAAGGGTTGCTTCCAAACTCTTTTGAAGCACTGACACCGATTAAAATCGGCAAGAAAATAAATGCCGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCATTGATCATa                                                                                   <  5:82105/139‑1 (MQ=255)
                                                         gggTTGCTTCCAAACTCTTTTGAAGCACTGACACCGATTAAAATCGGCAAGAAAATAAATGCCGCACTTGAGAACATGTCGagc                                                                                                                                      <  2:40020/84‑1 (MQ=255)
                                                         gggTTGCTTCCAAACTCTTTTGAAGCACTGACACCGATTAAAATCGGCAAGAAAATAAATGCCGCACTTGAGAACATGTCGagc                                                                                                                                      >  1:40020/1‑84 (MQ=255)
                                                         gggTTGCTTCCAAACTCTTTTGAAGCACTGACACCGATTAAAATCGGCAAGAAAATAAATGCCGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCa                                                                                                                       >  2:348344/1‑99 (MQ=255)
                                                         gggTTGCTTCCAAACTCTTTTGAAGCACTGACACCGATTAAAATCGGCAAGAAAATAAATGCCGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCa                                                                                                                       <  1:348344/99‑1 (MQ=255)
                                                         gggTTGCTTCCAAACTCTTTTGAAGCACTGACACCGATTAAAATCGGCAAGAAAATAAATGCCGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCATTGATCATAccc                                                                                >  2:285716/1‑138 (MQ=255)
                                                         gggTTGCTTCCAAACTCTTTTGAAGCACTGACACCGATTAAAATCGGCAAGAAAATAAATGCCGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCATTGATCATAccc                                                                                >  7:64166/1‑138 (MQ=255)
                                                           gTTGCTTCCAAACTCTTTTGAAGCACTGACACCGATTAAAATCGGCAAGAAAATAAATGCCGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCAAATGAAACGCATTGATCATAcaaaa                                                                              >  3:38736/1‑134 (MQ=255)
                                                            ttGCTTCCAAACTCTTTTGAAGCACTGACACCGATTAAAATCGGCAAGAAAATAAATGCCGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCATTGATCATACCCaata                                                                            <  3:159222/139‑1 (MQ=255)
                                                             tgctgcCAAACTCTTTTGAAGCGCTGACCCCGATTAGAATCGGCAAGAAAATGAATGCTGCACTTGAAAACATATCAAGCAGCTGCAAAAGCGCAGAATCCTTGCTCATCCAATGAAACGCATTGATCATACCCagcaa                                                                           >  2:209316/1‑135 (MQ=255)
                                                             tgattCCAAACTCTTTTGAAGCACTGACACAGATTAAAATCGTAAAGAAAATAAATGACTAACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCATTGATCATACCCaata                                                                            <  6:23507/135‑1 (MQ=255)
                                                              gCTTCCAAACTCTTTTGAAGCACTGACACCGATTAAAATCGGCAAGAAAATAAATGCCGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCATTGATCATa                                                                                   >  4:78250/1‑130 (MQ=255)
                                                              gCTTCCAAACTCTTTTGAAGCACTGACACCGATTAAAATCGGCAAGAAAATAAATGCCGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCATTGATCATa                                                                                   <  3:78250/130‑1 (MQ=255)
                                                              gCTTCCAAACTCTTTTGAAGCACTGACACCGATTAAAATCGGCAAGAAAATAAATGCCGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCATTGATCATACCCaataat                                                                          >  1:116744/1‑139 (MQ=255)
                                                                  ccAAACTCTTTTGAAGCGCTGACCCCGATTAGAATCGGCAAGAAAATGAATGCTGCACTTGAAAACATATCAAGCAGCTGCAAAAGCGCAGAATCCTTGCTCATCCAATGAAACGCATTGATCATACCCAGCAATcc                                                                        <  4:150166/137‑1 (MQ=255)
                                                                    aaaCTCTTTTGAAGCACTGACACCGATTAAAATCGGCAAGAAAATAAATGCCGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCATTGATCATACCCAATAATCCCa                                                                      >  5:117486/1‑137 (MQ=255)
                                                                    aaaCTCTTTTGAAGCACTGACACCGATTAAAATCGGCAAGAAAATAAATGCCGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCATTGATCATACCCAATAATCCCa                                                                      <  6:117486/137‑1 (MQ=255)
                                                                              gAAGCACTGACACCGATTAAAATCGGCAAGAAAATAAATGCCGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCATTGATCATACCCAATAATCCCATTAATAGGCCg                                                           <  1:242875/138‑1 (MQ=255)
                                                                              gAAGCACTGACACCGATTAAAATCGGCAAGAAAATAAATGCCGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCATTGATCATACCCAATAATCCCATTAATAGGCCg                                                           <  3:247917/138‑1 (MQ=255)
                                                                                  cACTGACACCGATTAAAATCGGCAAGAAAATAAATGCCGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCATTGATCATa                                                                                   <  6:309144/110‑1 (MQ=255)
                                                                                  cACTGACACCGATTAAAATCGGCAAGAAAATAAATGCCGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCATTGATCATa                                                                                   >  5:309144/1‑110 (MQ=255)
                                                                                  cACTGACACCGATTAAAATCGGCAAGAAAATAAATGCCGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCATTGATCATACCCAATAATCCCATTAATAGGCCGCTGGc                                                      <  1:258283/139‑1 (MQ=255)
                                                                                  cACTGACACCGATTAAAATCGGCAAGAAAATAAATGCCGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCATTGATCATACCCAATAATCCCATTAATAGGCCGCTGGc                                                      <  8:64166/139‑1 (MQ=255)
                                                                                  cACTGACACCGATTAAAATCGGCAAGAAAATAAATGCCGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCATTGATCATACCCAATAATCCCATTAATAGGCCGCTGGc                                                      <  5:41394/139‑1 (MQ=255)
                                                                                  cACTGACACCGATTAAAATCGGCAAGAAAATAAATGCCGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCATTGATCATACCCAATAATCCCATTAATAGGCCGCTGGc                                                      <  3:296811/139‑1 (MQ=255)
                                                                                    cTGACACCGATTAAAATCGGCAAGAAAATAAATGCCGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCATTGATCATACCCAATAATCCCATTAATAGGCCg                                                           <  2:274370/132‑1 (MQ=255)
                                                                                    cTGACACCGATTAAAATCGGCAAGAAAATAAATGCCGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCATTGATCATACCCAATAATCCCATTAATAGGCCg                                                           >  1:274370/1‑132 (MQ=255)
                                                                                    cTGACACCGATTAAAATCGGCAAGAAAATAAATGCCGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCATTGATCATACCCAATAATCCCATTAATAGGCCGCTGGCca                                                    >  8:330408/1‑139 (MQ=255)
                                                                                             aTTAAAATCGGCAAGAAAATAAATGCCGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAGACGCATTGATCATACCCAATAATCCCATTAATAGGCCGCTGGCCACAATAGCCg                                           >  2:288637/1‑139 (MQ=255)
                                                                                             aTTAAAATCGGCAAGAAAATAAATGCCGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCATTGATCATACCCAATAATCCCATTAATAGGCCGCTGGCCACAATAGCCg                                           >  1:74843/1‑139 (MQ=255)
                                                                                                aaaaTCGGCAAGAAAATAAATGCCGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCATTGATCATACCCAATAATCCCATTAATAGGCCGCTGGCCACAATAGCCggg                                         <  5:186850/138‑1 (MQ=255)
                                                                                                          aGAAAATAAATGCCGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCATTGATCATACCCAATAATCCCATTAATAGGCCGCTGGCCACAATAGCCGGGATGATTGGAAc                              <  1:78769/139‑1 (MQ=255)
                                                                                                          aGAAAATAAATGCCGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCATTGATCATACCCAATAATCCCATTAATAGGCCGCTGGCCACAATAGCCGGGATGATTGGAAc                              <  1:83997/139‑1 (MQ=255)
                                                                                                            aaaaTAAATGCCGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCATTGATCATACCCAATAATCCCATTAATAGGCCGCTGGCCACAATAGCCGGGATGATTGGAACaa                            <  8:180907/139‑1 (MQ=255)
                                                                                                              aataaatGCCGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCATTGATCATACCCAATAATCCCATTAATAGGCCGCTGGCCACAATAGCCGGGATGATTGGAACaaaa                          <  1:285716/139‑1 (MQ=255)
                                                                                                                      ccGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCATTGATCATa                                                                                   <  5:25396/74‑1 (MQ=255)
                                                                                                                      ccGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCATTGATCATa                                                                                   >  6:25396/1‑74 (MQ=255)
                                                                                                                      ccGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCATTGATCATACCCAATAATCCCATTAATAGGCCGCTGGCCACAATAGCCGGGATGATTGGAACAAAAATATTAg                   >  8:335765/1‑138 (MQ=255)
                                                                                                                      ccGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCATTGATCATACCCAATAATCCCATTAATAGGCCGCTGGCCACAATAGCCGGGATGATTGGAACAAAAATATAAg                   >  4:4530/1‑138 (MQ=255)
                                                                                                                      ccGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCATTGATCATACCCAATAATACCATTAATAGGCCGCTGGCCACAATAGCCGGGATGATTGGAACAAAAATATTAg                   >  7:75764/1‑138 (MQ=255)
                                                                                                                                      gTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCAt                                                                                           <  2:186385/50‑1 (MQ=255)
                                                                                                                                      gTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCAt                                                                                           >  1:186385/1‑50 (MQ=255)
                                                                                                                                      gTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCATTGATCATACCCAATAATCCCATTAATAGGCCGCTGGCCACAATAGCCGGGATGATTGGAACAAAAATATTAGAAAGCGTTTTCGCaaa   >  5:325890/1‑138 (MQ=255)
                                                                                                                                      gTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCATTGATCATACCCAATAATCCCATTAATAGGCCGCTGGCCACAATAGCCGGGATGATTGGAACAAAAATATTAGAAAGCGTTTTCGCaaa   >  3:131497/1‑138 (MQ=255)
                                                                                                                                      gTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCATTGATCATACCCAATAATCCCATTAATAGGCCGCTGGCCACAATAGCCGGGATGATTGGAACAAAAATATTAGAAAGCGTTTTCGCaaa   >  4:50986/1‑138 (MQ=255)
                                                                                                                                      gTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCATTGATCATACCCAATAATCCCATTAATAGGCCGCTGGCCACAATAGCCGGGATGATTGGAACAAAAATATTAGAAAGCGTTTTCGCAAAt  >  6:278526/1‑139 (MQ=255)
                                                                                                                                      gTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCATTGATCATACCCAATAATCCCATTAATAGGCCGCTGGCCACAATAGCCGGGATGATTGGAACAAAAATATTAGAAAGCGTTTTCGCAAAt  >  1:43970/1‑139 (MQ=255)
                                                                                                                                      gTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCATTGATCATACCCAATAATCCCATTAATAGGCCGCTGGCCACAATAGCCGGGATGATTGGAACAAAAATATTAGAAAGCGTTTTCGCAAAt  >  2:262921/1‑139 (MQ=255)
                                                                                                                                      gTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCATTGATCATACCCAATAATCCCATTAATAGGCCGCTGGCCACAATAGCCGGGATGATTGGAACAAAAATATTAGAAAGCGTTTTCGCAAAt  >  4:49236/1‑139 (MQ=255)
                                                                                                                                      |                                                                                                                                          
CATGGATTCAAAAGCGTCGGATGAATCATGATACCTCCGATGACCGCTCCCAAGTAAGGGTTGCTTCCAAACTCTTTTGAAGCACTGACACCGATTAAAATCGGCAAGAAAATAAATGCCGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCATTGATCATACCCAATAATCCCATTAATAGGCCGCTGGCCACAATAGCCGGGATGATTGGAACAAAAATATTAGAAAGCGTTTTCGCAAAT  >  CP006881/3722303‑3722575

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: