Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,969,051 +G coding (21/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,969,0481.G88.9% 192.0 / 14.5 54coding (18/1047 nt)wzzEEntobacterial Common Antigen (ECA) polysaccharide chain length modulation protein
Reads supporting (aligned to +/- strand):  ref base . (4/2);  new base G (26/22);  total (30/24)
Fisher's exact test for biased strand distribution p-value = 6.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

TGCTAGCATTCTTCCTCTATGGATATTGCATTAAGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAACGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGAAACCGGCCGAAGACGCTGAAAATGAACTGGATATTCGTGGGTTGTTTCGTACCTTGTGGGCTGGGAAGCTATGGATTATTGGCATGGGGCTGGCGTTTGCGTTAATCGCGCTGGCGTATACGTTTTTTGCTCGTC  >  NC_000913/3968898‑3969187
                                                                                                                                                       |                                                                                                                                             
tGCTAGCATTCTTCCTCTATGGATATTGCATTAAGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAACGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCt                                                                                                                                                >  1:531701/1‑151 (MQ=255)
tGCTAGCATTCTTCCTCTATGGATATTGCATTAAGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAACGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCt                                                                                                                                                >  2:557607/1‑151 (MQ=255)
  cTAGCATTCTTCCTCTATGGATATTGCATTAAGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAACGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCACTTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTgg                                                                                                                                            >  1:566323/1‑151 (MQ=255)
   tAGCATTCTTCCTCTATGGATATTGCATTAAGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAACGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTg                                                                                                                                              >  1:540797/1‑149 (MQ=255)
     gCATTCTTCCTCTATGGATATTGCATTAAGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAACTCAGACTGCGTAGAAATCGTGGTGGCAGCCCAAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTgggg                                                                                                                                          <  1:478330/150‑2 (MQ=255)
       cttcttcCTCTATGGATATTGCATTAAGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAACGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTTGGgaa                                                                                                                                       >  1:144085/2‑150 (MQ=255)
       aTTCTTCCTCTATGGATATTGCATTAAGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAACGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGgaaa                                                                                                                                      <  2:606823/151‑2 (MQ=255)
       aTTCTTCCTCTATGGATATTGCATTAAGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAACGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGgaa                                                                                                                                       >  1:106087/1‑150 (MQ=255)
         tcttcCTCTATGGATATTGCATTAAGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAACGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGGAAAcc                                                                                                                                     >  1:556194/1‑151 (MQ=255)
          cttcCTCTATGGATATTGCATTAAGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAACGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGGAAACCg                                                                                                                                    >  2:557656/1‑151 (MQ=255)
           ttcCTCTATGGATATTGCATTAAGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAACGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGGAAACCgg                                                                                                                                   >  1:342706/1‑151 (MQ=255)
            tcCTCTATGGATATTGCATTAAGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAACGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGGAAACCg                                                                                                                                    <  2:545996/149‑1 (MQ=255)
                  aTGGATATTGCATTAAGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAACGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGGAAACCGGCCGaa                                                                                                                              <  2:303888/149‑1 (MQ=255)
                    ggATATTGCATTAAGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAACGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGGAAACCGGCCGAAGAc                                                                                                                           <  1:141663/150‑1 (MQ=255)
                    ggATATTGCATTAAGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAACGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGGAAACCGGCCGAAGACg                                                                                                                          >  2:533142/1‑151 (MQ=255)
                    ggATATTGCATTAAGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAACGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGGAAACCGGCCGAAGACg                                                                                                                          <  1:314750/151‑1 (MQ=255)
                                aaGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAACGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGGAAACCGGCCGAAGACGCTGAAAATGAAc                                                                                                              >  1:41169/1‑151 (MQ=255)
                                aaGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAACGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGGAAACCGGCCGAAGACGCTGAAAATGAAc                                                                                                              >  2:435155/1‑151 (MQ=255)
                                    gTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAACGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGGAAACCGGCCGAAGACGCTGAAAATGAACTgg                                                                                                           >  2:385344/1‑150 (MQ=255)
                                       ccTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAACGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGGAAACCGGCCGAAGACGCTGAAAATGAACTGGatat                                                                                                       <  2:47857/151‑1 (MQ=255)
                                         tGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAACGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGGAAACCGGCCGAAGACGCTGAAAATGAACTGGATATTc                                                                                                     >  2:44778/1‑151 (MQ=255)
                                           gAAAGCTGCTCGCTTTATTCAGCGCGTAAAACGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGGAAACCGGCCGAAGACGCTGAAAATGAACTGGATATTCg                                                                                                    >  1:618963/1‑150 (MQ=255)
                                                     cGCTTTATTAAGCGCGTAAAACGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGGAAACCGGCCGAAGACGCTGAAAATGAACTGGATATTCGTGGgttgtt                                                                                           <  2:348094/149‑1 (MQ=255)
                                                     cGCTTTATTAAGCGCGTAAAACGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGGAAACCGGCCGAAGACGCTGAAAATGAACTGGATATTCGTGGGTTGTTTc                                                                                         <  2:40774/151‑1 (MQ=255)
                                                         ttattaAGCGCGTAAAACGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGGAAACCGGCCGAAGACGCTGAAAATGAACTGGATATTCGTGGGTTGTTTCGTAc                                                                                     >  2:238118/1‑151 (MQ=255)
                                                          tattaAGCGCGTAAAACGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGGAAACCGGCCGAAGACGCTGAAAATGAACTGGATATTCGTGGGTTGTTTCGTAcc                                                                                    >  2:619276/1‑151 (MQ=255)
                                                          tattaAGCGCGTAAAACGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGGAAACCGGCCGAAGACGCTGAAAATGAACTGGATATTCGTGGGTTGTTTCGTAcc                                                                                    >  1:112795/1‑151 (MQ=255)
                                                          tattaAGCGAGTAAAACGCAGACTGCGTAGAAACCGTGTTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACGCAACCAATGCCTGAGGAAACCGGCCGAAGACGCTGAAAATGAACTGGCTATTCGTGGGTTGTTTCGTAcc                                                                                    >  1:153590/1‑151 (MQ=255)
                                                           attaAGCGCGTAAAACGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGGAAACCGGCCGAAGACGCTGAAAATGAACTGGATATTCGTGGGTTGTTTCGTACCt                                                                                   >  1:274520/1‑151 (MQ=255)
                                                              aaGCGCGTAAAACGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGGAAACCGGCCGAAGACGCTGAAAATGAACTGGATATTCGTGGGTTGTTTCGTACCTtg                                                                                 <  1:77821/150‑1 (MQ=255)
                                                              aaGCGCGTAAAACGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGGAAACCGGCCGAAGACGCTGAAAATGAACTGGATATTCGTGGGTTGTTTCGTACCTtg                                                                                 <  2:339702/150‑1 (MQ=255)
                                                               aGCGCGTAAAACGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGGAAACCGGCCGAAGACGCTGAAAATGAACTGGATATTCGTGGGTTGTTTCGTACCTtgtg                                                                               <  2:84940/151‑1 (MQ=255)
                                                                   cgTAAAACGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGGAAACCGGCCGAAGACGCTGAAAATGAACTGGATATTCGTGGGTTGTTTCGTACCTTGTGGGCt                                                                           <  2:598526/151‑1 (MQ=255)
                                                                     tAAAACGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGGAAACCGGCCGAAGACGCTGAAAATGAACTGGATATTCGTGGGTTGTTTCGTACCTTGTGGGCTgg                                                                         >  1:93518/1‑151 (MQ=255)
                                                                        aaCGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGGAAACCGGCCGAAGACGCTGAAAATGAACTGGATATTCGTGGGTTGTTTCGTACCTTGTGGGCTGGGaa                                                                      <  1:337413/151‑1 (MQ=255)
                                                                           gCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGGAAACCGGCCGAAGACGCTGAAAATGAACTGGATATTCGTGGGTTGTTTCGTACCTTGTGGGCTGGGAAGCt                                                                   >  1:116936/1‑151 (MQ=255)
                                                                           gCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGGAAACCGGCCGAAGACGCTGAAAATGAACTGGATATTCGTGGGTTGTTTCGTACCTTGTGGGCTGGGAAGCt                                                                   <  1:98872/151‑1 (MQ=255)
                                                                                 tGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGGAAACCGGCCGAAGACGCTGAAAATGAACTGGATATTCGTGGGTTGTTTCGTACCTTGTGGGCTGGGAAGCTATGGat                                                             >  2:510749/1‑151 (MQ=255)
                                                                                  gCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGGAAACCGGCCGAAGACGCTGAAAATGAACTGGATATTCGTGGGTTGTTTCGTACCTTGTGGGCTGGGAAGCTATGGatt                                                            >  2:594501/1‑151 (MQ=255)
                                                                                    gTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGGAAACCGGCCGAAGACGCTGAAAATGAACTGGATATTCGTGGGTTGTTTCGTACCTTGTGGGCTGGGAAGCTATGGattat                                                          >  2:474662/1‑151 (MQ=255)
                                                                                      aGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGGAAACCGGCCGAAGACGCTGAAAATGAACTGGATATTCGTGGGTTGTTTCGTACCTTGTGGGCTGGGAAGCTATGGATTATTg                                                        <  1:585330/151‑1 (MQ=255)
                                                                                         aaTCGTGGTGGCAGCCCCAACTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGGAAACCGGCCGAAGACGCTGAAAATGAACTGGATATTCGTGGGTTGTTTCGTACCTTGTGGGCTGGGAAGCTATGGATTATTGGCa                                                     <  1:607264/151‑1 (MQ=255)
                                                                                             gtggtgGCAGCCCCAATTTAACCAAATAGATGAGGGTGTGATGACAGAACCAATGACTGGGCAAACCGGCCGAAGATGCTGAAACATAACTGGATATTCGTGGGGTGACTCGTACATCGCGGGCTGGGACGCTATGGCTTATTGGCAt                                                    >  1:667933/1‑148 (MQ=255)
                                                                                                gtggCAGCCCCAATTTAACCAAATAAATGAAGATGTGATGACACAACCAATGCCTGGGGAAACCGGCCGAAGACGCTGAAAATGAACTGGATATTCGTGGGTTGTTTCGTACCTTGTGGGCTGGGAAGCTATGGATTATTGGCATGGGGCt                                              <  2:336704/151‑1 (MQ=255)
                                                                                                      gCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGGAAACCGGCCGAAGACGCTGAAAATGAACTGGATATTCGTGGGTTGTTTCGTACCTTGTGGGCTGGGAAGCTATGGATTATTGGCATGGGGCTGGCGtt                                        >  2:420347/1‑151 (MQ=255)
                                                                                                                aaCCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGGAAACCGGCCGAAGACGCTGAAAATGAACTGGATATTCGTGGGTTGTTTCGTACCTTGTGGGCTGGGAAGCTATGGATTATTGGCATGGGGCTGGCGTTTGCGTTAATc                              <  1:625608/151‑1 (MQ=255)
                                                                                                                aaCCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGGAAACCGGCCGAAGACGCTGAAAATGAACTGGATATTCGTGGGTTGTTTCGTACCTTGTGGGCTGGGAAGCTATGGATTATTGGCATGGGGCTGGCGTTTGCGTTAATc                              >  2:658725/1‑151 (MQ=255)
                                                                                                                  ccAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGGAAACCGGCCGAAGACGCTGAAAATGAACTGGATATTCGTGGGTTGTTTCGTACCTTGTGGGCTGGGAAGCTATGGATTATTGGCATGGGGCTGGCGTTTGCGTTAATcgc                            <  1:452811/151‑1 (MQ=255)
                                                                                                                  ccAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGGAAACCGGCCGAAGACGCTGAAAATGAACTGGATATTCGTGGGTTGTTTCGTACCTTGTGGGCTGGGAAGCTATGGATTATTGGCATGGGGCTGGCGTTTGCGTTAATcgc                            <  2:595996/151‑1 (MQ=255)
                                                                                                                    aaataaatGAGGATGTGATGACACAACCAATGCCTGGGGAAACCGGCCGAAGACGCTGAAAATGAACTGGATATTCGTGGGTTGTTTCGTACCTTGTGGGCTGGGAAGCTATGGATTATTGGCATGGGGCTGGCGTTTGCGTTAATCgcgc                          >  2:417259/1‑151 (MQ=255)
                                                                                                                      ataaatGAGGATGTGATGACACAACCAATGGCTGTGGAAACCGGCAGAAGACGCTGAAAATGAACTGGATAGGCGTGGGTTGTTTCGTGCCTTGTGTGCTGGGAAGCTATGGATTATGGGCATGGGGCTGGCGTTTGCGTTAATCGCGCTg                        <  1:35401/151‑1 (MQ=255)
                                                                                                                        aaatGAGGATGTGATGACACAACCAATGCCTGGGGAAACCGGCCGAAGACGCTGAAAATGAACTGGATATTCGTGGGTTGTTTCGTACCTTGTGGGCTTGGAAGCTATGGATTATTGGCATGGGGCTGGCGTTTGCGTTAATCGCGCTGGc                      <  1:460275/151‑1 (MQ=255)
                                                                                                                        aaatGAGGATGTGATGACACAACCAATGCCTGGGGAAACCGGCCGAAGACGCTGAAAATGAACTGGATATTCGTGGGTTGTTTCGTACCTTGTGGGCTGGGAAGCTATGGATTATTGGCATGGGGCTGGCGTTTGCGTTAATCGCGCTGGc                      >  2:493928/1‑151 (MQ=255)
                                                                                                                          atGAGGATGTGATGACACAACCAATGCCTGGGGAAACCGGCCGACGACGCTGAAAATGAACTGGATATTCGTGGGTTGTTTCGTACCTTGTGGGCTGGGAAGCTATGGATTATTGGCATGGGGCTGGCGTTTGCGTTAATCGCGCTGGCg                     >  1:412370/1‑150 (MQ=255)
                                                                                                                                aTGTGATGACACAACCAATGCCTGGGGAAACCGGCCGAAGACGCTGAAAATGAACTGGATATTCGTGGGTTGTTTCGTACCTTGTGGGCTGGGAAGCTATGGATTATTGGCATGGGGCTGGCGTTTGCGTTAATCGCGCTGGCGTATACGt              >  2:2530/1‑151 (MQ=255)
                                                                                                                                  gtgATGACACAACCAATGCCTGGGGAAACCGGCCGAAGACGCTGAAAATGAACTGGATATTCGTGGGTTGTTTCGTACCTTGTGGGCTGGGAAGCTATGGATTATTGGCATGGGGCTGGCGTTTGCGTTAATCGCGCTGGCGTATACGttt            >  2:512830/1‑151 (MQ=255)
                                                                                                                                     atgaCACAACCAATGCCTGGGGAAACCGGCCGAAGACGCTGAAAATGAACTGGATATTCGTGGGTTGTTTCGTACCTTGTGGGCTGGGAAGCTATGGATTATTGGCATGGGGCTGGCGTTTGCGTTAATCGCGCTGGCGTATACGtttttt         <  2:634012/151‑1 (MQ=255)
                                                                                                                                       gaCACAACCAATGCCTGGGGAAACCGGCCGAAGACGCTGAAAATGAACTGGATATTCGTGGGTTGTTTCGTACCTTGTGGGCTGGGAAGCTATGGATTATTGGCATGGGGCTGGCGTTTGCGTTAATCGCGCTGGCGTATACGTTTTTTg        <  1:578416/150‑1 (MQ=255)
                                                                                                                                            aaCCAATGCCTGGGGAAACCGGCCGAAGACGCTGAAAATGAACTGGATATTCGTGGGTTGTTTCGTACCTTGTGGGCTGGGAAGCTATGGATTATTGGCATGGGGCTGGCGTTTGCGTTAATCGCGCTGGCGTATAAGTTTTTTGCTCGTc  >  1:363716/1‑151 (MQ=255)
                                                                                                                                                       |                                                                                                                                             
TGCTAGCATTCTTCCTCTATGGATATTGCATTAAGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAACGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGAAACCGGCCGAAGACGCTGAAAATGAACTGGATATTCGTGGGTTGTTTCGTACCTTGTGGGCTGGGAAGCTATGGATTATTGGCATGGGGCTGGCGTTTGCGTTAATCGCGCTGGCGTATACGTTTTTTGCTCGTC  >  NC_000913/3968898‑3969187

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