Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,382,820	T→C	100%	Y292H (TAT→CAT)	ompG →	outer membrane porin G

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,382,820	0	T	C	98.9%	285.0 / ‑4.6	87	Y292H (TAT→CAT)	ompG	outer membrane porin G
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (44/41);  total (45/42)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.43e-01

TTTTACGGTTATGATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTAGGCCGGATA  >  NC_000913/1382724‑1382916                                                                                                 |                                                                                                 ttttACGGTTATGATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATg                                                                                               <  1:2061655/100‑1 (MQ=255) ttttACGGTTATGATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATg                                                                                               >  1:1796085/1‑100 (MQ=255)     aCGGTTATGATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAgg                                                                                           <  1:2414986/100‑1 (MQ=255)     aCGGTTATGATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAgg                                                                                           <  1:1672182/100‑1 (MQ=255)      cGGTTCTGATTTCCAGAACGGACTTTCCGTATCGCTTGCATACGCGTTTGAGTGGCAGGATCACGACGAAGTCGACAGGGATAAATTCCCTCATGCACgt                                                                                          >  2:2184175/1‑100 (MQ=255)      cGGTTATGATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGgt                                                                                          >  2:2277093/1‑100 (MQ=255)        gTTATGATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGgtgt                                                                                        >  1:985992/1‑100 (MQ=255)        gTTATGATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGgtgt                                                                                        >  2:1320466/1‑100 (MQ=255)         ttATGATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTc                                                                                       >  1:1079151/1‑100 (MQ=255)         ttATGATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGACAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTc                                                                                       >  1:1667447/1‑100 (MQ=255)         ttATGATTTCCAGAACGGACTTTCCGTTTCGCGGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTc                                                                                       <  2:1796071/100‑1 (MQ=255)          tATGATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCg                                                                                      >  1:1381923/1‑100 (MQ=255)            tGATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGc                                                                                    >  1:2146865/1‑100 (MQ=255)             gATTTCCAGAACGGACTTTCCGTTTCGCTGGTATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCg                                                                                   >  1:1732371/1‑100 (MQ=255)             gATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCg                                                                                   >  1:2170719/1‑100 (MQ=255)              aTTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGCCGAAGGCGACAGTGATAAATTCCATTATCCAGGTGTCGGctt                                                                                  <  1:1735906/100‑3 (MQ=255)              aTTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCg                                                                                   <  2:687125/99‑1 (MQ=255)              aTTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGt                                                                                  >  1:2142025/1‑100 (MQ=255)              aTTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGt                                                                                  >  1:889411/1‑100 (MQ=255)              aTTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGt                                                                                  <  2:1333404/100‑1 (MQ=255)               tttCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTa                                                                                 >  2:1953146/1‑100 (MQ=255)                ttCCAGAACAGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTaa                                                                                >  1:992297/1‑100 (MQ=255)                 tCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTaaa                                                                               >  2:2224651/1‑100 (MQ=255)                  ccAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAAt                                                                              >  1:1890295/1‑100 (MQ=255)                    aGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGAGGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAAtt                                                                             <  2:530449/99‑1 (MQ=255)                     gAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTAc                                                                           >  1:2538507/1‑100 (MQ=255)                      aaCGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTATCGGCGTAAATTACt                                                                          >  1:2520839/1‑100 (MQ=255)                        cGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCg                                                                        >  2:2412431/1‑100 (MQ=255)                        cGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCg                                                                        <  1:2462218/100‑1 (MQ=255)                               tCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGAt                                                                 <  2:1265238/100‑1 (MQ=255)                               tCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGAt                                                                 >  2:1214686/1‑100 (MQ=255)                                      cGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGc                                                          <  2:1664370/100‑1 (MQ=255)                                      cGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGAAAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGTGc                                                          >  2:1984221/1‑100 (MQ=255)                                       gCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCt                                                         >  1:166035/1‑100 (MQ=255)                                       gCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCt                                                         >  1:2472976/1‑100 (MQ=255)                                         tGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTaa                                                       <  1:1522768/100‑1 (MQ=255)                                          ggAGTACGCGTTTGAGTGGCAGGAGCACGACGGAGGCGCCAGTGATAAATCCCATCATGCAGGTGTCGGCGTAAACTACTCGTTCTGATCATGGGCTaaa                                                      <  2:303950/100‑1 (MQ=255)                                              tACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGc                                                  >  2:230445/1‑100 (MQ=255)                                               aCGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGcc                                                 <  2:1388768/100‑1 (MQ=255)                                                cgcgTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCg                                                <  1:1939803/100‑1 (MQ=255)                                                  cgTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGa                                              <  2:2072193/100‑1 (MQ=255)                                                  cgTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGa                                              >  1:337683/1‑100 (MQ=255)                                                  cgTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGACGGAGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGa                                              >  2:1909631/1‑100 (MQ=255)                                                  cgTTTAAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCCTCCTGAAGGTGTCGGAGTAAATTACTCGTTCTGTTAATGGGCTCAATTGCCg                                                <  2:2137391/98‑1 (MQ=255)                                                   gTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGAt                                             >  1:1128735/1‑100 (MQ=255)                                                    tttGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATg                                            >  1:1627324/1‑100 (MQ=255)                                                     ttGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATgc                                           <  1:1995775/100‑1 (MQ=255)                                                       gAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATgcgg                                         >  1:234307/1‑100 (MQ=255)                                                        aGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATgcggc                                        <  1:970435/100‑1 (MQ=255)                                                        aGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATgcggc                                        >  1:1703717/1‑100 (MQ=255)                                                        aGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATgcggc                                        >  1:1703777/1‑100 (MQ=255)                                                         gTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATgcggcg                                       <  2:1597628/100‑1 (MQ=255)                                                           ggCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGcgcg                                     >  1:1270748/1‑100 (MQ=255)                                                           ggCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGcgcg                                     >  1:1459457/1‑100 (MQ=255)                                                           ggCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGcgcg                                     <  1:621041/100‑1 (MQ=255)                                                           ggCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGcgcg                                     <  1:1743993/100‑1 (MQ=255)                                                           ggCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGcgcg                                     >  2:2047372/1‑100 (MQ=255)                                                            gCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGa                                    <  2:961402/100‑1 (MQ=255)                                                             cAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAg                                   <  1:2349832/100‑1 (MQ=255)                                                              aGGATCACGACGAAGGCGACGGTGATAAATTCCATGATGAAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGTCAGATGAGGCGCGaat                                  >  2:2132721/1‑98 (MQ=255)                                                              aGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGt                                  <  1:848647/100‑1 (MQ=255)                                                                gATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTAc                                <  2:856078/100‑1 (MQ=255)                                                                  tCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACtt                              >  1:2561063/1‑100 (MQ=255)                                                                   cACGACGAAGGCGACAGTGATAAATTCCAACATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACtt                              <  1:2481504/99‑1 (MQ=255)                                                                     cgacgaAGGCGACAGTGATAAATTCCATCATGCAGGTGTAGGCGTAAATTACTAGTTCTGATAATGGGATAAATTGCAGGATGCGGCGCGAGTACTTTAt                           >  2:535381/1‑100 (MQ=255)                                                                          aaGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGAt                      <  2:2262615/100‑1 (MQ=255)                                                                          aaGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGAt                      <  2:2262603/100‑1 (MQ=255)                                                                          aaGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCCGCGCGAGTACTTTATCCGAt                      <  2:1495323/100‑1 (MQ=255)                                                                           aGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATc                     >  2:722877/1‑100 (MQ=255)                                                                           aGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATc                     <  2:2275330/100‑1 (MQ=255)                                                                           aGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATc                     <  1:444547/100‑1 (MQ=255)                                                                            ggCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCt                    <  2:605369/100‑1 (MQ=255)                                                                             gCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCta                   >  1:488091/1‑100 (MQ=255)                                                                               gACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCtata                 <  1:1949934/100‑1 (MQ=255)                                                                                aCAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATaa                <  2:166033/100‑1 (MQ=255)                                                                                   gTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATg             >  2:296177/1‑100 (MQ=255)                                                                                     gATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTa           <  1:1786654/100‑1 (MQ=255)                                                                                     gATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTa           <  1:788878/100‑1 (MQ=255)                                                                                     gATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTa           >  2:564156/1‑100 (MQ=255)                                                                                      aTAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTAg          >  1:154396/1‑100 (MQ=255)                                                                                         aaTTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTAGGcc       <  2:112014/100‑1 (MQ=255)                                                                                         aaTTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTAGGcc       >  2:49656/1‑100 (MQ=255)                                                                                          aTTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTAGGcc       <  2:2087792/99‑1 (MQ=255)                                                                                          aTTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTAGGCCg      >  2:10390/1‑100 (MQ=255)                                                                                           ttCCATCCTGAAGGTGTCAGCGTAACTTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTAGGCCgg     <  1:377358/100‑1 (MQ=255)                                                                                             ccATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGTCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTAGGCCGGAt   <  2:899128/100‑1 (MQ=255)                                                                                             ccATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTAGGCCGGAt   <  1:2098316/100‑1 (MQ=255)                                                                                              catcatGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTAGGCCGGAt   <  2:2335348/99‑1 (MQ=255)                                                                                              catcatGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTAGGCCGGATa  >  1:1336113/1‑100 (MQ=255)                                                                                              catcatGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTAGGCCGGATa  <  2:381641/100‑1 (MQ=255)                                                                                                 |                                                                                                 TTTTACGGTTATGATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTAGGCCGGATA  >  NC_000913/1382724‑1382916

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