Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	238,795	T→C	30.7%	pseudogene (308/357 nt)	yafU ←	pseudogene

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	238,795	0	T	C	30.7%	108.2 / 67.9	88	pseudogene (308/357 nt)	yafU	pseudogene
Reads supporting (aligned to +/- strand):  ref base T (29/32);  new base C (15/12);  total (44/44)
Fisher's exact test for biased strand distribution p-value = 6.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACATCGAAAGAAACTCTTTTAATTTAACATGTAAACGCATGGTTAATCCTCATATCACGGGTGGAGTGTTAAGAACATACATAAATGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCTGTTGCCGTTTTAGTCCATCTCTAATTGCATATTTTAATTTTTCTGATAAATGGCATTGAGCATCGATTTCATTTA  >  NC_000913/238697‑238894                                                                                                   |                                                                                                    aCATCGAAAGAAACTCTTTTAATTTAACATGTAAACGCATGGTTAATCCTCATATCACGGGTGGAGTGTTAAGAACATACATAAATGGAGTCATGTTTTc                                                                                                    >  1:2005236/1‑100 (MQ=255)  cATCGAAAGAAACTCTTTTAATTTAACATGTAAACGCATGGTTAATCCTCATATCACGGGTGGAGTGTTAAGAACATACATAAATGGAGTCATGTTTTcc                                                                                                   <  2:1686490/100‑1 (MQ=255)   aTCGAAAGAAACTCTTTTAATTTAACATGTAAACGCATGGTTAATCCTCATATCACGGGTGGAGTGTTAAGAACATACATAAATGGAGTCATGTTTTccc                                                                                                  <  1:863367/100‑1 (MQ=255)    tCGAAAGAAACTCTTTTAATTTAACATGTAAACGCATGGTTAATCCTCATATCACGGGTGGAGTGTTAAGAACATACATAAATGGAGTCATGTTTTCCCt                                                                                                 <  2:2193716/100‑1 (MQ=255)    tCGAAAGAAACTCTTTTAATTTAACATGTAAACGCATGGTTAATCCTCATATCACGGGTGGAGTGTTAAGAACATACATAAATGGAGTCATGTTTTCCCt                                                                                                 >  2:1706212/1‑100 (MQ=255)     cGAAAGAAACTCTTTTAATTTAACATGTAAACGCATGGTTAATCCTCATATCACGGGTGGAGTGTTAAGAACATACATAAATGGAGTCATGTTTTCCCt                                                                                                 >  1:1950183/1‑99 (MQ=255)      gaaagaaaCTCTTTTAATTTAACATGTAAACGCATGGTTAATCCTCATATCACGGGTGGAGTGTTAAGAACATACATAAATGGATTCATGTTTTCCCttt                                                                                               >  1:1128385/1‑100 (MQ=255)      gaaagaaaCTCTTTTAATTTAACATGAAAACGCATGGTTAATCCTCATATCACGGGTGGAGTGTTAAGAACATACATAAATGGAGTCATGTTTTCCCttt                                                                                               >  1:1366961/1‑100 (MQ=255)          gaaaCTCTTTTAATTTAACATGTAAACGCATGGTTAATCCTCATATCACGGGTGGAGTGTTAAGAACATACATAAATGGAGTCATGTTTTCCCTTTTcc                                                                                            <  1:629249/99‑1 (MQ=255)           aaaCTCTTTTAATTTAACATGTAAACGCATTGTTAATCCTCATATCACGGGTGGAGTGTTAAGAACATACATAAATGGAGTCATGTTTTCCCTTTTCCa                                                                                           <  2:381942/99‑1 (MQ=255)              ctctTTTAATTTAACATGTAAACGCATGGTTAATCCTCATATCACGGGTGGAGTGTTAAGAACATACATAAATGGAGTCATGTTTTCCCTTTTCCAttt                                                                                        <  2:5755/99‑1 (MQ=255)                ctTTTAATTTAACATGTAAACGCATGGTTCATCCTCATATCACGGGTGGAGTGTTAAGAACATACATAAATGGAGTCATGTTTCCCCTTTTCCATTTAt                                                                                      >  1:650383/1‑99 (MQ=255)                ctTTTAATTTAACATGTAAACGCATGGTTAATCCTCATATCACGGGTGGAGTGTTAAGAACATACATAAATGGAGTCATGTTTTCCCTTTTCCATTTATc                                                                                     <  2:1985446/100‑1 (MQ=255)                ctTTTAATTTAACATGTAAACGCATGGTTAATCCTCATATCACGGGTGGAGTGTTAAGAACATACATAAATGGAGTCATGTTTCCCCTTTTCCATTTATc                                                                                     >  1:1003656/1‑100 (MQ=255)                  tttaatttaaCATGTAAACGCATGGTTAATCCTCATATCACGGGTGGAGTGTTAAGAACATACATAAATGGAGTCATGTTTTCCCTTTTCCATTTATCaa                                                                                   <  1:1973750/100‑1 (MQ=255)                   ttaatttaaCATGTAAACGCATGGTTAATCCTCATATCACGGGTGGAGTGTTAAGAACATACATAAATGGAGTCATGTTTTCCCTTTTCCATTTATCaa                                                                                   <  2:1508317/99‑1 (MQ=255)                    taatttaaCATGTAAACGCATGGTTAATCCTCATATCACGGGTGGAGTGTTAAGAACATACATAAATGGAGTCATGTTTCCCCTTTTCCATTTATCAAGt                                                                                 >  2:284465/1‑100 (MQ=255)                         taaCATGTAAACGCATGGTTAATCCTCATATCACGGGTGGAGTGTTAAGAACATACATAAATGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCTg                                                                            <  1:905775/100‑1 (MQ=255)                          aaCATGTAAACGCATGGTTAATCCTCATATCTCGGGTGGAGTGTTAAGAACATACATAAATGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCTGt                                                                           <  2:1712820/100‑1 (MQ=255)                          aaCATGTAAACGCATGGTTAATCCTCATATCACGGGTGGAGTGTTAAGAACATACATAAATGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCTGt                                                                           <  1:250944/100‑1 (MQ=255)                          aaCATGTAAACGCATGGTTAATCCTCATATCACGGGTGGAGTGTTAAGAACATACATAAATGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCTGt                                                                           >  2:1632990/1‑100 (MQ=255)                          aaCATGTAAACGCATGGTTAATCCTCATATCACGGGTGGAGTGTTAAGAACATACATAAATGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCTGt                                                                           <  2:1089627/100‑1 (MQ=255)                          aaCATGTAAACGCATGGTTAATCCTCATATCACGGGTGGAGTGTTAAGAACATACATAAATGGAGTCATGTTTCCCCTTTTCCATTTATCAAGTTCCTGt                                                                           <  2:768160/100‑1 (MQ=255)                             aTGTAAACGCATGGTTAATCCTCATATCACGGGTCGAGTGTTAAGAACATACATAAATGGAGTCATGTTTCCCCTTTTCCATTTATCAAGTTCCTGTTGc                                                                        <  2:1178861/100‑1 (MQ=255)                              tGTAAACGCATGGTTAATCCTCATATCACGGGTGGAGTGTTAAGAACATACATAAATGGAGGCATGTTTCCTCTTTTCCATTTATCAAGTTACTGTTGc                                                                        >  1:1795270/1‑99 (MQ=255)                                tAAACGCATGGTTAATCCTCATATCACGGGTGGAGTGTTAAGAACATACATAAATGGAGTCATGTTTCCCCTTTTCCATTTATCAAGTTCCTGTTGccat                                                                     <  1:1090771/100‑3 (MQ=255)                                tAAACGCATGGTTAATCCTCATATCACGGGTGGAGTGTTAAGAACATACATAAATGGAGTCATGTTTCCCCTTTTCCATTTATCAAGTTCCTGTTGCCGt                                                                     >  2:1099643/1‑100 (MQ=255)                                tAAACGCATGGTTAATAATCATATCAAGGGTGGAGTGTTAAGAAAATAAATAAATGGATTTATGTTTTCCATTTTACAATTATAAAATTTCtttt                                                                          >  2:96781/1‑95 (MQ=255)                                   ccGCATGGTTAACCCTCATATCACGCGTGGAGTGTTAAGAACATACACAAACGCAGCCATCTTTTCCCTTTTCCATTTATCAAGTTCCTGCTGCCGtttt                                                                  >  1:2238503/2‑100 (MQ=255)                                   aCGCATGGTTAATCCTCATATCACGGGTGGAGTGTTAAGAACATACATAAATGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCTGTTGCCGtttt                                                                  <  2:1852724/100‑1 (MQ=255)                                   aCGCATGGTTAATCCTCATATCACGGGTGGAGTGTTAAGAACATACATAAATGGAGTCATGTTTCCCCTTTTCCATTTATCAAGTTCCTGTTGCCGtttt                                                                  <  1:882747/100‑1 (MQ=255)                                    cGCATGGTTAATCCTCATATCACGGGTGGAGTGTTAAGAACATACATAAATGGAGTCATGTTTCCCCTTTTCCATTTATCAAGTTCCTGTTGCCGTTTTa                                                                 >  1:167196/1‑100 (MQ=255)                                     gCATGGTTAATCCTCATATCACGGGTGGAGTGTTAAGAACATACATAAATGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCTGTTGCCGTTTTAg                                                                >  2:748688/1‑100 (MQ=255)                                     gCATGGTTAATCCTCATATCACGGGTGGAGTGTTAAGAACATACATAAATGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCTGTTGCCGTTTTAg                                                                <  1:519940/100‑1 (MQ=255)                                       aTGGTTAATCCTCATATCACGGGTGGAGTGTTAAGAACATACATAAATGGAGTCATGTTTCCCCTTTTCCATTTATCAAGTTCCTGTTGCCGTTTTAGTc                                                              <  1:1933054/100‑1 (MQ=255)                                         ggTTAATCCTCATATCACGGGTGGAGTGTTAAGAACATACATAAATGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCTGTTGCCGTTTTAGTCCa                                                            >  1:556714/1‑100 (MQ=255)                                         ggTTAATCCTCATATCACGGGTGGAGTGTTAAGAACATACATAAATGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCTGTTGCCGTTTTAGTCCa                                                            >  1:1658933/1‑100 (MQ=255)                                         ggTTAATCCTCATATCACGGGTGGAGTGTTAAGAACATACATAAATGGAGTCATGTTTCCCCTTTTCCATTTATCAAGTTCCTGTTGCCGTTTTAGTCCa                                                            >  1:273199/1‑100 (MQ=255)                                             aaTCCTCATATCACGGGTGGAGTGTTAAGAACATACATAAATGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCTGTTGCCGTTTTAGTCCAtctc                                                        <  2:2145826/100‑1 (MQ=255)                                             aaTCCTCATATCACGGGTGGAGTGTTAAGAACATACATAAATGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCTGTTGCCGTTTTAGTCCAtctc                                                        <  2:1584984/100‑1 (MQ=255)                                              aTCCTCATATCACGGGTGGAGTGTTAAGAACATACATAAATGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCTGTTGCCGTTTTAGTCCATctct                                                       >  1:2135893/1‑100 (MQ=255)                                                 cTCATATCACGGGTGGAGTGTTAAGAACATACATAAATGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCTGTTGCCGTTTTAGTCCATCTCTaa                                                     >  1:690714/1‑99 (MQ=255)                                                  tCATATCACGGGTGGAGTGTTAAGAACATACATAAATGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCTGTTGCCGTTTTAGTCCATCTCTAAtt                                                   <  1:1633001/100‑1 (MQ=255)                                                   cATATCACGGGTGTAGTGTTAAGAACATACATAAATGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCTGTTGCCGTTTTAGTCCATCTCTAATTg                                                  <  2:610718/100‑1 (MQ=255)                                                   cATATCACGGGTGGAGTGTTAAGAACATACATAAATGGAGTCATGTTTCCCCTTTTCCATTTATCAAGTTCCTGTTGCCGTTTTAGTCCATCTCTAATTg                                                  >  2:1204522/1‑100 (MQ=255)                                                     tatCACTGGTGTATTGTTAAGATCATACATATATGGAGTCATGTTTTCTTTTTTCCATTTATCAAGTTCCTGTTGCCGTTTTAGTCCATCTCTAATTGCa                                                <  2:1229718/100‑1 (MQ=255)                                                     tatCACGGGTGGAGTGTTAAGAACATACATAAATGGAGTCATGTTTTCCCTTTTCCGTTTATCAAGTTCCTGTTGCCGTTTTAGTCCATCTCTAATTGCa                                                >  2:1429105/1‑100 (MQ=255)                                                     tatCACGGGTGGAGTGTTAAGAACATACATAAATGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCTGTTGCCGTTTTAGTCCATCTCTAATTGCa                                                >  2:28063/1‑100 (MQ=255)                                                     tatCACGGGTGGAGTGTTAAGAACATACATAAATGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCTGTTGCCGTTTTAGTCCATCTCTAATTGCa                                                >  1:191759/1‑100 (MQ=255)                                                     tatCACGGGTGGAGTGTTAAGAACATACATAAATGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCTGTTGCCGTTTTAGTCCATCTCTAATTGCa                                                >  2:1452786/1‑100 (MQ=255)                                                     tatCACGGGTGGAGTGTTAAGAACATACATAAATGGAGTCATGTTTCCCCTTTTCCATTTATCAAGTTCCTGTTGCCGTTTTAGTCCATCTCTAATTGCa                                                <  1:1963719/100‑1 (MQ=255)                                                       tCACGGGTGGAGTGTTAAGAACATACATAAATGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCTGTTGCCGTTTTAGTCCATCTCTAATTGCata                                              >  1:336850/1‑100 (MQ=255)                                                        cACGGGTTGAGTGTTAAGAACATACATAAATGGAGTCATATTTTCCCTTTTCCATTTATCATGTTCCTGTTGCCGTTTTAGTCCATCTCTAATTTCata                                              <  2:1148233/99‑1 (MQ=255)                                                         aCGGGTGGAGTGTTAAGAACATACATAAATGGAGTCATGTTTCCCCTTTTCCATTTATCAAGTTCCTGTTGCCGTTTTAGTCCATCTCTAATTGCatat                                             <  2:1448051/99‑1 (MQ=255)                                                            ggTGGAGTGTTAAGAACATACATAAATGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCTGTTGCCGTTTTAGTCCATCTCTAATTGCATATTTTa                                         >  1:1551825/1‑100 (MQ=255)                                                            ggTGGAGTGTTAAGAACATACATAAATGGAGTCATGTTTCCCCTTTTCCATTTATCAAGTTCCTGTTGCCGTTTTAGTCCATCTCTAATTGCATATTTTa                                         >  2:738579/1‑100 (MQ=255)                                                                    gtTAAGAACATACATAAATGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCTGTTGCCGTTTTAGTCCATCTCTAATTGCATATTTTAATTTTTCt                                 >  1:66340/1‑100 (MQ=255)                                                                    gtTAAGAACATACATAAATGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCTGTTGCCGTTTTAGTCCATCTCTAATTGCATATTTTAATTTTTCt                                 <  1:1934445/100‑1 (MQ=255)                                                                    gtTAAGAACATACATAAATGGAGTCATGTTTCCCCTTTTCCATTTATCAAGTTCCTGTTGCCGTTTTAGTCCATCTCTAATTGCATATTTTAATTTTTCt                                 >  2:635523/1‑100 (MQ=255)                                                                     ttAAGAACATACATAAATGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCTGTTGCCGTTTTAGTCCATCTCTAATTGCATATTTTAATTTTTCTg                                >  1:1868821/1‑100 (MQ=255)                                                                       aaGAACATACATAAATGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCTGTTGCCGTTTTAGTCCATCTCTAATTGCATATTTTAATTTTTCTGAt                              <  1:1104533/100‑1 (MQ=255)                                                                        aGAACATACATAAATGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCTGTTGCCGTTTTAGTCCATCTCTAATTGCATATTTTAATTTTTCTGATa                             >  1:1867517/1‑100 (MQ=255)                                                                        aGAACATACATAAATGGAGTCATGTTTCCCCTTTTCCATTTATCAAGTTCCTGTTGCCGTTTTAGTCCATCTCTAATTGCATATTTTAATTTTTCTGATa                             <  2:225113/100‑1 (MQ=255)                                                                        aGAACATACATAAATGGAGTCATGTTTCCCCTTTTCCATTTATCAAGTTCCTGTTGCCGTTTTAGTCCATCTCTAATTGCATATTTTAATTTTTCTGATa                             >  2:1753292/1‑100 (MQ=255)                                                                        aGAACATACATAAATGGAGTCATGTTTCCCCTTTTCCATTTATCAAGTTCCTGTTGCCGTTTTAGTCCATCTCTAATTGCATATTTTAATTTTTCTGATa                             >  1:2196107/1‑100 (MQ=255)                                                                            catacataAATTGAGTCATGTTTTCTTTTTTCCATTTATCTTTTTTCTGTTTTTTTTTTATTCTATTTCTATTTGTATATTTTAATTTTTCTGAtn                             <  2:477798/96‑2 (MQ=255)                                                                            catacataAATGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCTGTTGCCGTTTTAGTCCATCTCTAATTGCATATTTTAATTTTTCTGATAAATg                         <  2:749026/100‑1 (MQ=255)                                                                            catacataAATGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCTGTTGCCGTTTTAGTCCATCTCTAATTGCATATTTTAATTTTTCTGATAAATg                         >  2:1151564/1‑100 (MQ=255)                                                                             atacataAATGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCTGTTGCCGTTTTAGTCCATCTCTAATTGCATATTTTAATTTTTCTGATAAATgg                        >  2:281662/1‑100 (MQ=255)                                                                              atcttaAATGGAGTCATGTTTCCCCTTTTCCATTTATCAAGTTCCTGTTGCCGTTTTAGTCCATCTCTAATTGCATATTTTAATTTTTCTGATAAATGGc                       <  2:1646243/96‑1 (MQ=255)                                                                                cataAATGGAGTCATGTTTCCCCTTTTCCATTTATCAAGTTCCTGTTGCCGTTTTAGTCCATCTCTAATTGCATATTTTAATTTTTCTGATAAATGGAAt                     >  2:822796/1‑100 (MQ=255)                                                                                  taAATGGAGTCATGTTTTCCCTTTTCAATTTATCAAGTTCCTGTTGCCGTTTTAGTCCATCTCTAATTGAATATTTTAATTTTTCTGATAAATGGCAtt                    >  2:1699104/1‑99 (MQ=255)                                                                                  taAATGGAGTCATGTTTCCCCTTTTCCATTTATCAAGTTCCTGTTGCCGTTTTAGTCCATCTCTAATTGCATATTTTAATTTTTCTGATAACTGGCATTg                   <  2:1795258/100‑1 (MQ=255)                                                                                  taAATGGAGTCATGTTTCCCCTTTTCCATTTATCAAGTTCCTGTTGCCGTTTTAGTCCATCTCTAATTGCATATTTTAATTTTTCTGATAAATGGCagat                   >  1:1646254/1‑97 (MQ=255)                                                                                   aaaTGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCTGTTGCCGTTTTAGTCCATCTCTAATTGCATATTTTAATTTTTCTGATAAATGGCATTGa                  <  1:1968920/100‑1 (MQ=255)                                                                                    aaTGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCTGTTGCCGTTTTAGTCCATCTCTAATTGCATATTTTAATTTTTCTGATAAATGGCATTGAg                 >  1:1723818/1‑100 (MQ=255)                                                                                    aaTGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCTGTTGCCGTTTTAGTCCATCTCTAATTGCATATTTTAATTTTTCTGATAAATGGCATTGAg                 >  1:1815797/1‑100 (MQ=255)                                                                                    aaTGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCTGTTGCCGTTTTAGTCCATCTCTAATTGCATATTTTAATTTTTCTGATAAATGGCATTGAg                 <  1:1497934/100‑1 (MQ=255)                                                                                    aaTGGAGTCATGTTTTCCCTTTTCCATGTATCAAGTTCCTGTTGCCGTTTTAGTCCATCTCTAATTGCATATTTTAATTTTTCTGATAAATGGCATTGAg                 <  1:184413/100‑1 (MQ=255)                                                                                      tGGAGTCATGTTTCCCCTTTTCCATTTATCAAGTTCCTGTTGCCGTTTTAGTCCATCTCTAATTGCATATTTTAATTTTTCTGATAAATGGCATTGAGCa               <  2:995051/100‑1 (MQ=255)                                                                                         aGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCTGTTGCCGTTTTAGTCCATCTCTAATTGCATATTTTAATTTTTCTGATAAATGGCATTGAGCATCg            <  2:390163/100‑1 (MQ=255)                                                                                         aGTCATGTTTCCCCTTTTCCATTTATCAAGTTCCTGTTGCCGTTTTAGTCCATCTCTAATTGCATATTTTAATTTTTCTGATAAATGGCATTGAGCATCg            <  2:318168/100‑1 (MQ=255)                                                                                          gTCATGTTTTCCCTTTTCCATTTATCAAGTTCCTGTTGCCGTTTTAGTCCATCTCTAATTGCATATTTTAATTTTTCTGATAAATGGCATTGAGCATCGa           >  2:592961/1‑100 (MQ=255)                                                                                          gTCATGTTTTCCCTTTTCCATTTATCAAGTTCCTGTTGCCGTTTTAGTCCATCTCTAATTGCATATTTTAATTTTTCTGATAAATGGCATTGAGCATCGa           <  2:1604568/100‑1 (MQ=255)                                                                                          gTCATGTTTTCCCTTTTCCATTTATCAAGTTCCTGTTGCCGTTTTAGTCCATCTCTAATTGCATATTTTAATTTTTCTGATAAATGGCATTGAGCATCGa           <  1:1454428/100‑1 (MQ=255)                                                                                          gTCATGTTTCCCCTTTTCCATTTATCAAGTTCCTGTTGCCGTTTTAGTCCATCTCTAATTGCATATTTTAATTTTTCTGATAAATGGCATTGAGCATCGa           >  1:1616135/1‑100 (MQ=255)                                                                                            cATGTTTTCCCTTTTCCATTTATCAAGTTCCTGTTGCCGTTTTAGTCCATCTCTAATTGCATATTTTAATTTTTCTGATAAATGGCCTTGAGCATCg            >  1:1674407/1‑97 (MQ=255)                                                                                            cATGTTTTCCCTTTTCCATTTATCAAGTTCCTGTTGCCGTTTTAGTCCATCTCTAATTGCATATTTTAATTTTTCTGATAAATGGCATTGAGCATCGAtt         <  2:896506/100‑1 (MQ=255)                                                                                                ttttCCCTTTTCCATTTATCAAGTTTTTGTTGCAGTTTTATTCCATCTCTAATTGCATATTTTAATTTTTCTGATAAATGGCATTGAGCATCGATTTCAt     <  2:545107/100‑1 (MQ=255)                                                                                                  ttCCCTTTTCCATTTATCAAGTTCCTGTTGCCGTTTTAGTCCATCTCTAATTGCATATTTTAATTTTTCTGATAAATGGCATTGAGCATCGATTTCAttt   <  2:2005222/100‑1 (MQ=255)                                                                                                   tCCCTTTTCCATTTATCAAGTTCCTGTTGCCGTTTTAGTCCATCTCTAATTGCATATTTTAATTTTTCTGATAAATGGCATTGAGCATCGATTTCATTTa  >  2:1141913/1‑100 (MQ=255)                                                                                                   |                                                                                                    ACATCGAAAGAAACTCTTTTAATTTAACATGTAAACGCATGGTTAATCCTCATATCACGGGTGGAGTGTTAAGAACATACATAAATGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCTGTTGCCGTTTTAGTCCATCTCTAATTGCATATTTTAATTTTTCTGATAAATGGCATTGAGCATCGATTTCATTTA  >  NC_000913/238697‑238894

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