Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 3,722,431 G→A 6.7% F152F (TTC→TTT 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,4310GA 6.7% 241.2 / 4.5 74F152F (TTC→TTTU712_19190PTS system sucrose‑specific EIIBC component
Reads supporting (aligned to +/- strand):  ref base G (28/41);  new base A (2/3);  total (30/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.82e-01

GCCAAACCCCATGGATTCAAAAGCGTCGGATGAATCATGATACCTCCGATGACCGCTCCCAAGTAAGGGTTGCTTCCAAACTCTTTTGAAGCACTGACACCGATTAAAATCGGCAAGAAAATAAATGCCGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCATTGATCATACCCAATAATCCCATTAATAGGCCGCTGGCCACAATAGCCGGGATGATTGGAACAAAAATATTAG  >  CP006881/3722294‑3722558
                                                                                                                                         |                                                                                                                               
gCCAAACCCCATGGATTCAAAAGCGTCGGATGAATCATGATACCTCCGATGACCGCTCCCAAGTAAGGGTTGCTTCCAAACTCTTTTGAAGCACTGACACCGATTAAAATCGGCAAGAAAATAAATGCCGCACTTgaga                                                                                                                                <  3:173796/139‑1 (MQ=255)
   aaaCCCCATGGATTCAAAAGCGTCGGATGAATCATGATACCTCCGATGACCGCTCCCAAGTAAGGGTTGCTTCCAAACTCTTTTGAAGCACTGACACCGATTAAAATCGGCAAGAAAATAAATGCCGCACTTGAGAACa                                                                                                                             <  3:177270/139‑1 (MQ=255)
         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                                                                           <  3:294479/139‑1 (MQ=255)
                                                     cGCTCCCAAGTAAGGGTTGCTTCCAAACTCTTTTGAAGCACTGACACCGATTAAAATCGGCAAGAAAATAAATGCCGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCa                                                                           <  2:294128/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                                                                  <  5:82105/139‑1 (MQ=255)
                                                              gTAAGGGTTGCTTCCAAACTCTTTTGAAGCACTGACACCGATTAAAATCGGCAAGAAAATAAATGCCGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCATTGATCATa                                                                  <  8:4465/139‑1 (MQ=255)
                                                              gTAAGGGTTGCTTCCAAACTCTTTTGAAGCACTGACACCGATTAAAATCGGCAAGAAAATAAATGCCGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCATTGATCATa                                                                  <  2:228706/139‑1 (MQ=255)
                                                                  gggTTGCTTCCAAACTCTTTTGAAGCACTGACACCGATTAAAATCGGCAAGAAAATAAATGCCGCACTTGAGAACATGTCGagc                                                                                                                     >  1:40020/1‑84 (MQ=255)
                                                                  gggTTGCTTCCAAACTCTTTTGAAGCACTGACACCGATTAAAATCGGCAAGAAAATAAATGCCGCACTTGAGAACATGTCGagc                                                                                                                     <  2:40020/84‑1 (MQ=255)
                                                                  gggTTGCTTCCAAACTCTTTTGAAGCACTGACACCGATTAAAATCGGCAAGAAAATAAATGCCGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCa                                                                                                      <  1:348344/99‑1 (MQ=255)
                                                                  gggTTGCTTCCAAACTCTTTTGAAGCACTGACACCGATTAAAATCGGCAAGAAAATAAATGCCGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCa                                                                                                      >  2:348344/1‑99 (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)
                                                                      tGCTTCCAAACTCTTTTGAAGCACTGACACCGATTAAAATCGGCAAGAAAATAAATGCCGCACTTgaga                                                                                                                                >  1:331965/1‑69 (MQ=255)
                                                                      tGCTTCCAAACTCTTTTGAAGCACTGACACCGATTAAAATCGGCAAGAAAATAAATGCCGCACTTgaga                                                                                                                                <  2:331965/69‑1 (MQ=255)
                                                                       gCTTCCAAACTCTTTTGAAGCACTGACACCGATTAAAATCGGCAAGAAAATAAATGCCGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCATTGATCATa                                                                  <  3:78250/130‑1 (MQ=255)
                                                                       gCTTCCAAACTCTTTTGAAGCACTGACACCGATTAAAATCGGCAAGAAAATAAATGCCGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCATTGATCATa                                                                  >  4:78250/1‑130 (MQ=255)
                                                                       gCTTCCAAACTCTTTTGAAGCACTGACACCGATTAAAATCGGCAAGAAAATAAATGCCGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCATTGATCATACCCaataat                                                         >  1:116744/1‑139 (MQ=255)
                                                                           ccAAACTCTTTTGAAGCGCTGACCCCGATTAGAATCGGCAAGAAAATGAATGCTGCACTTGAAAACATATCAAGCAGCTGCAAAAGCGCAGAATCCTTGCTCATCCAATGAAACGCATTGATCATACCCAGCAATcc                                                       <  4:150166/137‑1 (MQ=255)
                                                                             aaaCTCTTTTGAAGCACTGACACCGATTAAAATCGGCAAGAAAATAAATGCCGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCATTGATCATACCCAATAATCCCa                                                     <  6:117486/137‑1 (MQ=255)
                                                                             aaaCTCTTTTGAAGCACTGACACCGATTAAAATCGGCAAGAAAATAAATGCCGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCATTGATCATACCCAATAATCCCa                                                     >  5:117486/1‑137 (MQ=255)
                                                                                       gAAGCACTGACACCGATTAAAATCGGCAAGAAAATAAATGCCGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCATTGATCATACCCAATAATCCCATTAATAGGCCg                                          <  3:247917/138‑1 (MQ=255)
                                                                                       gAAGCACTGACACCGATTAAAATCGGCAAGAAAATAAATGCCGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCATTGATCATACCCAATAATCCCATTAATAGGCCg                                          <  1:242875/138‑1 (MQ=255)
                                                                                           cACTGACACCGATTAAAATCGGCAAGAAAATAAATGCCGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCATTGATCATa                                                                  >  5:309144/1‑110 (MQ=255)
                                                                                           cACTGACACCGATTAAAATCGGCAAGAAAATAAATGCCGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCATTGATCATa                                                                  <  6:309144/110‑1 (MQ=255)
                                                                                           cACTGACACCGATTAAAATCGGCAAGAAAATAAATGCCGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCATTGATCATACCCAATAATCCCATTAATAGGCCGCTGGc                                     <  5:41394/139‑1 (MQ=255)
                                                                                           cACTGACACCGATTAAAATCGGCAAGAAAATAAATGCCGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCATTGATCATACCCAATAATCCCATTAATAGGCCGCTGGc                                     <  8:64166/139‑1 (MQ=255)
                                                                                           cACTGACACCGATTAAAATCGGCAAGAAAATAAATGCCGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCATTGATCATACCCAATAATCCCATTAATAGGCCGCTGGc                                     <  1:258283/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:83997/139‑1 (MQ=255)
                                                                                                                   aGAAAATAAATGCCGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCATTGATCATACCCAATAATCCCATTAATAGGCCGCTGGCCACAATAGCCGGGATGATTGGAAc             <  1:78769/139‑1 (MQ=255)
                                                                                                                     aaaaTAAATGCCGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCATTGATCATACCCAATAATCCCATTAATAGGCCGCTGGCCACAATAGCCGGGATGATTGGAACaa           <  8:180907/139‑1 (MQ=255)
                                                                                                                       aataaatGCCGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCATTGATCATACCCAATAATCCCATTAATAGGCCGCTGGCCACAATAGCCGGGATGATTGGAACaaaa         <  1:285716/139‑1 (MQ=255)
                                                                                                                               ccGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCATTGATCATa                                                                  >  6:25396/1‑74 (MQ=255)
                                                                                                                               ccGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCATTGATCATa                                                                  <  5:25396/74‑1 (MQ=255)
                                                                                                                               ccGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCATTGATCATACCCAATAATCCCATTAATAGGCCGCTGGCCACAATAGCCGGGATGATTGGAACAAAAATATTAg  >  8:335765/1‑138 (MQ=255)
                                                                                                                               ccGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCATTGATCATACCCAATAATCCCATTAATAGGCCGCTGGCCACAATAGCCGGGATGATTGGAACAAAAATATAAg  >  4:4530/1‑138 (MQ=255)
                                                                                                                               ccGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCATTGATCATACCCAATAATACCATTAATAGGCCGCTGGCCACAATAGCCGGGATGATTGGAACAAAAATATTAg  >  7:75764/1‑138 (MQ=255)
                                                                                                                                         |                                                                                                                               
GCCAAACCCCATGGATTCAAAAGCGTCGGATGAATCATGATACCTCCGATGACCGCTCCCAAGTAAGGGTTGCTTCCAAACTCTTTTGAAGCACTGACACCGATTAAAATCGGCAAGAAAATAAATGCCGCACTTGAGAACATGTCGAGCAGCTGCAGCAGCGCAGAATCCTTGCTCATCCAATGAAACGCATTGATCATACCCAATAATCCCATTAATAGGCCGCTGGCCACAATAGCCGGGATGATTGGAACAAAAATATTAG  >  CP006881/3722294‑3722558

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: