Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,969,767	T→C	V246A (GTA→GCA)	wzzE →	enterobacterial common antigen polysaccharide co‑polymerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,969,767	0	T	C	100.0%	224.8 / NA	63	V246A (GTA→GCA)	wzzE	enterobacterial common antigen polysaccharide co‑polymerase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (25/38);  total (25/38)

CGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGTACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCGTCCAATGCTTCAGGCTCGACTGGAAAATTTACAGGCCGTCGGTCCGGCCTTTGATCTCGACTATGATCAGAATCGGGCCATGTTAAACA  >  NC_000913/3969631‑3969901                                                                                                                                         |                                                                                                                                       cgTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGCAcc                                                                                                                                     >  2:628452/1‑140 (MQ=255)     ccATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGCACCTGcc                                                                                                                                 <  1:1070465/140‑1 (MQ=255)     ccATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGCACCTGcc                                                                                                                                 <  1:628452/140‑1 (MQ=255)     ccATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGCACCTGcc                                                                                                                                 <  1:36499/140‑1 (MQ=255)      cATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGCACCTGCCg                                                                                                                                <  1:841524/140‑1 (MQ=255)      cATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGCACCTGCCg                                                                                                                                <  1:597677/140‑1 (MQ=255)         ccAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGCAc                                                                                                                                      >  1:226375/1‑131 (MQ=255)         ccAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGCAc                                                                                                                                      <  2:226375/131‑1 (MQ=255)               gAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGCACCTGCCGAGGAATTAc                                                                                                                       <  1:414250/140‑1 (MQ=255)                aaaGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGCACCTGCCGAGGAATTAcc                                                                                                                      <  2:166769/140‑1 (MQ=255)                     tCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGCACCTGCCGAGGAATTACCTGAtt                                                                                                                 <  1:67954/140‑1 (MQ=255)                     tCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGCACCTGCCGAGGAATTACCTGAtt                                                                                                                 <  1:992473/140‑1 (MQ=255)                          tGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGCACCTGCCGAGGAATTACCTGATTCAGaa                                                                                                            <  1:1017205/140‑1 (MQ=255)                           gAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGCACCTGCCGAGGAATTACCTGATTCAGaaa                                                                                                           <  2:752439/140‑1 (MQ=255)                            aaGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGCACCTGCCGAGGAATTACCTGATTCAGAAAt                                                                                                          >  2:1044684/1‑140 (MQ=255)                               cGTCAGGAAGAGGTGGCGAAAGCCATCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGCACCTGCCGAGGAATTAcc                                                                                                                      <  1:361196/125‑1 (MQ=255)                               cGTCAGGAAGAGGTGGCGAAAGCCATCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGCACCTGCCGAGGAATTAcc                                                                                                                      >  2:361196/1‑125 (MQ=255)                                   aGGAAGAGGTGGCGAAAGCCATCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGCACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTg                                                                                                   >  1:236141/1‑140 (MQ=255)                                        gagGTGGCGAAAGCCATCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGCACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTgg                                                                                              <  1:134790/140‑1 (MQ=255)                                         agGTGGCGAAAGCCATCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGCACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTggg                                                                                             >  1:487457/1‑140 (MQ=255)                                         agGTGGCGAAAGCCATCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGCACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTggg                                                                                             >  1:338141/1‑140 (MQ=255)                                           gTGGCGAAAGCCATCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGCACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCg                                                                                           <  1:669591/140‑1 (MQ=255)                                             ggCGAAAGCCATCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGCACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCGTc                                                                                         <  1:23015/140‑1 (MQ=255)                                               cGAAAGCCATCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGCACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCGTCCa                                                                                       >  2:597032/1‑140 (MQ=255)                                               cGAAAGCCATCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGCACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCGTCCa                                                                                       >  2:798819/1‑140 (MQ=255)                                                   aGCCATCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGCACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTgg                                                                                              >  2:410043/1‑129 (MQ=255)                                                   aGCCATCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGCACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTgg                                                                                              <  1:410043/129‑1 (MQ=255)                                                     ccATCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGCACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCGTCCAATGCtt                                                                                 >  1:139136/1‑140 (MQ=255)                                                        tCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGCACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCGTCCAATGCTTCAg                                                                              >  2:981366/1‑140 (MQ=255)                                                           aCGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGCACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCGTCCAATGCTTCAGGCt                                                                           <  1:981366/140‑1 (MQ=255)                                                            cGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGCACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCGTCCAATGCTTCAGGCTc                                                                          <  1:798819/140‑1 (MQ=255)                                                                cgccgGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGCACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCGTCCAATGCTTCAGGCTCGACt                                                                      <  2:601443/140‑1 (MQ=255)                                                                cgccgGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGCACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCGTCCAATGCTTCAGGCTCGACt                                                                      <  2:236141/140‑1 (MQ=255)                                                                  ccgGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGCACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCGTCCAATGCTTCAGGCTCGACTgg                                                                    >  1:266604/1‑140 (MQ=255)                                                                  ccgGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGCACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCGTCCAATGCTTCAGGCTCGACTgg                                                                    >  2:48018/1‑140 (MQ=255)                                                                  ccgGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGCACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCGTCCAATGCTTCAGGCTCGACTgg                                                                    >  1:127441/1‑140 (MQ=255)                                                                            aGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGCACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCGTCCAATGCTTCAGGCTCGACTGGAAAATTTACa                                                          <  2:211490/140‑1 (MQ=255)                                                                             gCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGCACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCGTCCAATGCTTCAGGCTCGACTGGTAAATTTACAg                                                         >  1:841325/1‑140 (MQ=255)                                                                                 tGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGCACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCGTCCAATGCTTCAGGCTCGACTGGAAAATTTACAGGCCg                                                     >  2:691210/1‑140 (MQ=255)                                                                                       ggCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGCACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCGTCCAATGCTTCAGGCTCGACTGGAAAATTTACAGGCCGTCGGTc                                               <  2:487457/140‑1 (MQ=255)                                                                                       ggCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGCACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCGTCCAATGCTTCAGGCTCGACTGGAAAATTTACAGGCCGTCGGTc                                               <  1:913937/140‑1 (MQ=255)                                                                                        gcgcTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGCACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCGTCCAATGCTTCAGGCTCGACTGGAAAATTTACAGGCCGTCGGTcc                                              >  1:980438/1‑140 (MQ=255)                                                                                        gcgcTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGCACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCGTCCAATGCTTCAGGCTCGACTGGAAAATTTACAGGCCGTCGGTcc                                              <  1:438303/140‑1 (MQ=255)                                                                                         cgcTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGCACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCGTCCAATGCTTCAGGCTCGACTGGAAAATTTACAGGCCGTCGGTCCg                                             <  1:48018/140‑1 (MQ=255)                                                                                              aaaaTTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGCACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCGTCCAATGCTTCAGGCTCGACTGGAAAATTTACAGGCCGTCGGTCCGGCCtt                                        <  2:841325/140‑1 (MQ=255)                                                                                                aaTTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGCACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCGTCCAATGCTTCAGGCTCGACTGGAAAATTTACAGGCCGTCGGTCCGGCCTTTg                                      <  2:266604/140‑1 (MQ=255)                                                                                                    gCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGCACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCGTCCAATGCTTCAGGCTCGACTGGAAAATTTACAGGCCGTCGGTCCGGCCTTTGAtct                                  <  2:720989/140‑1 (MQ=255)                                                                                                        agcagcATAATATTTCGCGCAGTGCGACAGATGCACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCGTCCAATGCTTCAGGCTCGACTGGAAAATTTACAGGCCGTCGGTCCGGCCTTTGATCTCGAc                              >  1:71137/1‑140 (MQ=255)                                                                                                         gcagcaTAATATTTCGCGCAGTGCGACAGATGCACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCGTCCAATGCTTCAGGCTCGACTGGAAAATTTACAGGCCGTCGGTCCGGCCTTTGATCTCGACt                             <  2:127441/140‑1 (MQ=255)                                                                                                          cagcaTAATATTTCGCGCAGTGCGACAGATGCACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCGTCCAATGCTTCAGGCTCGACTGGGAAATTTACAGGCCGTCGGTCCGGCCTTTGATCTCGACTa                            >  2:613238/1‑140 (MQ=255)                                                                                                             caTAATATTTCGCGCAGTGCGACAGATGCACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCGTCCAATGCTTCAGGCTCGACTGGAAAATTTACAGGCCGTCGGTCCGGCCTTTGATCTCGACTATGa                         >  1:118751/1‑140 (MQ=255)                                                                                                              ataataTTTCGCGCAGTGCGACAGATGCACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCGTCCAATGCTTCAGGCTCGACTGGAAAATTTACAGGCCGTCGGTCCGGCCTTTGATCTCGACTATGAt                        <  1:667701/140‑1 (MQ=255)                                                                                                              ataataTTTCGCGCAGTGCGACAGATGCACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCGTCCAATGCTTCAGGCTCGACTGGAAAATTTACAGGCCGTCGGTCCGGCCTTTGATCTCGACTATGAt                        <  1:307715/140‑1 (MQ=255)                                                                                                                    ttccgcgCTGTGGGACAGATGCCCCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCGTCCAATGCTTCAGGCTCGACTGGAAAATTTACAGGCCGTCGGTCCGGCCTTTGATCTCGACTATGATCAGAAt                  <  2:537335/137‑1 (MQ=255)                                                                                                                       cgcgCAGTGCGACAGATGCACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCGTCCAATGCTTCAGGCTCGACTGGAAAATTTACAGGCCGTCGGTCCGGCCTTTGATCTCGACTATGATCAGAATCgg               >  2:38577/1‑140 (MQ=255)                                                                                                                       cgcgCAGTGCGACAGATGCACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCGTCCAATGCTTCAGGCTCGACTGGAAAATTTACAGGCCGTCGGTCCGGCCTTTGATCTCGACTATGATCAGAATCgg               >  2:299149/1‑140 (MQ=255)                                                                                                                        gcgcAGTGCGACAGATGCACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCGTCCAATGCTTCAGGCTCGACTGGGAAATTTACAGGCCGTCGGTCCGGCCTTTGATCTCGACTATGATCAGAATCggg              <  1:613238/140‑1 (MQ=255)                                                                                                                        gcgcAGTGCGACAGATGCACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCGTCCAATGCTTCAGGCTCGACTGGAAAATTTACAGGCCGTCGGTCCGGCCTTTGATCTCGACTATGATCAGAATCggg              <  1:691210/140‑1 (MQ=255)                                                                                                                        gcgcAGTGCGACAGATGCACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCGTCCAATGCTTCAGGCTCGACTGGAAAATTTACAGGCCGTCGGTCCGGCCTTTGATCTCGACTATGATCAGAATCggg              <  1:299149/140‑1 (MQ=255)                                                                                                                        ccggaGTGAGACAGATGCACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCGTGCAATGCTTCAGGCTCGACTGGAAAATTTACAGGCCGTCGGTCCGGCCTTTGATCTCGACTATGATCAGAATCggg              <  2:139136/136‑1 (MQ=255)                                                                                                                          gcAGTGCGACAGATGCACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCGTCCAATGCTTCAGGCTCGACTGGAAAATTTACAGGCCGTCGGTCCGGCCTTTGATCTCGACTATGATCAGAATCGGGcc            <  1:597032/140‑1 (MQ=255)                                                                                                                            aGTGCGACAGATGCACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCGTCCAATGCTTCAGGCTCGACTGGAAAATTTACAGGCCGTCGGTCCGGCCTTTGATCTCGACTATGAt                        <  1:518388/126‑1 (MQ=255)                                                                                                                            aGTGCGACAGATGCACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCGTCCAATGCTTCAGGCTCGACTGGAAAATTTACAGGCCGTCGGTCCGGCCTTTGATCTCGACTATGAt                        >  2:518388/1‑126 (MQ=255)                                                                                                                            aGTGCGACAGATGCACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCGTCCAATGCTTCAGGCTCGACTGGAAAATTTACAGGCCGTCGGTCCGGCCTTTGATCTCGACTATGATCAGAATCGGGCCAt          >  2:838830/1‑140 (MQ=255)                                                                                                                                    aGATGCACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCGTCCAATGCTTCAGGCTCGACTGGAAAATTTACAGGCCGTCGGTCCGGCCTTTGATCTCGACTATGATCAGAATCGGGCCATGTTAAaca  <  2:843191/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                                       CGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGTACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCGTCCAATGCTTCAGGCTCGACTGGAAAATTTACAGGCCGTCGGTCCGGCCTTTGATCTCGACTATGATCAGAATCGGGCCATGTTAAACA  >  NC_000913/3969631‑3969901

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