Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,703,693	C→T	58.1%	A206T (GCG→ACG)	ydgJ ←	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,703,693	0	C	T	58.1%	32.9 / 96.1	86	A206T (GCG→ACG)	ydgJ	putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (22/14);  new base T (36/14);  total (58/28)
Fisher's exact test for biased strand distribution p-value = 3.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.13e-01

ACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCAAGAAGATGTGGTGCTAAATCGTACCAGA  >  NC_000913/1703597‑1703791                                                                                                 |                                                                                                   aCCGTGCTGACTCATCGGCTGCCAGAATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGTTcc                                                                                                 <  2:2812102/100‑1 (MQ=255)   cGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCgg                                                                                               >  2:1582208/1‑100 (MQ=255)     tGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGc                                                                                             <  1:3578085/100‑1 (MQ=255)     tGCTGACTCAGCGGCTGCCAGCATGGTACCGCGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGc                                                                                             <  1:2859732/100‑1 (MQ=255)        tGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGTTCCGGGCCGt                                                                                          >  1:1297769/1‑100 (MQ=255)         gACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGTTCCGGGCCGTa                                                                                         >  2:180561/1‑100 (MQ=255)         gACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTa                                                                                         >  2:3412918/1‑100 (MQ=255)         gACTCAGAGGCAGCCAGCATGGTACCGTGTAAAATGACTCGCCGCACCGGATAGGACAAGATGGCATGGAAATAATCGGTCGACTGCGCCCCGGGCCGTa                                                                                         >  2:2166096/1‑100 (MQ=255)            tCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGTTCCGGGCCGTAACt                                                                                      >  1:2615745/1‑100 (MQ=255)              aGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGTTCCGGGCCGTAACtgt                                                                                    >  2:742342/1‑100 (MQ=255)              aGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACtgt                                                                                    >  1:2307351/1‑100 (MQ=255)               ttccgatctaGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATTATCGGTCGACTGCGCTCCGGGCCGTAACtgtg                                                                                   <  1:3926645/91‑1 (MQ=255)               gCGGCTGCCCGCATGGGACCGTGTAAAGTGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGTTCCGGGCCGTAACtgtg                                                                                   <  2:1971889/100‑1 (MQ=255)               gCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGTTCCGGGCCGTAACtgtg                                                                                   >  1:1590708/1‑100 (MQ=255)                cGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCAGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGc                                                                                  >  2:1734221/1‑100 (MQ=255)                 ggCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGTTCCGGGCCGTAACtg                                                                                     >  2:326898/1‑96 (MQ=255)                 ggCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGcc                                                                                 >  2:446749/1‑100 (MQ=255)                 ggCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGcc                                                                                 >  1:1643547/1‑100 (MQ=255)                     gCCAGCATGGTACCGCGTAAAATGACTCGCCGCTGTGGATAGGACAAGATAGCGTGGAAATAAGCGGTCGACTGCGCGCCGGGCCGCAACTGTGCCaaa                                                                              >  2:725340/1‑99 (MQ=255)                      ccAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGTTCCGGGCCGTAACTGTGCCAAATc                                                                            >  1:2501341/1‑100 (MQ=255)                      ccAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGTTCCGGGCCGTAACTGTGCCAAATc                                                                            >  1:618011/1‑100 (MQ=255)                      ccAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATc                                                                            >  1:3476095/1‑100 (MQ=255)                       cAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCt                                                                           <  1:1239540/100‑1 (MQ=255)                       cAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCt                                                                           >  2:1934529/1‑100 (MQ=255)                        aGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACtgtgagatcggaa                                                                          >  2:3926619/1‑91 (MQ=255)                             ggTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGTTCCGGGACGTAACTGTTCCAAATCTACCGTc                                                                     >  1:2614842/1‑100 (MQ=255)                              gTACCGTGTAAAATGACTCGCCGGTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTAcc                                                                        <  2:3148493/96‑1 (MQ=255)                              gTACCGTGTAAAATGACTCGCCGCTGTGCATAGGACAAGATGGCGTGGAAACAATCGGTCGACTGCGCTACGGGCCGTAACTGTGCCAAATCTACCGTc                                                                     >  2:2683572/1‑99 (MQ=255)                               tACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGTTCCGGGCCGTAACTGTGCCAAATCTACCGTCAt                                                                   >  2:2028222/1‑100 (MQ=255)                               tACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGTTCCGGGCCGTAACTGTGCCAAATCTACCGTCAt                                                                   >  1:3941484/1‑100 (MQ=255)                               tACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGTTCCGGGCCGTAACTGTGCCAAATCTACCGTCAt                                                                   >  1:3942333/1‑100 (MQ=255)                                 ccGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGTTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGc                                                                 <  2:43870/100‑1 (MQ=255)                                 ccGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGc                                                                 <  2:258793/100‑1 (MQ=255)                                  cGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGTTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCt                                                                >  1:4100812/1‑100 (MQ=255)                                     gtAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGTTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGAc                                                             >  1:2563572/1‑100 (MQ=255)                                      tAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGAcc                                                            >  2:3061791/1‑100 (MQ=255)                                        aaaTGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCgg                                                          >  2:3016219/1‑100 (MQ=255)                                            gACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGTTCCGGGCCGTAACTGTGCCAAACCTACCGTCATGCTGACCGGTaaa                                                      >  1:363473/1‑100 (MQ=255)                                              cTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGTTCGGGGACGTAACTGTGCACAATCTACCGTTATGCTGACGGGTaa                                                       >  2:3734491/1‑97 (MQ=255)                                              cTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGTTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAAcc                                                    >  1:1434578/1‑100 (MQ=255)                                               tCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGTTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTa                                                        <  2:3085026/95‑1 (MQ=255)                                                cgccgcTGTGGATAGGACGAGATGGCGTGGAAATAATCGGTCGACTGCGTTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTaa                                                       <  2:1011682/95‑1 (MQ=255)                                                cgccgcTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGTTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCaa                                                  >  1:1165238/1‑100 (MQ=255)                                                 gccgcTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGTTCCGGGCCGTAACTGTGCCAAATCTTCCGTCATGCTGACCGGTAAACCaaa                                                 >  2:2157218/1‑100 (MQ=255)                                                   cgcTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGTTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCaaaaa                                               >  1:1269337/1‑100 (MQ=255)                                                    gcTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAg                                              >  1:1066405/1‑100 (MQ=255)                                                    gcTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAg                                              >  1:3396049/1‑100 (MQ=255)                                                    gcTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGATAAACCaaaa                                                <  2:417378/98‑1 (MQ=255)                                                     cTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGTTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAg                                              <  2:1524937/99‑1 (MQ=255)                                                     cTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGc                                             <  1:376928/100‑1 (MQ=255)                                                     cTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGc                                             <  1:376909/100‑1 (MQ=255)                                                       gtgGATTGGACAAGATGGCGTGGATATAATCGGTCGACTGCGTGCCGGGCCGTAACTGTGCCAAATCTACCGGCAGGCTGACCGTTAAACCAAAAAGCGt                                           <  2:4119149/100‑1 (MQ=255)                                                       gtgGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGt                                           >  1:3064163/1‑100 (MQ=255)                                                         ggATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGTTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTaa                                         >  1:519433/1‑100 (MQ=255)                                                          gATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGTTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGc                                             <  2:3819876/95‑1 (MQ=255)                                                           agaGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGTTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATg                                       >  1:2795751/3‑100 (MQ=255)                                                           aTAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGTTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATg                                       >  1:1238371/1‑100 (MQ=255)                                                           aTAGAACAAGATAGCGCGGAAATAAACGCTCGACTGCGTTCCGGCACGTAACTGTGCCAAATCTACCGTCATGCTGAAAAAAAAAAAAAAAAGCGTAAtn                                       <  2:2600175/100‑2 (MQ=255)                                                             aGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGTTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGc                                     >  1:2680877/1‑100 (MQ=255)                                                              ggACAAGATGGCGTGGAAATAATCGGTCGACTGCGTTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGcc                                    <  1:943219/100‑1 (MQ=255)                                                              ggACAAGATGGCGTGGAAATAATCGGTCGACTGCGTTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGcc                                    >  2:1511256/1‑100 (MQ=255)                                                              ggACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGcc                                    <  1:1289880/100‑1 (MQ=255)                                                              ggACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGcc                                    >  1:1820180/1‑100 (MQ=255)                                                              ggACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGcc                                    >  1:2313369/1‑100 (MQ=255)                                                               gACAAGATGGCGTGGAAATAATCGGTCGACTGCGTTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGcc                                    >  2:715806/1‑99 (MQ=255)                                                                  aaGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGAt                                <  2:2119086/100‑1 (MQ=255)                                                                   aGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATc                               >  2:3321717/1‑100 (MQ=255)                                                                      tGGGGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATGTCCCGTCATGCTGAAAAAAAAAAAAAAAAGCGTAATGGCCTg                                  <  2:3465812/94‑1 (MQ=255)                                                                       ggCGTGGAAATAATCGGTCGACTGCGTTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCaaga                           <  2:2345246/100‑1 (MQ=255)                                                                       ggCGTGGAAATAATCGGTCGACTGCGTTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCaaga                           >  1:1309184/1‑100 (MQ=255)                                                                         cGTGGAAATAATCGGTCGACTGCGTTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCaaaaag                         >  1:2946526/1‑100 (MQ=255)                                                                         cGTGAAAATAATCGGTCGACTGGGTTGCGGGCCGTACCTGTGCCAAATCTACCGTCATGCTTACCGGTAAACCAAAAAGCGTAATGGCCTGATCaagaa                          <  2:3189488/99‑1 (MQ=255)                                                                          gTGGAAATAATCGGTCGACTGCGTTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCAagaaga                        >  1:2558974/1‑100 (MQ=255)                                                                           tGGAAATAATCGGTCGACTGCGTTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCAAGAAGAt                       >  1:3689215/1‑100 (MQ=255)                                                                           tGGAAATAATCGGTCGACTGCGTTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCAAGAAGAt                       >  2:405520/1‑100 (MQ=255)                                                                             gAAATAATCGGTCGACTGCGTTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCAAGAAGAtg                      <  2:390697/99‑1 (MQ=255)                                                                                ataatCGGTCGACTGCGTTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCAAGAAGATgtggt                  >  2:1223943/1‑100 (MQ=255)                                                                                  natCGGTCGACTGCGTTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCAAGAAGATGTGGTGc                >  2:1710084/2‑100 (MQ=255)                                                                                  aatCGGTCGACTGCGTTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCAAGAAGATGTGGTGc                >  2:3739968/1‑100 (MQ=255)                                                                                    tCGGTCGACTGCGTTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCAAGAAGATGTGGTGCTa              <  1:3592437/100‑1 (MQ=255)                                                                                      ggTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCAAGAAGATGTGGTGCTaaa            <  1:581495/100‑1 (MQ=255)                                                                                       gTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCAAGAAGATGTGGTGCTAAAt           <  2:2697563/100‑1 (MQ=255)                                                                                         cGACTGCGTTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCAAGAAGATGTGGTGCTAAATCg         >  1:981038/1‑100 (MQ=255)                                                                                         cGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCAAGAAGATGTGGTGCTAAATCg         >  1:1251443/1‑100 (MQ=255)                                                                                         cGACTGCGCTCCGGGCCGTAACTGTGCCAAATCAACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCAAGAAGATGTGATGCTAAATCg         >  2:2710265/1‑100 (MQ=255)                                                                                            cTGCGTTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCAAGAAGATGTGGTGCTAAATCGTAc      >  1:1865337/1‑100 (MQ=255)                                                                                              gcgtTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCAAGAAGATGTGGTGCTAAATCGTACCa    <  2:2614826/100‑1 (MQ=255)                                                                                              gcgtTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCAAGAAGATGTGGTGCTAAATCGTACCa    >  2:2219243/1‑100 (MQ=255)                                                                                              gcgcTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCAAGAAGATGTGGTGCTAAATCGTACCa    >  1:3512486/1‑100 (MQ=255)                                                                                                gcTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCAAGAAGATGTGGTGCTAAATCGTACCAGa  >  2:3939700/1‑100 (MQ=255)                                                                                                 |                                                                                                   ACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCAAGAAGATGTGGTGCTAAATCGTACCAGA  >  NC_000913/1703597‑1703791

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