Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,544,598	A→G	N99D (AAT→GAT)	fimC →	periplasmic chaperone

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,544,598	0	A	G	98.8%	307.1 / ‑4.4	85	N99D (AAT→GAT)	fimC	periplasmic chaperone
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (39/45);  total (39/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

TATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGG  >  NC_000913/4544502‑4544691                                                                                                 |                                                                                              tATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATa                                                                                            >  2:1057323/1‑100 (MQ=255) tATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGAGGGGCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATa                                                                                            <  1:1812369/100‑1 (MQ=255)  aTTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATAc                                                                                           >  2:420726/1‑100 (MQ=255)     tAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCtt                                                                                        >  1:1214716/1‑100 (MQ=255)     tAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCtt                                                                                        >  2:1758276/1‑100 (MQ=255)         tCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGt                                                                                    <  2:1700827/100‑1 (MQ=255)          cAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTa                                                                                   <  1:319554/100‑1 (MQ=255)           aaTCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTAt                                                                                  >  1:1786060/1‑100 (MQ=255)                tGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTg                                                                             <  2:1310128/100‑1 (MQ=255)                 gggTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGa                                                                            >  2:771430/1‑100 (MQ=255)                  ggTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGAt                                                                           <  2:1214710/100‑1 (MQ=255)                  ggTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGAt                                                                           >  1:1587404/1‑100 (MQ=255)                  ggTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGAt                                                                           <  1:1620520/100‑1 (MQ=255)                  ggTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGAt                                                                           <  1:954801/100‑1 (MQ=255)                  ggTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGAt                                                                           >  1:407097/1‑100 (MQ=255)                   gTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATg                                                                          >  1:1682476/1‑100 (MQ=255)                   gTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATg                                                                          >  2:582102/1‑100 (MQ=255)                   gTGGAAAATGCCGATAGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTTTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGAt                                                                           <  2:305435/99‑1 (MQ=255)                    tGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGc                                                                         >  1:769861/1‑100 (MQ=255)                      gAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGcan                                                                       <  1:822279/100‑2 (MQ=255)                       aaaaTGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGcaac                                                                      <  1:1305913/100‑1 (MQ=255)                        aaaTGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACTTATTCTTGATGcaaca                                                                     >  1:405559/1‑100 (MQ=255)                           tGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAAt                                                                  <  1:1057328/100‑1 (MQ=255)                            gCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATa                                                                 <  2:1847044/100‑1 (MQ=255)                            gCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATa                                                                 <  1:1758290/100‑1 (MQ=255)                             ccGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATaa                                                                >  2:937749/1‑100 (MQ=255)                              cGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAAc                                                               <  1:449719/100‑1 (MQ=255)                              cGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATAAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAAc                                                               <  2:103334/100‑1 (MQ=255)                               gATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAAcc                                                              >  1:1634563/1‑100 (MQ=255)                                   gtgtAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAAtt                                                          >  2:265735/1‑100 (MQ=255)                                    tgtAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTg                                                         <  2:397951/100‑1 (MQ=255)                                    tgtAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGAAAAAAAGAAAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTg                                                         <  1:855576/100‑1 (MQ=255)                                      tAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGcc                                                       >  1:344434/1‑100 (MQ=255)                                      tAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGcc                                                       <  1:651605/100‑1 (MQ=255)                                          ggATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTAACGTAGTCTTGATGCAACAAATAACCAATTGCCACAg                                                   >  2:1375777/1‑100 (MQ=255)                                            aTGGTCGATTTATCGTGACGCCTCCTCTGTTTGAAAAAAAAAAAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAgg                                                  <  2:1363555/99‑1 (MQ=255)                                              ggTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGAcc                                               >  2:1803076/1‑100 (MQ=255)                                              ggTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGAcc                                               >  1:1509925/1‑100 (MQ=255)                                               gTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCg                                              >  1:737177/1‑100 (MQ=255)                                                 cGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCggg                                            <  2:1786045/100‑1 (MQ=255)                                                  ggtttATCGTGACACCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGa                                           >  1:1383485/3‑100 (MQ=255)                                                  gTTTTATCGTGACGCCTCCTCTGTTTGCAATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGa                                           >  2:888260/1‑100 (MQ=255)                                                     ttATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAg                                        <  2:393689/100‑1 (MQ=255)                                                     ttATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAg                                        <  1:1011228/100‑1 (MQ=255)                                                      tATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGt                                       <  2:407096/100‑1 (MQ=255)                                                      tATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGt                                       >  2:1673590/1‑100 (MQ=255)                                                       aTCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGtt                                      <  1:1235010/100‑1 (MQ=255)                                                          gTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTAt                                   >  1:1753318/1‑100 (MQ=255)                                                           tGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTAtt                                  <  1:771433/100‑1 (MQ=255)                                                              cgccgcCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACTTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCt                                <  1:1346952/99‑1 (MQ=255)                                                               gCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTgg                              <  2:1587393/100‑1 (MQ=255)                                                               gCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTgg                              >  2:778320/1‑100 (MQ=255)                                                               gCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTgg                              <  2:1753304/100‑1 (MQ=255)                                                               gCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTgg                              <  2:1579475/100‑1 (MQ=255)                                                                cctcctCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGa                             >  1:105695/1‑100 (MQ=255)                                                                cctcctCTGTTTGCGAAAAAAAGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTg                               <  2:1776080/98‑1 (MQ=255)                                                                 ctcctcTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGAt                            <  2:1634550/100‑1 (MQ=255)                                                                 ctcctcTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGAt                            >  1:1102302/1‑100 (MQ=255)                                                                 ctcctcTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGAt                            >  1:1553291/1‑100 (MQ=255)                                                                  tcctcTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATg                           <  2:1383476/100‑1 (MQ=255)                                                                   cctcTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGa                          >  2:586742/1‑100 (MQ=255)                                                                     tctGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAAc                        >  1:119077/1‑100 (MQ=255)                                                                     tctGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAAc                        >  1:802693/1‑100 (MQ=255)                                                                      ctGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACg                       <  1:420727/100‑1 (MQ=255)                                                                      ctGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACg                       <  1:1673603/100‑1 (MQ=255)                                                                      ctGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACg                       <  1:582104/100‑1 (MQ=255)                                                                      ctGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACg                       <  2:601074/100‑1 (MQ=255)                                                                       tGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAAc                        >  2:1218387/1‑98 (MQ=255)                                                                       tGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGt                      <  2:1894856/100‑1 (MQ=255)                                                                       tGTTTGCAAAGAAAGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACg                       <  2:1698229/99‑1 (MQ=255)                                                                        gTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGtt                     >  2:784059/1‑100 (MQ=255)                                                                        gTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGtt                     >  2:1495431/1‑100 (MQ=255)                                                                         tttGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTa                    >  1:878491/1‑100 (MQ=255)                                                                         tttGCGATGAAAAGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTa                    <  2:1608538/100‑1 (MQ=255)                                                                             caaTGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGc                <  1:888263/98‑1 (MQ=255)                                                                             cGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGc                <  2:769858/100‑1 (MQ=255)                                                                               aTGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGa              <  1:265736/100‑1 (MQ=255)                                                                                    gggAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCg         <  1:262059/100‑1 (MQ=255)                                                                                      gAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTc       <  1:1132260/100‑1 (MQ=255)                                                                                       aaaaaaaGAGGATACCTTACGTGTTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCa      >  2:1185816/1‑100 (MQ=255)                                                                                        aaaaaaGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCaa     <  2:405558/100‑1 (MQ=255)                                                                                         aaaaaGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAAt    >  1:1852356/1‑100 (MQ=255)                                                                                          aaaaGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAGCGTTAAAGCGATTCCGTCAATg   >  1:303772/1‑100 (MQ=255)                                                                                          aaaaGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATg   <  2:119076/100‑1 (MQ=255)                                                                                           aaaGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATgg  >  2:1979293/1‑100 (MQ=255)                                                                                                 |                                                                                              TATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGG  >  NC_000913/4544502‑4544691

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