Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,810,756	G→A	intergenic (+77/‑14)	ygaH → / → mprA	putative L‑valine exporter, norvaline resistance protein/transcriptional repressor of microcin B17 synthesis and multidrug efflux

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,810,756	0	G	A	94.1%	261.7 / 5.2	86	intergenic (+77/‑14)	ygaH/mprA	putative L‑valine exporter, norvaline resistance protein/transcriptional repressor of microcin B17 synthesis and multidrug efflux
Reads supporting (aligned to +/- strand):  ref base G (3/2);  new base A (35/45);  total (38/48)
Fisher's exact test for biased strand distribution p-value = 6.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.63e-01

AAGTGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCT  >  NC_000913/2810660‑2810852                                                                                                 |                                                                                                 aaGTGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAg                                                                                                <  1:483211/99‑1 (MQ=255)  aGTGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGt                                                                                              <  2:1447904/100‑1 (MQ=255)   gTGATGGCTATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTc                                                                                             <  1:1927127/100‑1 (MQ=255)   gTGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTc                                                                                             >  1:260399/1‑100 (MQ=255)   gTGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTc                                                                                             <  1:1762541/100‑1 (MQ=255)   gTGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTc                                                                                             <  2:173482/100‑1 (MQ=255)   gTGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTc                                                                                             <  1:532322/100‑1 (MQ=255)     gATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCAt                                                                                           >  1:1011970/1‑100 (MQ=255)      aTGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATa                                                                                          >  2:1499410/1‑100 (MQ=255)       tGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATAc                                                                                         <  1:1751631/100‑1 (MQ=255)        ggCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATAcc                                                                                        <  2:139761/100‑1 (MQ=255)        ggCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATAcc                                                                                        >  2:69710/1‑100 (MQ=255)         gCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATAccc                                                                                       <  2:1439798/100‑1 (MQ=255)          cGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCa                                                                                      >  1:1035872/1‑100 (MQ=255)              tatataACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATg                                                                                  <  1:493457/100‑1 (MQ=255)               atataACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATgg                                                                                 >  2:1274115/1‑100 (MQ=255)                tataACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGa                                                                                <  1:240160/100‑1 (MQ=255)                 ataACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGAt                                                                               >  1:1484555/1‑100 (MQ=255)                  taACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATa                                                                              >  1:709388/1‑100 (MQ=255)                   aaCCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAg                                                                             >  2:695089/1‑100 (MQ=255)                   aaCCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAg                                                                             <  1:1093961/100‑1 (MQ=255)                    aCCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGt                                                                            <  1:143709/100‑1 (MQ=255)                    aCCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGt                                                                            <  2:1600757/100‑1 (MQ=255)                      ccACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGtt                                                                           >  1:1272518/1‑99 (MQ=255)                       cacaAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCg                                                                         <  2:1444019/100‑1 (MQ=255)                        acaAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCg                                                                         >  2:139858/1‑99 (MQ=255)                         caAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGt                                                                        >  2:1926403/1‑99 (MQ=255)                          aaGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGttt                                                                      <  1:618427/100‑1 (MQ=255)                           aGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTa                                                                     >  1:758750/1‑100 (MQ=255)                           aGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTa                                                                     >  2:7796/1‑100 (MQ=255)                            gAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTAc                                                                    <  2:995832/100‑1 (MQ=255)                             aaTCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACg                                                                   >  1:1285625/1‑100 (MQ=255)                              aTCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTAc                                                                    >  1:595190/1‑98 (MQ=255)                                 aTTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCa                                                               >  2:1087203/1‑100 (MQ=255)                                    tttCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTg                                                            <  2:1475203/100‑1 (MQ=255)                                       cTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGaan                                                         <  1:1516605/100‑2 (MQ=255)                                       cTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAAc                                                         >  1:114975/1‑100 (MQ=255)                                        tAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACa                                                        >  2:82471/1‑100 (MQ=255)                                        tAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGCTCATACGCAAATGGATAGTTCGTTTACGCCCATTGAACa                                                        >  1:1365860/1‑100 (MQ=255)                                        tAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGATCATACACAAATGGATAGTTCGTTTACGCCCATTGAACa                                                        >  1:1366048/1‑100 (MQ=255)                                         aaaaCAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACaa                                                       >  1:627355/1‑100 (MQ=255)                                         aaaaCAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACaa                                                       >  1:880013/1‑100 (MQ=255)                                          aaaCAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACaaa                                                      <  2:401131/100‑1 (MQ=255)                                          aaaCAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCACATTGAACaaa                                                      >  1:566035/1‑100 (MQ=255)                                           aaCAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAAt                                                     <  1:1263780/100‑1 (MQ=255)                                             cAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGc                                                   <  1:1556335/100‑1 (MQ=255)                                                 aCATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATCAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTaaa                                               <  1:408729/100‑1 (MQ=255)                                                   aTTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACATATGCTAAAAt                                             <  1:1368954/100‑1 (MQ=255)                                                   aTTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAAt                                             <  2:665350/100‑1 (MQ=255)                                                       aCTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTcg                                         <  1:1777110/100‑1 (MQ=255)                                                       aCTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTcg                                         <  1:1198185/100‑1 (MQ=255)                                                       aCTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTcg                                         >  2:1466496/1‑100 (MQ=255)                                                       aCTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTcg                                         <  2:1365868/100‑1 (MQ=255)                                                       aCTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTcg                                         <  2:1285632/100‑1 (MQ=255)                                                        cTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTcgc                                        >  2:1566435/1‑100 (MQ=255)                                                         tttATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTcgcg                                       <  2:758751/100‑1 (MQ=255)                                                          ttATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCgcgc                                      >  1:745263/1‑100 (MQ=255)                                                          ttATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCgcgc                                      <  2:485384/100‑1 (MQ=255)                                                          ttATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCgcgc                                      <  1:1659558/100‑1 (MQ=255)                                                             tttGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAg                                   <  1:1473330/100‑1 (MQ=255)                                                              ttGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAgc                                  <  2:231598/100‑1 (MQ=255)                                                              ttGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAgc                                  <  1:1760703/100‑1 (MQ=255)                                                               tGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAgcc                                 >  1:496714/1‑100 (MQ=255)                                                                gTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAgccg                                <  1:1550236/100‑1 (MQ=255)                                                                     tGTCGCTACTATATCGGCTGTACTTAAGCAGGTCATACGCACATGGATAGTTCGCTTACGCCCATTGAACAAATGTTAAAATTTCGCGcc                                     <  1:379568/90‑1 (MQ=255)                                                                        cGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAg                        >  2:1263875/1‑100 (MQ=255)                                                                        cGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAg                        >  1:430724/1‑100 (MQ=255)                                                                         gTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGa                       <  1:1913966/100‑1 (MQ=255)                                                                          ttACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGAt                      >  2:1080280/1‑100 (MQ=255)                                                                           tACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGAtt                     <  2:579575/100‑1 (MQ=255)                                                                           tACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGAtt                     >  1:1379156/1‑100 (MQ=255)                                                                             cTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGAtttt                   >  2:737980/1‑100 (MQ=255)                                                                              tataTCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTc                  <  2:1719832/100‑1 (MQ=255)                                                                               atatCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTcc                 <  2:481216/100‑1 (MQ=255)                                                                               atatCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTcc                 <  2:1289632/100‑1 (MQ=255)                                                                               atatCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTcc                 <  2:752837/100‑1 (MQ=255)                                                                               atatCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTcc                 <  1:82471/100‑1 (MQ=255)                                                                               atatCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTcc                 <  1:573723/100‑1 (MQ=255)                                                                                   cGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAACTTCCGCGCCAGCCGCCACGAAGATTTTCCTTAt             <  1:1281049/100‑1 (MQ=255)                                                                                    ggCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTagc            >  1:1183717/1‑98 (MQ=255)                                                                                      cTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAg          >  2:1549531/1‑100 (MQ=255)                                                                                       tGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAgg         >  2:1174781/1‑100 (MQ=255)                                                                                       tGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAgg         >  2:964725/1‑100 (MQ=255)                                                                                       tGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAgg         <  2:496716/100‑1 (MQ=255)                                                                                         aaaTTAATAATGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGTGCCAGCCGCCACGAAGATTTTCCTTATCAGgag       <  1:1274110/100‑1 (MQ=255)                                                                                           aataataaGGTCAGACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGAt     >  2:1319406/3‑100 (MQ=255)                                                                                             taataaGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATcc   >  1:114202/1‑100 (MQ=255)                                                                                             taataaGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATcc   >  2:784645/1‑100 (MQ=255)                                                                                             taataaGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATc    >  1:1660867/1‑99 (MQ=255)                                                                                              aaTAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCt  >  1:912351/1‑100 (MQ=255)                                                                                                 |                                                                                                 AAGTGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCT  >  NC_000913/2810660‑2810852

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