Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,970,032 (G)6→7 coding (1002/1047 nt) wzzE → Entobacterial Common Antigen (ECA) polysaccharide chain length modulation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,970,0261.G93.8% 220.7 / 5.6 64coding (996/1047 nt)wzzEEntobacterial Common Antigen (ECA) polysaccharide chain length modulation protein
Reads supporting (aligned to +/- strand):  ref base . (3/1);  new base G (27/33);  total (30/34)
Fisher's exact test for biased strand distribution p-value = 3.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCCATGTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTG‑TC‑GGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACG  >  NC_000913/3969888‑3970159
                                                                                                                                            |                                                                                                                                      
ggCCATGTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTG‑Tc                                                                                                                                         >  1:643918/1‑139 (MQ=255)
ggCCATGTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTG‑Tc                                                                                                                                         >  1:1234076/1‑139 (MQ=255)
ggCCATGTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTG‑Tc                                                                                                                                         >  2:253304/1‑139 (MQ=255)
                cTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTG‑TCGGGGGGG‑CTGATCggg                                                                                                                        >  1:187832/1‑139 (MQ=255)
                cTGAATGTTGGTCAAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCGTTG‑TCGGGGGGG‑CTGATCggg                                                                                                                        >  1:771448/1‑139 (MQ=255)
                   aaTGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTG‑TC‑GGGGGG‑CTGATCGGGGCTg                                                                                                                    >  1:777993/1‑139 (MQ=255)
                     tGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTG‑TC‑GGGGGG‑CTGATCGGGGCTGgt                                                                                                                  <  2:777993/139‑1 (MQ=255)
                       ttGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTG‑TCGGGGGGG‑CTGATCGGGGCTGgtg                                                                                                                 <  2:1234076/139‑1 (MQ=255)
                                    ggATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTG‑TCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACccg                                                                                                    >  1:516847/1‑139 (MQ=255)
                                     gATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTG‑TC‑GGGGGGCCTGATCGGGGCTGGTGTCGCATTAACccgc                                                                                                   >  2:766640/1‑139 (MQ=255)
                                         cgcgTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTG‑TCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAAcc                                                                                                      >  1:84476/1‑132 (MQ=255)
                                         cgcgTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTG‑TCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAAcc                                                                                                      <  2:84476/132‑1 (MQ=255)
                                          gcgTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTG‑TCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTg                                                                                              >  1:1069538/1‑139 (MQ=255)
                                           cgTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTG‑TCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGc                                                                                             >  2:759488/1‑139 (MQ=255)
                                           cgTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTG‑TCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGc                                                                                             >  1:1296764/1‑139 (MQ=255)
                                           cgTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTG‑TCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGc                                                                                             >  1:552209/1‑139 (MQ=255)
                                            gTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTG‑TCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCt                                                                                            <  1:253304/139‑1 (MQ=255)
                                            gTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTG‑TCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCt                                                                                            <  2:516847/139‑1 (MQ=255)
                                            gTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTG‑TCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCt                                                                                            <  1:616625/139‑1 (MQ=255)
                                            gTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTG‑TCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCt                                                                                            <  1:766640/139‑1 (MQ=255)
                                             ttttCCGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCCGATGATTATGCGGGGCATTG‑TCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTc                                                                                           <  1:819087/139‑1 (MQ=255)
                                               ttCAGACCTTTCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTG‑TCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGa                                                                                         <  2:902071/139‑1 (MQ=255)
                                                         tCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTTCTGATGATTATGCGGGGCATTG‑TCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAacaaca                                                                               >  1:1168937/1‑139 (MQ=255)
                                                         tCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTG‑TCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAacaaca                                                                               <  2:601751/139‑1 (MQ=255)
                                                         tCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTG‑TCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAacaaca                                                                               >  1:670900/1‑139 (MQ=255)
                                                                   cGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTG‑TCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGt                                                                     >  2:741723/1‑139 (MQ=255)
                                                                       cGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTG‑TCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGa                                                                                         <  2:1159230/115‑1 (MQ=255)
                                                                       cGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTG‑TCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGa                                                                                         >  1:1159230/1‑115 (MQ=255)
                                                                             aagaaCCGGTAAAACGCGATAGCCCACGCCGTGCCTTCCTGATGATTATGCGGGGCATTGTTG‑GGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAg                                                           >  2:436546/1‑139 (MQ=255)
                                                                                aaCCGGGAAAACGCGATAGCCGACGTCGTGCCTTCCTGATGATTATGCGGGGCATTG‑TCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTCCTGCGGTGAGCGCAAAGgcg                                                        <  2:771448/139‑1 (MQ=255)
                                                                                  ccGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTG‑TCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACaa                                                                                 >  1:1040399/1‑112 (MQ=255)
                                                                                  ccGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTG‑TCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACaa                                                                                 <  2:1040399/112‑1 (MQ=255)
                                                                                  ccGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTG‑TCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACActg                                                                            <  1:137566/117‑1 (MQ=255)
                                                                                  ccGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTG‑TCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACActg                                                                            >  2:137566/1‑117 (MQ=255)
                                                                                  ccGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTG‑TCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCt                                                      >  1:954555/1‑139 (MQ=255)
                                                                                   cGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTG‑TCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTc                                                     <  1:436546/139‑1 (MQ=255)
                                                                                       aaaaCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTG‑TCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTcgccg                                                 <  2:643918/139‑1 (MQ=255)
                                                                                            gcgATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTG‑TCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTAt                                            <  2:1296764/139‑1 (MQ=255)
                                                                                               aTAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTG‑TCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCg                                         >  2:1110304/1‑139 (MQ=255)
                                                                                                   cccACGTCGTGCCTTCCTGATGATTATGCGGGGCATTG‑TCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAaga                                     >  2:407989/1‑139 (MQ=255)
                                                                                                          cgtGCCTTCCTGATGATTATGCGGGGCATTG‑TCGGGGGGG‑CTGTTCGGGGTTGGTGTCGCATTAACCCGCCGTTGCTCGAAAGAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGa                              <  2:308496/139‑1 (MQ=255)
                                                                                                          cgtGCCTTCCTGATGATTATGCGGGGCATTG‑TCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGa                              <  2:954555/139‑1 (MQ=255)
                                                                                                          cgtGCCTTCCTGATGATTATGCGGGGCATTG‑TCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGa                              <  1:1110304/139‑1 (MQ=255)
                                                                                                          cgtGCCTTCCTGATGATTATGCGGGGCATTG‑TCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGa                              <  1:1138521/139‑1 (MQ=255)
                                                                                                          cgtGCCTTCCTGATGATTATGCGGGGCATTG‑TCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGa                              <  2:1069538/139‑1 (MQ=255)
                                                                                                            tGCCTTCCTGATGATTATGCGGGGCATTG‑TCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGAtg                            >  2:542137/1‑139 (MQ=255)
                                                                                                                  ccTGATGATTATGCGGGGCATTG‑TCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATaa                                                                                    <  1:797932/77‑2 (MQ=255)
                                                                                                                  ccTGATGATTATGCGGGGCATTG‑TCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATaa                                                                                    >  2:797932/1‑76 (MQ=255)
                                                                                                                      atgatTATGCGGGGCATTG‑TCGGGGGGG‑CTGATCGGGGCTGgtg                                                                                                                 >  1:1131055/1‑44 (MQ=255)
                                                                                                                      atgatTATGCGGGGCATTG‑TCGGGGGGG‑CTGATCGGGGCTGgtg                                                                                                                 <  2:1131055/44‑1 (MQ=255)
                                                                                                                      atgatTATGCGGGGCATTG‑TCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTact                  <  2:670900/139‑1 (MQ=255)
                                                                                                                        gatTATGCGGGGCATTG‑TCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactga                <  2:1168937/139‑1 (MQ=255)
                                                                                                                         atTATGCGGGGCATTG‑TCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactgac               <  2:187832/139‑1 (MQ=255)
                                                                                                                         atTATGCGGGGCATTG‑TCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactgac               <  1:542137/139‑1 (MQ=255)
                                                                                                                           tATGCGGGGCATTG‑TCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAAt                                                                                      <  2:1191405/66‑1 (MQ=255)
                                                                                                                           tATGCGGGGCATTG‑TCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAAt                                                                                      >  1:1191405/1‑66 (MQ=255)
                                                                                                                           tATGCGGGGCATTG‑TCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactgactg             >  2:272788/1‑139 (MQ=255)
                                                                                                                             tgcggggCATTG‑TCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGa                                                                                         <  2:433672/58‑1 (MQ=255)
                                                                                                                             tgcggggCATTG‑TCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGa                                                                                         >  1:433672/4‑61 (MQ=255)
                                                                                                                             tgcggggCATTG‑TCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTa           >  1:44516/4‑139 (MQ=255)
                                                                                                                                ggggCATTG‑TCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCg                               <  1:298953/116‑1 (MQ=255)
                                                                                                                                ggggCATTG‑TCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCg                               >  2:298953/1‑116 (MQ=255)
                                                                                                                                  ggCATTG‑TCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATc                                >  1:1261191/1‑113 (MQ=255)
                                                                                                                                  ggCATTG‑TCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATc                                <  2:1261191/113‑1 (MQ=255)
                                                                                                                                    cATTG‑TCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTAt                                            <  1:283076/99‑1 (MQ=255)
                                                                                                                                    cATTG‑TCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTAt                                            >  2:283076/1‑99 (MQ=255)
                                                                                                                                      ttG‑TCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTAcg  <  1:768237/139‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                      
GGCCATGTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTG‑TC‑GGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACG  >  NC_000913/3969888‑3970159

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑