Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,121,597	(T)5→4	coding (151/1332 nt)	hslU ←	molecular chaperone and ATPase component of HslUV protease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,121,593	0	T	.	97.3%	322.8 / 2.6	73	coding (155/1332 nt)	hslU	molecular chaperone and ATPase component of HslUV protease
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base . (41/30);  total (42/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.21e-01

TTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTC  >  NC_000913/4121494‑4121687                                                                                                    |                                                                                               ttGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGatat                                                                                                <  1:1844366/100‑1 (MQ=255) ttGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGatat                                                                                                <  1:1843694/100‑1 (MQ=255) ttGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGatat                                                                                                <  2:1169183/100‑1 (MQ=255) ttGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGatat                                                                                                <  1:726724/100‑1 (MQ=255)  tGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATAtt                                                                                               >  2:1164191/1‑100 (MQ=255)   ggTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATAttt                                                                                              >  1:210637/1‑100 (MQ=255)    gTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATAtttt                                                                                             >  1:732113/1‑100 (MQ=255)    gTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATAtttt                                                                                             >  2:680390/1‑100 (MQ=255)      cGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATAttttcg                                                                                           <  2:1187188/100‑3 (MQ=255)        cTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATAttttcggg                                                                                         <  2:76963/100‑5 (MQ=255)        cTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATAttttcggg                                                                                         <  2:1739551/100‑5 (MQ=255)        cTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATAttttcggg                                                                                         <  1:1852076/100‑5 (MQ=255)         ttCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCgggg                                                                                       <  2:277733/100‑1 (MQ=255)         ttCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCgggg                                                                                       <  2:723004/100‑1 (MQ=255)         ttCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCgggg                                                                                       <  1:1908723/100‑1 (MQ=255)         ttCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCgggg                                                                                       >  2:1553662/1‑100 (MQ=255)           cAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTc                                                                                     >  1:579043/1‑100 (MQ=255)            aaCTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCa                                                                                    <  2:548429/100‑1 (MQ=255)             aCTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCAc                                                                                   <  1:1221384/100‑1 (MQ=255)             aCTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCAc                                                                                   <  2:154168/100‑1 (MQ=255)              cTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGACGATCCTACCGATA‑TATTCGGGGTCAc                                                                                   >  2:740824/1‑99 (MQ=255)                 tgatgaACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCa                                                                               >  2:206638/1‑100 (MQ=255)                   atgaACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATg                                                                             >  2:293842/1‑100 (MQ=255)                    tgaACTGCGCATTCGCCAGCTTAGCCAGACGACGTGCGATTTCAGTTGTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATgg                                                                            <  1:903004/100‑1 (MQ=255)                    tgaACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATgg                                                                            >  1:350584/1‑100 (MQ=255)                    tgaACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATgg                                                                            <  2:1895263/100‑1 (MQ=255)                    tgaACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATgg                                                                            >  1:1706255/1‑100 (MQ=255)                     gaACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGgc                                                                           <  1:1576364/100‑1 (MQ=255)                      agcGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTGCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGgc                                                                           <  1:545840/97‑1 (MQ=255)                        cGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCa                                                                        <  2:66896/100‑1 (MQ=255)                        cGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCa                                                                        <  2:38853/100‑1 (MQ=255)                         ggCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAg                                                                       >  2:1224474/1‑100 (MQ=255)                             cATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGctct                                                                   >  1:1732889/1‑100 (MQ=255)                              nttCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCtt                                                                  >  2:2038791/2‑100 (MQ=255)                                tCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTc                                                                 >  2:553129/1‑99 (MQ=255)                                 cGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGt                                                               >  1:1980850/1‑100 (MQ=255)                                 cGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGt                                                               >  2:741355/1‑100 (MQ=255)                                 cGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGt                                                               >  1:800737/1‑100 (MQ=255)                                 cGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGt                                                               >  1:1980900/1‑100 (MQ=255)                                  gCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGtt                                                              >  1:699232/1‑100 (MQ=255)                                   ccAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTg                                                             <  2:1330017/100‑1 (MQ=255)                                     aGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAg                                                           >  1:1279954/1‑100 (MQ=255)                                      gCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGc                                                          >  2:2220760/1‑100 (MQ=255)                                       cTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCt                                                         >  2:1090124/1‑100 (MQ=255)                                        ttAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTg                                                        <  2:1728731/100‑1 (MQ=255)                                          aGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGc                                                       >  1:135458/1‑100 (MQ=255)                                            ccAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATg                                                    <  2:960224/100‑1 (MQ=255)                                            ccAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATg                                                    <  1:1257109/100‑1 (MQ=255)                                             cAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGc                                                   <  1:647538/100‑1 (MQ=255)                                             cAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGc                                                   >  1:1912939/1‑100 (MQ=255)                                              aGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCg                                                  >  1:440533/1‑100 (MQ=255)                                               gacgacGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGa                                                 >  1:828239/1‑100 (MQ=255)                                                acgacgGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGAc                                                >  1:200723/1‑100 (MQ=255)                                                 cgacgGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGAc                                                >  1:432493/1‑99 (MQ=255)                                                 cgacgGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑‑TTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGc                                              >  1:1212105/1‑100 (MQ=255)                                                   acgGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGcc                                             <  1:2220226/100‑1 (MQ=255)                                                      ggCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGc                                          >  1:304107/1‑100 (MQ=255)                                                           tttCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTa                                     >  1:368974/1‑100 (MQ=255)                                                           tttCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTa                                     <  2:1160104/100‑1 (MQ=255)                                                            ttCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTAc                                    >  1:1398940/1‑100 (MQ=255)                                                               aGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCa                                 <  1:755924/100‑1 (MQ=255)                                                               aGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTCTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCa                                 >  1:574156/1‑100 (MQ=255)                                                                gTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCag                                >  1:7864/1‑100 (MQ=255)                                                                    tACCGACACCGGTCGGGCCGATCATCAGGATA‑TGTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCGGCATGCGACGCCAGCGGTTACGCAGAGCa                            >  1:1334811/1‑100 (MQ=255)                                                                      ccGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAAt                          >  2:420352/1‑100 (MQ=255)                                                                        gACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCg                        <  2:1033819/100‑1 (MQ=255)                                                                        gACACCGGTCGGGCCGATCATCAGGATA‑TCTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCg                        <  1:1339303/100‑1 (MQ=255)                                                                         acacCGGTCGGGCCGATCATCAGGATA‑TTGTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGc                       <  1:1340458/100‑1 (MQ=255)                                                                          cacCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGcc                      >  2:1023693/1‑100 (MQ=255)                                                                            ccGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCcac                    >  2:1266315/1‑100 (MQ=255)                                                                            ccGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCcac                    <  2:200723/100‑1 (MQ=255)                                                                             cGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCcaca                   >  2:388340/1‑100 (MQ=255)                                                                             cGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCcaca                   <  2:1279956/100‑1 (MQ=255)                                                                             cGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCcaca                   <  2:1114057/100‑1 (MQ=255)                                                                                  gggCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACg              >  2:1999065/1‑100 (MQ=255)                                                                                   ggCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGc             >  2:1154529/1‑100 (MQ=255)                                                                                   ggCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGc             >  2:689348/1‑100 (MQ=255)                                                                                       gATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTgg         >  1:2157357/1‑100 (MQ=255)                                                                                       gATCATCAGGATA‑‑TTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGc        <  2:1212107/100‑1 (MQ=255)                                                                                        atcatcAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGc        <  2:1990778/100‑1 (MQ=255)                                                                                         tcatcaGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCg       <  2:936292/100‑1 (MQ=255)                                                                                          catcaGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGt      <  1:1136714/100‑1 (MQ=255)                                                                                            tcaGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTg    >  2:1775000/1‑100 (MQ=255)                                                                                              aGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGt   <  2:984569/100‑1 (MQ=255)                                                                                              aGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTc  >  2:6468/1‑100 (MQ=255)                                                                                                    |                                                                                               TTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTC  >  NC_000913/4121494‑4121687

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