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.0%	340.6 / ‑2.4	99	N99D (AAT→GAT)	fimC	periplasmic chaperone
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base G (48/49);  total (50/49)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.47e-01

TATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATA  >  NC_000913/4544502‑4544694                                                                                                 |                                                                                                 tATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATa                                                                                               >  1:645910/1‑100 (MQ=255) tATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATa                                                                                               >  1:1348903/1‑100 (MQ=255)  aTTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATAc                                                                                              <  2:2313717/100‑1 (MQ=255)  aTTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATAc                                                                                              >  2:1669130/1‑100 (MQ=255)  aTTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATAc                                                                                              >  2:1615823/1‑100 (MQ=255)     tAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCtt                                                                                           <  1:2231526/100‑1 (MQ=255)      atctCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCGCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTa                                                                                          <  1:1649354/97‑1 (MQ=255)       aTTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTAc                                                                                         >  1:155960/1‑100 (MQ=255)        ttccgatctGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACg                                                                                        <  1:1004134/92‑1 (MQ=255)         tCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGt                                                                                       >  2:1224820/1‑100 (MQ=255)          cAATCCTGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCTATGAAGGGAAAAAAAGAGGATACCTTACGTa                                                                                      >  2:2183975/1‑100 (MQ=255)          cAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTaagat                                                                                      >  2:1649349/1‑96 (MQ=255)          cAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTa                                                                                      <  2:292533/100‑1 (MQ=255)          cAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTa                                                                                      >  2:2122233/1‑100 (MQ=255)           aaTCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTAt                                                                                     >  1:82934/1‑100 (MQ=255)           aaTCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTAt                                                                                     >  1:1187625/1‑100 (MQ=255)            aTCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCATTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTAtt                                                                                    >  2:1440218/1‑100 (MQ=255)             tCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGG‑‑AAAAAGAGGATACCTTACGTATTCtt                                                                                 >  2:1122321/1‑100 (MQ=255)              cATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAGAAAAGAGGATACCTTACGTATTCt                                                                                  <  2:1016005/100‑1 (MQ=255)               ttGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCtt                                                                                 <  2:1251542/99‑1 (MQ=255)               aTGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCtt                                                                                 <  2:155958/100‑1 (MQ=255)                tGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCtta                                                                                >  1:1251545/1‑99 (MQ=255)                tGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGAAAAAAAAAAAAAAAAAGAGGATACCTTACGTATTCTTg                                                                                <  2:82932/100‑1 (MQ=255)                 gggTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACgagatcggaa                                                                               >  2:1004132/1‑91 (MQ=255)                  ggTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATAAAGGGAAAAAAAGAGGATACCTTACGTATTCTTg                                                                                <  2:943784/98‑1 (MQ=255)                  ggTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTGATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGAt                                                                              >  1:1649904/1‑100 (MQ=255)                   gTGGAAAATGCTGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATg                                                                             <  1:1172969/100‑1 (MQ=255)                   gTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGaag                                                                             >  1:1356146/1‑98 (MQ=255)                   gTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATg                                                                             <  2:1235147/100‑1 (MQ=255)                   gTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATg                                                                             >  1:2328868/1‑100 (MQ=255)                    tGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGa                                                                               <  2:1356143/97‑1 (MQ=255)                    tGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGc                                                                            >  1:2308990/1‑100 (MQ=255)                      gAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGcaa                                                                          >  2:2148260/1‑100 (MQ=255)                       aaaaTGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGcaac                                                                         <  1:1224823/100‑1 (MQ=255)                        aaaTGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGcaaca                                                                        <  1:1384668/100‑1 (MQ=255)                           tGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAAt                                                                     <  1:2229981/100‑1 (MQ=255)                           tGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAAt                                                                     <  2:443075/100‑1 (MQ=255)                              cGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAAc                                                                  <  1:2283723/100‑1 (MQ=255)                              cGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAAc                                                                  <  2:1187622/100‑1 (MQ=255)                                 tGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCaa                                                               <  1:230026/100‑1 (MQ=255)                                     gtAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGGTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGc                                                           >  1:1653938/1‑100 (MQ=255)                                       aaaGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCca                                                         <  2:296507/100‑1 (MQ=255)                                         aGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCcaca                                                       >  2:1739833/1‑100 (MQ=255)                                         aGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCcaca                                                       <  2:1782974/100‑1 (MQ=255)                                         aGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCcaca                                                       >  2:879596/1‑100 (MQ=255)                                           gATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAgg                                                     <  2:1655428/100‑1 (MQ=255)                                             tGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGAc                                                   >  1:1910003/1‑100 (MQ=255)                                               gTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGAcc                                                  >  2:2183159/1‑99 (MQ=255)                                                 tgTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCggg                                               <  2:1197897/99‑1 (MQ=255)                                                  gTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCGATTGCCACAGGACCGGGa                                              >  2:2044647/1‑100 (MQ=255)                                                  gTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGa                                              >  1:1197900/1‑100 (MQ=255)                                                    tttATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGaaa                                            <  1:32060/100‑1 (MQ=255)                                                     ttATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAg                                           <  2:558970/100‑1 (MQ=255)                                                      tATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGt                                          >  1:2277385/1‑100 (MQ=255)                                                        tCGTGACGCCTCCTCTTTTTGCGATGAAGGAAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGttt                                        <  2:628121/100‑1 (MQ=255)                                                        tCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTTCCACAGGACCGGGAAAGttt                                        <  2:1427624/100‑1 (MQ=255)                                                        tCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGttt                                        >  2:561868/1‑100 (MQ=255)                                                        atctGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGttt                                        <  2:1481459/97‑1 (MQ=255)                                                          gTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTAt                                      >  1:1127068/1‑100 (MQ=255)                                                           tGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTAtt                                     >  2:944503/1‑100 (MQ=255)                                                            gACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTagat                                    >  1:1481463/1‑97 (MQ=255)                                                               gCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTgg                                 <  2:1649899/100‑1 (MQ=255)                                                               gCCTCCTCTGTTTGCGATGAAAGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTgg                                 <  2:506017/100‑1 (MQ=255)                                                                cctcctCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGa                                >  1:1411546/1‑100 (MQ=255)                                                                cctcctCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGa                                <  2:1127066/100‑1 (MQ=255)                                                                cctcctCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGa                                >  1:1988771/1‑100 (MQ=255)                                                                 ctcctcTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGAt                               >  1:608416/1‑100 (MQ=255)                                                                 ctcctcTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGAt                               <  1:170008/100‑1 (MQ=255)                                                                 ctcctcTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGAt                               <  1:1800463/100‑1 (MQ=255)                                                                 ctcctcTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGAt                               >  1:1063421/1‑100 (MQ=255)                                                                     tctGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAAc                           >  2:954963/1‑100 (MQ=255)                                                                      ctGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACg                          >  2:65389/1‑100 (MQ=255)                                                                      ctGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACg                          <  1:1440222/100‑1 (MQ=255)                                                                      ctGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACg                          <  1:1669135/100‑1 (MQ=255)                                                                      ctGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACg                          <  2:394285/100‑1 (MQ=255)                                                                      ctGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACg                          <  1:225539/100‑1 (MQ=255)                                                                      ctGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACg                          <  2:1056892/100‑1 (MQ=255)                                                                        gTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGtt                        <  2:2328859/100‑1 (MQ=255)                                                                        gTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGtt                        <  2:1771603/100‑1 (MQ=255)                                                                            gCGAAGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTTATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTGATTCTGGATGAACGTTaaa                     >  2:192338/1‑99 (MQ=255)                                                                              gATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCg                  >  1:2065945/1‑100 (MQ=255)                                                                                  aaGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTc              <  1:1243349/100‑1 (MQ=255)                                                                                  aaGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTc              >  1:1161525/1‑100 (MQ=255)                                                                                  aaGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAGCGTTAAAGCGATTc              >  2:1593890/1‑100 (MQ=255)                                                                                   aGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTcc             >  1:504242/1‑100 (MQ=255)                                                                                    gggAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCg            >  1:1408185/1‑100 (MQ=255)                                                                                      gAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTc          >  1:1210201/1‑100 (MQ=255)                                                                                      gAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTc          <  1:1054082/100‑1 (MQ=255)                                                                                      gAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTc          >  2:636832/1‑100 (MQ=255)                                                                                      gAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTc          <  2:608415/100‑1 (MQ=255)                                                                                      gAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTc          <  1:954965/100‑1 (MQ=255)                                                                                       aaaaaaaGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCa         >  1:1242632/1‑100 (MQ=255)                                                                                        aaaaaaGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCaa        >  2:517654/1‑100 (MQ=255)                                                                                        aaaaaaGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCaa        <  2:2065938/100‑1 (MQ=255)                                                                                         aaaaaGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAAt       >  1:252274/1‑100 (MQ=255)                                                                                         aaaaaGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAAt       >  2:1833146/1‑100 (MQ=255)                                                                                          aaaaGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATg      <  1:2072190/100‑1 (MQ=255)                                                                                           aaaGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATgg     <  2:1348900/100‑1 (MQ=255)                                                                                            aaGAGGATACCTTACGTATTCTGGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGa    <  2:473113/100‑1 (MQ=255)                                                                                              gaggATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATa  <  1:1615828/100‑1 (MQ=255)                                                                                                 |                                                                                                 TATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATA  >  NC_000913/4544502‑4544694

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