Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,098,076 C→T 5.4% A76T (GCC→ACC)  wzzB ← regulator of length of O‑antigen component of lipopolysaccharide chains

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,098,0760CT 5.4% 230.6 / 2.2 74A76T (GCC→ACC) wzzBregulator of length of O‑antigen component of lipopolysaccharide chains
Reads supporting (aligned to +/- strand):  ref base C (31/39);  new base T (2/2);  total (33/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.73e-01

AAGGTTCGATGGTAAGTTTTTCTCGTTCTTCCTGATTATCCAGCGTTTCTGCTAATGCAGAGAAGGCAGAACTGAAGCGACCAATTAACGTCTCCTGCAAATCCGATACTTTCGGTGCAGCCTGACCATAGATAACATTCATGGCATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCTGTTGACGTCCATTTCTCCTTCGCTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACAA  >  NC_000913/2097932‑2098223
                                                                                                                                                |                                                                                                                                                   
aaGGTTCGATGGTAAGTTTTTCTCGTTCTTCCTGATTATCCAGCGTTTCTGCTAATGCAGAGAAGGCAGAACTGAAGCGACCAATTAACGTCTCCTGCAAATCCGATACTTTCGGTGCAGCCTGACCATAGATAACATTCATGGCATTg                                                                                                                                                 <  1:118560/149‑1 (MQ=255)
aaGGTTCGATGGTAAGTTTTTCTCGTTCTTCCTGATTATCCAGCGTTTCTGCTAATGCAGAGAAGGCAGAACTGAAGCGACCAATTAACGTCTCCTGCAAATCCGATACTTTCGGTGCAGCCTGACCATAGATAACATTCATGGCATTg                                                                                                                                                 <  6:225832/149‑1 (MQ=255)
  ggTTCGATGGTAAGTTTTTCTCGTTCTTCCTGATTATCCAGCGTTTCTGCTAATGCAGAGAAGGCAGAACTGAAGCGACCAATTAACGTCTCCTGCAAATCCGATACTTTCGGTGCAGCCTGACCATAGATAACATTCATGGCATTGtt                                                                                                                                               <  8:266607/149‑1 (MQ=255)
  ggTTCGATGGTAAGTTTTTCTCGTTCTTCCTGATTATCCAGCGTTTCTGCTAATGCAGAGAAGGCAGAACTGAAGCGACCAATTAACGTCTCCTGCAAATCCGATACTTTCGGTGCAGCCTGACCATAGATAACATTCATGGCATTGtt                                                                                                                                               <  6:202722/149‑1 (MQ=255)
      cGATGGTAAGTTTTTCTCGTTCTTCCTGATTATCCAGCGTTTCTGCTATTGCAGAGAAGGCAGAACTGAAGCGACCAATTAACGTCTCCTGCAAATCCGATACTTTCGGTGCAGCCTGACCATAGATAACATTCATGGCATTGTtata                                                                                                                                            <  5:130015/148‑1 (MQ=255)
      cGATGGTAAGTTTTTCTCGTTCTTCCTGATTATCCAGCGTTTCTGCTAATGCAGAGAAGGCAGATCTGAAGCGACCAATTAACGTCTCCTGCAAATCCGATACTTTCGGTGCAGCCTGACCATAGATAACATTCATGGCATTGTTATAg                                                                                                                                           >  5:187911/1‑149 (MQ=255)
      cGATGGTAAGTTTTTCTCGTTCTTCCTGATTATCCAGCGTTTCTGCTAATGCAGAGAAGGCAGAACTGAAGCGACCAATTAACGTCTCCTGCAAATCCGATACTTTCGGTGCAGCCTGACCATAGATAACATTCATGGCATTGTtata                                                                                                                                            <  2:240412/148‑1 (MQ=255)
      cGATGGTAAGTTTTTCTCGTTCTTCCTGATTATCCAGCGTTTCTGCTAATGCAGAGAAGGCAGAACTGAAGCGACCAATTAACGTCTCCTGCAAATCCGATACTTTCGGTGCAGCCTGACCATAGATAACATTCATGGCATTGTtata                                                                                                                                            <  1:193352/148‑1 (MQ=255)
           gTAAGTTTTTCTCGTTCTTCCTGATTATCCAGCGTTTCTGCTAATGCAGAGAAGGCAGAACTGAAGCGACCAATTAACGTCTCCTGCAAATCCGATACTTTCGGTGCAGCCTGACCATAGATAACATTCATGGCATTGTTATAg                                                                                                                                           <  3:161919/144‑1 (MQ=255)
           gTAAGTTTTTCTCGTTCTTCCTGATTATCCAGCGTTTCTGCTAATGCAGAGAAGGCAGAACTGAAGCGACCAATTAACGTCTCCTGCAAATCCGATACTTTCGGTGCAGCCTGACCATAGATAACATTCATGGCATTGTTATAg                                                                                                                                           >  4:161919/1‑144 (MQ=255)
                  tttCTCGTTCTTCCTGATTATCCAGCGTTTCTGCTAATGCAGAGAAGGCAGAACTGAAGCGACCAATTAACGTCTCCTGCAAATCCGATACTTTCGGTGCAGCCTGACCATAGATAACATTCATGGCATTGTTATAGCCAGCAATTTg                                                                                                                                <  1:357749/148‑1 (MQ=255)
                  tttCTCGTTCTTCCTGATTATCCAGCGTTTCTGCTAATGCAGAGAAGGCAGAACTGAAGCGACCAATTAACGTCTCCTGCAAATCCGATACTTTCGGTGCAGCCTGACCATAGATAACATTCATGGCATTGTTATAGCCAGCAATTTg                                                                                                                                <  2:142034/148‑1 (MQ=255)
                    tctcGTTCTTCCTGATTATCCAGCGTTTCTGCTAATGCAGAGAAGGCAGAACTGAAGCGACCAATTAACGTCTCCTGCAAATCCGATACTTTCGGTGCAGCCTGACCATAGATAACATTCATGGCATTGTTATAGCCAGCAATTTGcc                                                                                                                              <  5:318806/148‑1 (MQ=255)
                         ttcttcCTGATTATCCAGCGTTTCTGCTAATGCACAGAAGGCAGAACTGAAGCGACCAATTAACGTCTCCTGCAAAACCGATACTTTCGGTGCAGCCTGACCATACAAAACATTAAAGGCATTGATATAGCCAGAAATTTGCCCCACGt                                                                                                                        >  1:377178/1‑149 (MQ=255)
                             tcCTGATTATCCAGCGTTTCTGCTAATGCAGAGAAGGCAGAACTGAAGCGACCAATTAACGTCTCCTGCAAATCCGATACTTTCGGTGCAGCCTGACCATAGATAACATTCATGGCATTGTTATAGCCAGCAATTTGCCCCACGTCggg                                                                                                                    >  7:299887/1‑149 (MQ=255)
                              ccTGATTATCCAGCGTTTCTGCTAATGCAGAGAAGGCAGAACTGAAGCGACCAATTAACGTCTCCTGCAAATCCGATACTTTCGGTGCAGCCTGACCATAGATAACATTCATGGCATTGTTATAGCCAGCAATTTGCCCCACGTCGGGc                                                                                                                   >  6:384425/1‑149 (MQ=255)
                                tGATTATCCAGCGTTTCTGCTAATGCAGAGAAGGCAGAACTGAAGCGACCAATTAACGTCTCCTGCAAATCCGATACTTTCGGTGCAGCCTGACCATAGATAACATTCATGGCATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCt                                                                                                                  >  6:199458/1‑148 (MQ=255)
                                 gATTATCCAGCGTTTCTGCTAATGCAGAGAAGGCAGAACTGAAGCGACCAATTAACGTCTCCTGCAAATCCGATACTTTCGGTGCAGCCTGACCATAGATAACATTCATGGCATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCt                                                                                                                  >  7:51399/1‑147 (MQ=255)
                                      tCCAGCGTTTCTGCTAATGCAGAGAAGGCAGAACTGAAGCGACCAATTAACGTCTCCTGCAAATCCGATACTTTCGGTGCAGCCTGACCATAGATAACATTCATGGCATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGAGTGa                                                                                                            <  1:190455/148‑1 (MQ=255)
                                             tttCTGCTAATGCAGAGAAGGCAGAACTGAAGCGACCAATTAACGTCTCCTGCAAATCCGATACTTTCGGTGCAGCCTGACCATAGATAACATTCATGGCATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCt                                                                                                    <  2:403970/149‑1 (MQ=255)
                                                cTGCTAATGCAGAGAAGGCAGAACTGAAGCGACCAATTAACGTCTCCTGCAAATCCGATACTTTCGGTGCAGCCTGACCATAGATAACATTCATGGCATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCTGtt                                                                                                 >  7:234114/1‑149 (MQ=255)
                                                     aaTGCAGAGAAGGCAGAACTGAAGCGACCAATTAACGTCTCCTGCAAATCCGATACTTTCGGTGCAGCCTGACCATAGATAACATTCATGGCATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCTGTTGACGt                                                                                            >  8:366643/1‑149 (MQ=255)
                                                           gagaAGGCAGAACTGAAGCGACCAATTAACGTCTCCTGCAAATCCGATACTTTCGGTGCAGACTGACCATAGATAACATTCATTGCATTGTTATAGCCAGCAATTTGACCAACGTCTGTCTGAGTGATAATTGATGTTGACGTCCAttt                                                                                      >  2:345558/1‑149 (MQ=255)
                                                                 gCAGAACTGAAGCGACCAATTAACGTCTCCTGCAAATCCGATACTTTCGGTGCAGCCTGACCATAGATAACATTCATGGCATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCTGTTGACGTCCATTTCTCCtt                                                                                <  3:316001/149‑1 (MQ=255)
                                                                 gCAGAACTGAAGCGACCAATAAACGTCTCCTCCAAATACGATACTTTCGGTGCAGCCTGACCATAGATAACATTCATGGCATTGTTATAGACAGCAATTAGCCCCACGTAGGGCTGACTGATAAATGCTGATGACGACCATTTCTCCtt                                                                                >  8:306464/1‑149 (MQ=255)
                                                                   aGAACTGAAGCGACCAATTAACGTCTCCTGCAAATCCGATACTTTCGGTGCAGCCTGACCATAGATAACATTCATGGCATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCTGTTGACGTCCATTTCTCCTTCg                                                                              <  4:34167/149‑1 (MQ=255)
                                                                   aGAACTGAAGCGACCAATTAACGTCTCCTGCAAATCCGATACTTTCGGTGCAGCCTGACCATAGATAACATTCATGGCATTGTTATAGCCAGCAATTTGCCCCACGGCGGGCTGAGTGATAATTGCTGTTGACGTCCATTTCTCCTTCg                                                                              >  6:211193/1‑149 (MQ=255)
                                                                    gatcttaaGCGACCAATTAACGTCTCCTGCAAATCCGTTACTATCGGTGCAGTCTGACCATTGATAACATACCTGGCATTGTTATAGCCAGAAATTTGACCCACGTCGGGCTGAGTGATAATTTCTGTTGACCGCCATTTCTCCTTCGc                                                                             <  6:378007/143‑1 (MQ=255)
                                                                        tGAAGCGACCAATTAACGTCTCCTGCAAATCCGATACTTTCGGTGCAGCCTGACCATAGATAACATTCATGGCATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCTGTTGACGTCCATTTCTCCTTCGCTACt                                                                         >  5:382096/1‑149 (MQ=255)
                                                                             cGACCAATTAACGTCTCCTGCAAATCCGATACTTTCGGTGCAGCCTGACCATAGATAACATTCATGGCATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCTGTTGACGTCCATTTCTCCTTCGCTACTGCCaa                                                                    <  5:384425/149‑1 (MQ=255)
                                                                              gACCAATTAACGTCTCCTGCAAATCCGATCCTTTCGGTGCAGCCTGACCATAGATAACCTTCATGGCATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCTGTTGACGTCCATTTCTCCTTCGCTACTGCCaaa                                                                   <  4:201018/149‑1 (MQ=255)
                                                                              gACCAATTAACGTCTCCTGCAAATCCGATACTTTCGGTGCAGCCTGACCATAGATAACATTCATGGCATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCTGTTGACGTCCATTTCTCCTTCGCTACTGCCaaa                                                                   >  4:173237/1‑149 (MQ=255)
                                                                                  aaTTAACGTCTCCTGCAAATCCGATACTTTCGGTGCAGCCTGACCATAGATAACATTCATGGCATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCTGTTGACGTCCATTTCTCCTTCGCTACTGCCAAATATc                                                               >  4:187274/1‑149 (MQ=255)
                                                                                  aaTTAACGTCTACTGCAAATCCGATACTTTCTGTGCAGCCTGACCATAGATAACATTCATGGCATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCTGTTGACGTCCATTTCTCCTTCGCTACTGCCAAATATc                                                               <  2:377178/149‑1 (MQ=255)
                                                                                   aTTAACGTCTCCTGCAAATCCGATACTTTCGGTGCAGCCTGACCATAGATAACATTCATGGCATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCTGTTGACGTCCATTTCTCCTTCGCTACTGCCAAATATc                                                               >  2:331963/1‑148 (MQ=255)
                                                                                       aagTCTCCTGCAAATCCGATACTTTCGGTGCAGCCTGACCATAGATATCATTCATGGCATTGTTATAGCCAGAAATTTACCCCACGTCGGGCTGAGTGATAATTGCTGTTGACGTCCATTTCTCCTTCGCTACTGCCAAATATCCAATa                                                          <  5:219387/147‑1 (MQ=255)
                                                                                       aCGTCTCCTGCAAATCCGATACTTTCGGTGCAGCCTGACCATAGATAACATTCATGGCATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCTGTTGACGTCCATTTCTCCTTCGCTACTGCCAAATATCCaa                                                            >  7:224828/1‑147 (MQ=255)
                                                                                       aCGTCTCCTGCAAATCCGATACTTTCGGTGCAGCCTGACCATAGATAACATTCATGGCATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCTGTTGACGTCCATTTCTCCTTCGCTACTGCCAAATATCCAAt                                                           <  5:199458/148‑1 (MQ=255)
                                                                                       aCGTCTCCTGCAAATCCGATACTTTCGGTGCAGCCTGACCATAGATAACATTCATGGCATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCTGTTGACGTCCATTTCTCCTTCGCTACTGCCAAATATCCAATa                                                          <  1:331963/149‑1 (MQ=255)
                                                                                       aCGTCTCCTGCAAATCCGATACTTTCGGTGCAGCCTGACCATAGATAACATTCATGGCATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCTGTTGACGTCCATTTCTCCTTCGCTACTGCCAAATATCCAATa                                                          <  1:391190/149‑1 (MQ=255)
                                                                                         gTCTCCTGCAAATCCGATACTTTCGGTGCAGCCTGACCATAGATAACATTCATGGCATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCTGTTGACGTCCATTTCTCCTTCGCTACTGCCAAATATCCAATAGc                                                        >  2:352036/1‑149 (MQ=255)
                                                                                           ctcCTGCAAATCCGATACTTTCGGTGCAGCCTGACCATAGATAACATTCATGGCATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCTGTTGACGTCCATTTCTCCTTCGCTACTGCCAAATATCCAATAGCCa                                                      <  2:27039/149‑1 (MQ=255)
                                                                                            tcCTGCAAATCCGATACTTTCGGTGCAGCCTGACCATAGATAACATTCATGGCATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCTGTTGACGTCCATTTCTCCTTCGCTACTGCCAAATATCCAATAGCCAg                                                     >  2:19369/1‑149 (MQ=255)
                                                                                               tGCAAATCCGATACTTTCGGTGCAGCCTGACCATAGATAACATTCATGGCATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCTGTTGACGTCCATTTCTCCTTCGCTACTGCCAAATATCCAATAGCCAggg                                                   >  6:338029/1‑148 (MQ=255)
                                                                                                gCAAATCCGATACTTTCGGTGCAGCCTGACCATAGATAACATTCATGGCATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCTGTTGACGTCCATTTCTCCTTCGCTACTGCCAAATATCCAATAGCCAGGGc                                                  >  4:357319/1‑148 (MQ=255)
                                                                                                gCAAATCCGATACTTTCGGTGCAGCCTGACCATAGATAACATTCATGGCATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCTGTTGACGTCCATTTCTCCTTCGCTACTGCCAAATATCCAATAGCCAGGGc                                                  >  3:69691/1‑148 (MQ=255)
                                                                                                     tCCGATACTTTCGGTGCAGCCTGACCATAGATAACATTCATGGCATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCt                                                                                                    <  1:115096/93‑1 (MQ=255)
                                                                                                     tCCGATACTTTCGGTGCAGCCTGACCATAGATAACATTCATGGCATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCt                                                                                                    >  2:115096/1‑93 (MQ=255)
                                                                                                       ccaTACTTTCGGTGCAGCCTGACCATAGATAACATTCATGGCATTGTTATAGCCAGCAATTTGCCCCACGTCGGTCTGAGTGATAATTGCTGTTGACGTCCATTTCTCCTTCGCTACTGCCAAATATCCAATAGCCAGGGCAATAGCca                                          <  4:89204/147‑1 (MQ=255)
                                                                                                       cGATACTTTCGGTGCAGCCTGACCATAGATAACATTCATGGCATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCTGTTGACGTCCATTTCTCCTTCGCTACTTCCAAATATCCAATAGCCAGGGCAATAGCca                                          >  2:363348/1‑149 (MQ=255)
                                                                                                       cGATACTTTCGGTGCAGCCTGACCATAGATAACATTCATGGCATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCTGTTGACGTCCATTTCTCCTTCGCTACTGCCAAATATCCAATAGCCAGGGCAATAGCca                                          <  8:303935/149‑1 (MQ=255)
                                                                                                        gATACTTTCGGTGCAGCCTGACCATAGATAACATTCATGGTATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCTGTTGACGTCCATTTCTCCTTCGCTACTGCCAAATATCCAATAGCCAGGGCAATAGCca                                          >  4:99861/1‑148 (MQ=255)
                                                                                                        gATACTTTCGGTGCAGCCTGACCATAGATAACATTCATGGTATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCTGTTGACGTCCATTTCTCCTACGCTACTGCCAAATATCCAATAGCCAGGGCAATAGCca                                          >  8:391317/1‑148 (MQ=255)
                                                                                                              ttCGGTGCAGCCTGACCATAGATAACATTCATGGCATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCTGTTGACGTCCATTTCTCCTTAGCTACTGCCAAATATCCAATAGCCAGGGCAATACCCACAATGAc                                   >  1:189863/1‑149 (MQ=255)
                                                                                                               tCGGTGCAGCCTGACCATAGATAACATTCATGGCATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCTGTTGACGTCCATTTCTCCTTCGCTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGAc                                   <  1:19369/148‑1 (MQ=255)
                                                                                                               tCGGTGCAGCCTGACCATAGATAACATTCATGGCATTGTTATAGCCAGCAATGTGCCCCACGTCGGGCTGAGAGATAATTGCTGTTGACGTCCATTTCTCCTTCGCTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGa                                    <  6:187911/147‑1 (MQ=255)
                                                                                                                cGGTGCAGCCTGACCATAGATAACATTCATGGTATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCTGTTGACGTCCATTTCTCCTTCGCTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACg                                  <  3:99861/148‑1 (MQ=255)
                                                                                                                cGGTGCAGCCTGACCATAGATAACATTCATGGTATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCTGTTGACGTCCATTTCTCCTACGCTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACg                                  <  7:391317/148‑1 (MQ=255)
                                                                                                                      aGCCTGACCATAGATAACATTCATGGCATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCTGTTGACGTCCATTTCTCCTTCGCTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGAAatga                           <  4:69691/149‑1 (MQ=255)
                                                                                                                         cTGACCATAGATAACATTCATGGCATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCTGTTGACGTCCATTTCTCCTTCGCTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGAAAtgatga                        <  6:382096/149‑1 (MQ=255)
                                                                                                                         cTGACCATAGATAACATTCATGGCATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCTGTTGACGTCCATTTCTCCTTCGCTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGAAAtgatga                        <  1:147303/149‑1 (MQ=255)
                                                                                                                            aCCATAGATAACATTCATGGCATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCTGTTGACGTCCATTTCTCCTTCGCTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGAAATGATGAtt                      <  8:361144/148‑1 (MQ=255)
                                                                                                                              cATAGATAACATTCATGGCATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCTGTTGACGTCCATTTCTCCTTCGCTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGAAATGATGATTGt                    >  8:194252/1‑148 (MQ=255)
                                                                                                                              cATAGATAACATTCATGGCATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCTGTTGACGTCCATTTCTCCTTCGCTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGAAATGATGATTGTc                   <  1:28928/149‑1 (MQ=255)
                                                                                                                               aTAGATAACATTCATGGCATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCTGTTGACGTCCATTTCTCCTTCGCTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGAAATGATGATTGTCa                  >  1:262049/1‑149 (MQ=255)
                                                                                                                                tAGATAACATTCATGGCATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCTGTTGACGTCCATTTCTCCTTCGCTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATTAAGGAAATGATGATAGACAt                 >  8:210334/1‑149 (MQ=255)
                                                                                                                                   aTAACATTCATGGCATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCTGTTGACGTCCATTTCTCCTTCGCTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGAAATGATGATTGTCATCtt              <  8:51399/149‑1 (MQ=255)
                                                                                                                                     aaCATTCATGGCATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCTGTTGACGTCCATTTCTCCTTCGCTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGAAATGATGATTGTCATCTTgc            >  7:162085/1‑149 (MQ=255)
                                                                                                                                      aCATTCATGGCATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCTGTTGACGTCCATTTCTCCTTCGCTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGAAATGATGATTGTCATCTTgc            <  3:16528/148‑1 (MQ=255)
                                                                                                                                        aTTCATGGCATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCTGTTGACGTCCATTTCTCCTTCGCTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGAAATGATGATTGTCATCTTgccac         >  3:283446/1‑149 (MQ=255)
                                                                                                                                            aTGGCATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCTGTTGACGTCCATTTCTCCTTCGCTACTGCCAAATATCCAATCGCCAGGGCAATAGCCACAATGACGGAAATGATGATTGTCATCTTgccacgcc      <  4:244234/148‑1 (MQ=255)
                                                                                                                                              ggCATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCTGTTGACGTCCATTTCTCCTTCGCTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGaa                               <  1:33149/121‑1 (MQ=255)
                                                                                                                                              ggCATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCTGTTGACGTCCATTTCTCCTTCGCTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGaa                               >  2:33149/1‑121 (MQ=255)
                                                                                                                                               gCATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCTGTTGACGTCCATTTCTCCTTCGCTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACaa  <  7:194252/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
AAGGTTCGATGGTAAGTTTTTCTCGTTCTTCCTGATTATCCAGCGTTTCTGCTAATGCAGAGAAGGCAGAACTGAAGCGACCAATTAACGTCTCCTGCAAATCCGATACTTTCGGTGCAGCCTGACCATAGATAACATTCATGGCATTGTTATAGCCAGCAATTTGCCCCACGTCGGGCTGAGTGATAATTGCTGTTGACGTCCATTTCTCCTTCGCTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACAA  >  NC_000913/2097932‑2098223

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