Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,969,424	A→G	S132G (AGC→GGC)	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,424	0	A	G	100.0%	251.5 / NA	71	S132G (AGC→GGC)	wzzE	enterobacterial common antigen polysaccharide co‑polymerase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (32/39);  total (32/39)

TGCCGACCAACCATCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTA  >  NC_000913/3969300‑3969555                                                                                                                             |                                                                                                                                    tGCCGACCAACCATCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACGGCAAAGCCGATgcgg                                                                                                                      <  2:596462/140‑1 (MQ=255)  gCCGACCAACCATCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACGGCAAAGCCGATgcggc                                                                                                                     <  2:843542/140‑1 (MQ=255)             aTCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACGGCAAAGCCGATGCGGCGTTGCTGGATg                                                                                                          <  1:127786/140‑1 (MQ=255)               cGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACGGCAAAGCCGATGCGGCGTTGCTGGATGaa                                                                                                        <  1:189340/140‑1 (MQ=255)                     tGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACGGCAAAGCCGATGCGGCGTTGCTGGATGAAATGAtt                                                                                                  >  2:646529/1‑140 (MQ=255)                     tGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACGGCAAAGCCGATGCGGCGTTGCTGGATGAAATGAtt                                                                                                  >  2:45392/1‑140 (MQ=255)                         cGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACGGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTaaca                                                                                              <  2:304228/140‑1 (MQ=255)                           aaGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACGGCAAAGCCGATGCGGCGTTGCTGGATg                                                                                                          <  2:23984/126‑1 (MQ=255)                           aaGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACGGCAAAGCCGATGCGGCGTTGCTGGATg                                                                                                          >  1:23984/1‑126 (MQ=255)                           aaGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACGGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTaaca                                                                                              >  1:851708/1‑138 (MQ=255)                           aaGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACGGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTaaca                                                                                              <  2:851708/138‑1 (MQ=255)                            aGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACGGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAacaaca                                                                                           >  1:993788/1‑140 (MQ=255)                              ccTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACGGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATc                                                                                         <  1:14484/140‑1 (MQ=255)                              ccTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACGGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATc                                                                                         <  1:774675/140‑1 (MQ=255)                                     agagTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACGGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATc                                                                                         <  2:581775/133‑1 (MQ=255)                                     agagTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACGGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATc                                                                                         >  1:581775/1‑133 (MQ=255)                                      gagTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACGGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTAt                                                                                 >  1:857158/1‑140 (MQ=255)                                        gTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACGGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTaac                                                                                               <  2:246666/124‑1 (MQ=255)                                        gTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACGGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTaac                                                                                               >  1:246666/1‑124 (MQ=255)                                        gTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACGGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATcc                                                                               >  1:578554/1‑140 (MQ=255)                                         tttGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACGGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATcc                                                                                        >  1:941414/1‑130 (MQ=255)                                         tttGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACGGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATcc                                                                                        <  2:941414/130‑1 (MQ=255)                                         tttGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACGGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATccc                                                                              >  2:538232/1‑140 (MQ=255)                                         tttGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACGGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATccc                                                                              >  1:741667/1‑140 (MQ=255)                                         tttGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACGGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATccc                                                                              >  2:672610/1‑140 (MQ=255)                                         tttGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACGGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATccc                                                                              >  2:168672/1‑140 (MQ=255)                                           tGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACGGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCg                                                                            <  1:263502/140‑1 (MQ=255)                                           tGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACGGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCg                                                                            <  1:571580/140‑1 (MQ=255)                                                   aaCTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACGGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACttt                                                                    >  1:352568/1‑140 (MQ=255)                                                          cTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACGGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCcgcg                                                             <  1:944201/140‑1 (MQ=255)                                                          cTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACGGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCcgcg                                                             <  2:302298/140‑1 (MQ=255)                                                          cTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACGGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCcgcg                                                             >  2:718054/1‑140 (MQ=255)                                                          cTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACGGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCcgcg                                                             <  1:513904/140‑1 (MQ=255)                                                          cTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACGGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCcgcg                                                             <  1:203568/140‑1 (MQ=255)                                                            cGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACGGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCcgcgcg                                                           <  2:578554/140‑1 (MQ=255)                                                             gTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACGGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCcgcg                                                             <  1:276594/137‑1 (MQ=255)                                                             gTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACGGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCcgcg                                                             >  2:276594/1‑137 (MQ=255)                                                              tGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACGGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGt                                                         <  1:195855/140‑1 (MQ=255)                                                                  aTACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACGGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAAt                                                     >  2:738519/1‑140 (MQ=255)                                                                    aCCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACGGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGa                                                   <  1:538232/140‑1 (MQ=255)                                                                      ccGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACGGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACa                                                 >  1:297596/1‑140 (MQ=255)                                                                           gagagTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACGGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACGCGCGCGGTCAATGACAGCGTg                                            >  2:910062/1‑140 (MQ=255)                                                                              agTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACGGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAg                                         >  2:855362/1‑140 (MQ=255)                                                                               gTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACGGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGc                                        <  1:718054/140‑1 (MQ=255)                                                                                ttCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACGGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCt                                       >  2:959288/1‑140 (MQ=255)                                                                                ttCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACGGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCt                                       >  1:300938/1‑140 (MQ=255)                                                                                 tCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACGGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCtt                                      <  1:910062/140‑1 (MQ=255)                                                                                       tGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACGGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGcc                                <  2:297596/140‑1 (MQ=255)                                                                                         cAAACCGACTATTACAAACAGCGGATGGTGGGCAACGGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGa                              <  2:741667/140‑1 (MQ=255)                                                                                                 cTATTACAAACAGCGGATGGTGGGCAACGGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCgcgc                      <  2:300938/140‑1 (MQ=255)                                                                                                     tACAAACAGCGGATGGTGGGCAACGGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGa                  <  1:45392/140‑1 (MQ=255)                                                                                                          aCAGCGGATGGTGGGCAACGGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGa                  >  2:572353/1‑135 (MQ=255)                                                                                                          aCAGCGGATGGTGGGCAACGGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGa                  <  1:572353/135‑1 (MQ=255)                                                                                                          aCAGCGGATGGTGGGCAACGGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCta             <  1:881015/140‑1 (MQ=255)                                                                                                           cAGCGGATGGTGGGCAACGGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCtaa            <  1:884692/140‑1 (MQ=255)                                                                                                              cGGATGGTGGGCAACGGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCt                    >  2:618622/1‑129 (MQ=255)                                                                                                              cGGATGGTGGGCAACGGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCt                    <  1:618622/129‑1 (MQ=255)                                                                                                              cGGATGGTGGGCAACGGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCtaataa         >  2:296456/1‑140 (MQ=255)                                                                                                              cGGATGGTGGGCAACGGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCtaataa         >  1:345270/1‑140 (MQ=255)                                                                                                              cGGATGGTGGGCAACGGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCtaataa         >  1:453631/1‑140 (MQ=255)                                                                                                              cGGATGGTGGGCAACGGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCtaataa         >  2:119625/1‑140 (MQ=255)                                                                                                              cGGATGGTGGGCAACGGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCtaataa         >  1:709978/1‑140 (MQ=255)                                                                                                              cGGATGGTGGGCAACGGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCtaataa         <  1:738519/140‑1 (MQ=255)                                                                                                                gATGGTGGGCAACGGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAAcc       <  2:345270/140‑1 (MQ=255)                                                                                                                  tggtggGCAACGGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTg     <  2:709978/140‑1 (MQ=255)                                                                                                                   ggtggGCAACGGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCGGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGt    <  1:461555/140‑1 (MQ=255)                                                                                                                   ggtggGCAACGGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGt    <  1:739799/140‑1 (MQ=255)                                                                                                                      gggCAACGGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCtaataa         >  2:804789/1‑132 (MQ=255)                                                                                                                      gggCAACGGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCtaataa         <  1:804789/132‑1 (MQ=255)                                                                                                                        gCAACGGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTa  <  1:589839/137‑1 (MQ=255)                                                                                                                        gCAACGGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTa  >  2:589839/1‑137 (MQ=255)                                                                                                                             |                                                                                                                                    TGCCGACCAACCATCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACATCCAGTTTATCCCCGGAGACTTTACCCGCGCGGTCAATGACAGCGTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTA  >  NC_000913/3969300‑3969555

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: ‑