Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,970,032 (G)6→7 coding (1002/1047 nt) wzzE → enterobacterial common antigen polysaccharide co‑polymerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,970,0261.G89.0% 244.6 / 22.1 73coding (996/1047 nt)wzzEenterobacterial common antigen polysaccharide co‑polymerase
Reads supporting (aligned to +/- strand):  ref base . (0/8);  new base G (49/16);  total (49/24)
Fisher's exact test for biased strand distribution p-value = 5.47e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

GGCCATGTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTC‑‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTA  >  NC_000913/3969888‑3970157
                                                                                                                                             |                                                                                                                                    
ggCCATGTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTC‑‑g                                                                                                                                    >  2:334715/1‑140 (MQ=255)
     tGTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTC‑‑gggggg                                                                                                                               >  1:1126011/1‑140 (MQ=255)
     tGTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTC‑‑gggggg                                                                                                                               >  2:1039756/1‑140 (MQ=255)
      gTTAAACCCCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTC‑‑ggggggg                                                                                                                              <  2:1536567/140‑2 (MQ=255)
      gTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTC‑‑ggggggg                                                                                                                              <  2:1335142/140‑2 (MQ=255)
       ttAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTC‑‑gggggggc                                                                                                                             <  2:197063/140‑3 (MQ=255)
           acccccTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGAt                                                                                                                          <  2:1512397/137‑1 (MQ=255)
           acacCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTAAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGAt                                                                                                                          <  1:698196/140‑1 (MQ=255)
               ccTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCggg                                                                                                                      >  2:915715/1‑140 (MQ=255)
               ccTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCggg                                                                                                                      >  2:1211752/1‑140 (MQ=255)
               ccTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTACTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCggg                                                                                                                      >  1:631363/1‑140 (MQ=255)
                cTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCgggg                                                                                                                     >  1:16175/1‑140 (MQ=255)
                 tGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGc                                                                                                                    >  1:102688/1‑140 (MQ=255)
                 tGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGc                                                                                                                    >  1:13967/1‑140 (MQ=255)
                   aaTGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTg                                                                                                                  >  2:810647/1‑140 (MQ=255)
                     tGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTTCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGgt                                                                                                                <  2:915630/140‑1 (MQ=255)
                      gTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGgtg                                                                                                               <  1:14603/140‑1 (MQ=255)
                         ggTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCg                                                                                                            >  1:1268403/1‑140 (MQ=255)
                         ggTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCg                                                                                                            >  1:432944/1‑140 (MQ=255)
                                     gATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACccgcc                                                                                                >  1:1539299/1‑140 (MQ=255)
                                         cgcgTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGttt                                                                                            >  2:1539704/1‑139 (MQ=255)
                                         cgcgTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTg                                                                                            >  1:875938/1‑140 (MQ=255)
                                          gcgTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGc                                                                                           >  2:21369/1‑140 (MQ=255)
                                          gcgTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGc                                                                                           >  1:820979/1‑140 (MQ=255)
                                          gcgTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGc                                                                                           >  2:930121/1‑140 (MQ=255)
                                           cgTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCt                                                                                          >  2:449040/1‑140 (MQ=255)
                                           cgTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCt                                                                                          >  1:1441268/1‑140 (MQ=255)
                                           cgTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCt                                                                                          >  2:1239635/1‑140 (MQ=255)
                                           cgTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCt                                                                                          >  2:1237219/1‑140 (MQ=255)
                                            tttttCAGCCCTATCGCTTTTTGCGTACGCCGGAAGAACCGGTAAACCGCGATAGCCCCCGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTc                                                                                         <  2:13967/139‑1 (MQ=255)
                                              tttCAGCCCTACCGCTTTTTGCGTCCGCCGGAAGACCCGGTAAAACGCGATGGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTGG‑‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGa                                                                                       <  2:1273086/140‑1 (MQ=255)
                                              tttCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCGCGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGa                                                                                       <  2:1438673/140‑1 (MQ=255)
                                                      cTATCGCTATTTGCGTACGCGGGAAGACCCGGTAAACCGCGATAGCCCCCGTCGTGCCTTCCTGATGATTTTGGGGGGCATTGTGG‑‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCaa                                                                               <  2:16175/140‑1 (MQ=255)
                                                      cTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCaa                                                                               >  1:1006226/1‑140 (MQ=255)
                                                       tATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAac                                                                              >  1:1153003/1‑140 (MQ=255)
                                                        aTCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAATTAGCAaca                                                                             >  2:412085/1‑140 (MQ=255)
                                                        aTCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAaca                                                                             >  2:1424820/1‑140 (MQ=255)
                                                        aTCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAaca                                                                             >  1:1356391/1‑140 (MQ=255)
                                                        aTCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAaca                                                                             >  1:790227/1‑140 (MQ=255)
                                                        aTCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAaca                                                                             >  1:172355/1‑140 (MQ=255)
                                                         tCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAacac                                                                            >  1:19447/1‑140 (MQ=255)
                                                                ttGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCg                                                                     >  1:1074362/1‑140 (MQ=255)
                                                                  gcgtcggCCGGAAGACCCGGAAAAGCGCGAGGGCCCCCGTGGTGCCTTCTTGATGATGGAGGGGGGCTTTGTCGGGGGGGTTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGt                                                                   <  2:712297/134‑1 (MQ=255)
                                                                  gCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTGG‑‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGt                                                                   <  2:1301459/140‑1 (MQ=255)
                                                                     tgcGCCGGAAGACCGGGTAAAACGCGATAGCCCTCGTCGTGCCTTCCTGATGATTATGTGGGGCTTTGTGG‑‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAg                                                                <  2:820979/138‑1 (MQ=255)
                                                                            gaagacCCGGAAAACCGCGATAGCCCGCGTCGTCCCTTCCTGATGTTTATGGGGGGCATTGGGG‑‑GGGGGGCTGATGGGGGCGGGTGTCGCATTACCCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAg                                                         <  2:471590/140‑1 (MQ=255)
                                                                            gaagaaCCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAg                                                         >  2:1104746/1‑140 (MQ=255)
                                                                             aagaaCCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAgg                                                        >  1:1549563/1‑140 (MQ=255)
                                                                             aagaaCCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAgg                                                        >  2:1479962/1‑140 (MQ=255)
                                                                                 cccGGAAAAACGCGATAGCCAACGTCGTGCCTTCCTGATGTTTTTGGGGGGTTTTGGCGGGGGGGCTGATCGGGGCGGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCt                                                    <  2:631363/139‑1 (MQ=255)
                                                                                 aCCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTGG‑‑GGGGGGTTGATCGGGGCTGGTGTCGCATTACCCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCt                                                    <  2:1222889/140‑1 (MQ=255)
                                                                                  ccGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTc                                                   >  1:1571958/1‑140 (MQ=255)
                                                                                  ccGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTc                                                   >  1:1199938/1‑140 (MQ=255)
                                                                                      aaaacCGCGATGGCCCGGGTGGGGCTTTCCTGATGATTATGGGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTACCCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTcgccg                                               <  2:790227/139‑1 (MQ=255)
                                                                                      aaaaaCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTcgccg                                               <  1:334715/139‑1 (MQ=255)
                                                                                           cgcgATGGCCCCCGTCGTGCCTTCTTGATGATTATGTGGGGCATTGTGG‑‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTAt                                          <  2:1054156/140‑1 (MQ=255)
                                                                                                           gtGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGAtg                          <  1:613813/140‑1 (MQ=255)
                                                                                                            tGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGAtgt                         >  1:1492286/1‑140 (MQ=255)
                                                                                                            tGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGAtgt                         >  2:902756/1‑140 (MQ=255)
                                                                                                             gCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGAtgtg                        <  1:1237219/140‑1 (MQ=255)
                                                                                                              ccTTCTTGATGATTATGGGGGGCTTTGTGG‑‑GGGGGGCTGATCGGGGCTGGTGTCGCATTACCCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGa                       <  2:347744/140‑1 (MQ=255)
                                                                                                              ccTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGa                       >  2:1115286/1‑140 (MQ=255)
                                                                                                               cTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGaa                      <  1:1140890/140‑1 (MQ=255)
                                                                                                                  ccTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGt                   >  2:754983/1‑140 (MQ=255)
                                                                                                                   cTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTa                  <  1:915715/140‑1 (MQ=255)
                                                                                                                   cTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTa                  <  1:930121/140‑1 (MQ=255)
                                                                                                                   cTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTa                  >  2:109913/1‑140 (MQ=255)
                                                                                                                     gatgatTAGGTGGGGCTTTGGGG‑‑GGGGGGCTGAGGGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTact                <  2:1006226/140‑1 (MQ=255)
                                                                                                                      atgttTATGTGGGGCTTGGTCGGGGGGGCTGAGCGGGGCTGGTGTCGCATTACCCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactg               <  2:1153003/140‑1 (MQ=255)
                                                                                                                       tgatTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactga              <  1:449040/140‑1 (MQ=255)
                                                                                                                        gatTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactgac             <  2:432944/140‑1 (MQ=255)
                                                                                                                           tATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGt          >  1:1332605/1‑140 (MQ=255)
                                                                                                                           tATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGt          >  1:1355646/1‑140 (MQ=255)
                                                                                                                           tATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGt          >  1:577437/1‑140 (MQ=255)
                                                                                                                           tATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGt          >  1:191142/1‑140 (MQ=255)
                                                                                                                             tgtgGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTAt        >  1:56659/1‑140 (MQ=255)
                                                                                                                             tgtgGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTAt        >  2:620394/1‑140 (MQ=255)
                                                                                                                             tgtgGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTAt        >  1:680645/1‑140 (MQ=255)
                                                                                                                             tgtgGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTAt        >  2:1527737/1‑140 (MQ=255)
                                                                                                                                ggggCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTg     >  2:226424/1‑140 (MQ=255)
                                                                                                                                   gCATTGTCGGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGt   >  2:1193855/1‑140 (MQ=255)
                                                                                                                                   gCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTa  >  2:738519/1‑140 (MQ=255)
                                                                                                                                             |                                                                                                                                    
GGCCATGTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTC‑‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTA  >  NC_000913/3969888‑3970157

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