Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,708,468	A→T	5.8%	E517V (GAA→GTA) ‡	rsxC →	SoxR iron‑sulfur cluster reduction factor component; putative membrane‑associated NADH oxidoreductase of electron transport complex

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,708,468	0	A	T	5.8%	246.4 / 5.2	69	E517V (GAA→GTA) ‡	rsxC	SoxR iron‑sulfur cluster reduction factor component; putative membrane‑associated NADH oxidoreductase of electron transport complex
Reads supporting (aligned to +/- strand):  ref base A (37/28);  new base T (2/2);  total (39/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.48e-01

GCGAAAAAGCGGCTCGCCTTGAACGACATAAGAGCGCAGCCGTTCAACCTGCAGCCAAAGATAAAGATGCGATTGCTGCCGCTCTGGCGCGGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGC  >  NC_000913/1708324‑1708615                                                                                                                                                 |                                                                                                                                                    gcgAAAAAGCGGCTCGCCTTGAACGACATAAGAGCGCAGCCGTTCAACCTGCAGCCAAAGATAAAGATGCGATTGCTGCCGCTCTGGCGCGGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACg                                                                                                                                                  <  2:21054/148‑1 (MQ=255)   gAAAAAGCGGCTCGCCTTGAACGACATAAGAGCGCAGCCGTTCAACCTGCAGCCAAAGATAAAGATGCGATTGCTGCCGCTCTGGCGCGGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGccc                                                                                                                                               <  4:682628/149‑1 (MQ=255)   gAAAAAGCGGCTCGCCTTGAACGACATAAGAGCGCAGCCGTTCAACCTGCAGCCAAAGATAAAGATGCGATTGCTGCCGCTCTGGCGCGGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGcc                                                                                                                                                <  2:499421/148‑1 (MQ=255)    aaaaaGCGGCTCGCCTTGAACGACATAAGAGCGCAGCCGTTCAACCTGCAGCCAAAGATAAAGATGCGATTGCTGCCGCTCTGGCGCGGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGccc                                                                                                                                               >  2:582450/1‑148 (MQ=255)    aaaaaGCGGCTCGCCTTGAACGACATAAGAGCGCAGCCGTTCAACCTGCAGCCAAAGATAAAGATGCGATTGCTGCCGCTCTGGCGCGGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGccc                                                                                                                                               >  1:741036/1‑148 (MQ=255)     aaaaGCGGCTCGCCTTGAACGACATAAGAGCGCAGCCGTTCAACCTGCAGCCAAAGATAAAGATGCGATTGCTGCCGCTCTGGCGCGGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGccc                                                                                                                                               >  4:222331/1‑147 (MQ=255)     aaaaGCGGCTCGCCTTGAACGACATAAGAGCGCAGCCGTTCAACCTGCAGCCAAAGATAAAGATGCGATTGCTGCCGCTCTGGCGCGGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGccc                                                                                                                                               <  3:222331/147‑1 (MQ=255)          cGGCTCGCCTTGAACGACATAAGAGCGCAGCCGTTCAACCTGCAGCCAAAGATAAAGATGCGATTGCTGCCGCTCTGGCGCGGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATaa                                                                                                                                         <  1:491618/148‑1 (MQ=255)            gCTCGCCTTGAACGACATAAGAGCGCAGCCGTTCAACCTGCAGCCAAAGATAAAGATGCGATTGCTGCCGCTCTGGCGCGGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAg                                                                                                                                      >  3:619815/1‑149 (MQ=255)                      aaCGACATAAGAGCGCAGCCGTTCAACCTGCAGCCAAAGATAAAGATGCGATTGCTGCCGCTCTGGCGCGGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAattat                                                                                                                             >  2:441897/1‑148 (MQ=255)                        cGACATAAGAGCGCAGCCGTTCAACCTGCAGCCAAAGATAAAGATGCGATTGCTGCCGCTCTGGCGCGGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGTTCGCCCgcttta                                                                                                                                         >  1:102140/1‑129 (MQ=255)                        cGACATAAGAGCGCAGCCGTTCAACCTGCAGCCAAAGATAAAGATGCGATTGCTGCCGCTCTGGCGCGGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGTTCGCCCgcttta                                                                                                                                         <  2:102140/134‑6 (MQ=255)                        cGACATAAGAGCGCAGCCGTTCAACCTGCAGCCAAAGATAAAGATGCGATTGCTGCCGCTCTGGCGCGGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTgc                                                                                                                          <  3:197285/149‑1 (MQ=255)                        cGACATAAGAGCGCAGCCGTTCAACCTGAAGCCAAAGATAAAGATGCGATTGCTGCCGCTCTGGCGCGGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTgc                                                                                                                          <  3:871656/149‑1 (MQ=255)                         gACATAAGAGCGCAGCCGTTCAACCTGCAGCCAAAGATAAAGATGCGATTGCTGCCGCTCTGGCGCGGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTgca                                                                                                                         <  3:84100/149‑1 (MQ=255)                            aTAAGAGCGCAGCCGTTCAACCTGCAGCCAAAGATAAAGATGCGATTGCTGCCGCTCTGGCGCGGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGTTCGCCCgct                                                                                                                                            <  3:832549/127‑3 (MQ=255)                            aTAAGAGCGCAGCCGTTCAACCTGCAGCCAAAGATAAAGATGCGATTGCTGCCGCTCTGGCGCGGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGTTCGCCCgct                                                                                                                                            >  4:832549/1‑125 (MQ=255)                              aaGAGCGCAGCCGTTCAACCTGCAGCCAAAGATAAAGATGCGATTGCTGCCTCTCTGGCGCGGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTTCAATTATTGCAGCAc                                                                                                                     >  2:263887/1‑148 (MQ=255)                                    gcAGCCGTTCAACCTGCAGCCAAAGATAAAGATGCGATTGCTGCCGCTCTGGCGCGGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGc                                                                                                              >  4:849396/1‑149 (MQ=255)                                      aGCCGTTCAACCTGCAGCCAAAGATAAAGATGCGATTGCTGCCGCTCTGGCGCGGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGccc                                                                                                            >  3:758323/1‑149 (MQ=255)                                          gTTCAACCTGCAGCCAAAGATAAAGATGCGATTGCTGCCGCTCTGGCGCGGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTaa                                                                                                        >  4:710164/1‑149 (MQ=255)                                             cAACCTGCAGCCAAAGATAAAGATGCGATTGCTGCCGCTCTGGCGCGGGTGAAAGAGAAACAGGCTCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGCATAACAGTGCAATTATTGCAGTACGGGAAGCCCGTAAAg                                                                                                      <  1:263887/148‑1 (MQ=255)                                             cAACCTGCAGCCAAAGATAAAGATGCGATTGCTGCCGCTCTGGCGCGGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAg                                                                                                      <  1:411786/148‑1 (MQ=255)                                                 cTGCAGCCAAAGATAAAGATGCGATTGCTGCCGCTCTGGCGCGGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCaa                                                                                                 <  4:394328/149‑1 (MQ=255)                                                  tGCAGCCAAAGATAAAGATGCGATTGCTGCCGCTCTGGCGCGGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAat                                                                                                                                >  4:415865/1‑117 (MQ=255)                                                  tGCAGCCAAAGATAAAGATGCGATTGCTGCCGCTCTGGCGCGGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAat                                                                                                                                <  3:415865/117‑1 (MQ=255)                                                       ccAAAGATAAAGATGCGATTGCTGCCGCTCTGGCGCGGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCaga                                                                                           >  4:427108/1‑149 (MQ=255)                                                        cAAAGATAAAGATGCGATTGCTGCCGCTCTGGCGCGGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCagag                                                                                          <  4:323419/149‑1 (MQ=255)                                                          aagataaagatGCGATTGCTGCCGCTCTGGCGCGGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGc                                                                                         >  2:67008/1‑148 (MQ=255)                                                              taaagatGCGATTGCTGCCGCTCTGGCGCGGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGaaa                                                                                     >  1:907802/1‑148 (MQ=255)                                                                aagatGCGATTGCTGCCGCTCTGGCGCGGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAg                                                                                  <  4:215851/149‑1 (MQ=255)                                                                    tGCGATTGCTGCCGCTCTGGCGCGGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCa                                                                               >  2:129954/1‑148 (MQ=255)                                                                      cGATTGCTGCCGCTCTGGCGCGGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGa                                                                             <  2:614175/148‑1 (MQ=255)                                                                      cGATTGCTGCCGCTCTGGCGCGGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGa                                                                             <  1:441897/148‑1 (MQ=255)                                                                           gctgcCGCTCTGGCGCGGGTGGAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTgca                                                                       <  4:171564/149‑1 (MQ=255)                                                                           gctgcCGCTCTGGCGCGGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTgc                                                                        <  1:129954/148‑1 (MQ=255)                                                                              gcCGCTCTGGCGCGGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTgcagca                                                                    >  4:203469/1‑149 (MQ=255)                                                                                  ctctGGCGCGGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACt                                                                >  3:54997/1‑149 (MQ=255)                                                                                   tctGGCGCGGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACt                                                                >  2:39804/1‑148 (MQ=255)                                                                                        cgcgGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAacga                                                           <  2:907802/148‑1 (MQ=255)                                                                                         gcgGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAacgac                                                          <  1:2362/148‑1 (MQ=255)                                                                                           gggTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACgc                                                        <  1:39804/148‑1 (MQ=255)                                                                                             gTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACgcagc                                                     >  4:610151/1‑149 (MQ=255)                                                                                                 aaGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAAc                                                  >  2:219237/1‑148 (MQ=255)                                                                                                   gagaAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCg                                                >  2:192323/1‑148 (MQ=255)                                                                                                     gaAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGtt                                              >  1:201621/1‑148 (MQ=255)                                                                                                     gaAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGtt                                              >  2:317293/1‑148 (MQ=255)                                                                                                      aaaCAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTg                                             <  1:582450/148‑1 (MQ=255)                                                                                                       aaCAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCt                                           >  3:428698/1‑149 (MQ=255)                                                                                                         cAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGa                                         <  4:54997/149‑1 (MQ=255)                                                                                                          aGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGAt                                        <  3:849396/149‑1 (MQ=255)                                                                                                                aGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTa                                                               >  2:809931/1‑120 (MQ=255)                                                                                                                aGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTa                                                               <  1:809931/120‑1 (MQ=255)                                                                                                                aGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGt                                  >  3:442382/1‑149 (MQ=255)                                                                                                                     acacAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAAc                             >  3:875617/1‑149 (MQ=255)                                                                                                                            cTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGt                       <  1:67008/148‑1 (MQ=255)                                                                                                                               ttGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGaa                   >  4:530218/1‑149 (MQ=255)                                                                                                                                tgtgATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGaa                   >  1:803587/1‑148 (MQ=255)                                                                                                                                 gtgATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAag                  >  1:653908/1‑148 (MQ=255)                                                                                                                                  tgATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAagca                >  3:731316/1‑149 (MQ=255)                                                                                                                                   gATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAagcag               <  3:710164/149‑1 (MQ=255)                                                                                                                                     ttAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTAGAagcagc              >  2:786391/1‑148 (MQ=255)                                                                                                                                     ttAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGACAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCt             >  4:536145/1‑149 (MQ=255)                                                                                                                                      tAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCt             >  2:546894/1‑148 (MQ=255)                                                                                                                                      tAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCt             >  2:409655/1‑148 (MQ=255)                                                                                                                                          gCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGc        >  4:902633/1‑149 (MQ=255)                                                                                                                                             ggCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTaaaa                              >  4:221602/1‑124 (MQ=255)                                                                                                                                             ggCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTaaaa                              <  3:221602/124‑1 (MQ=255)                                                                                                                                                gAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCgcgc  >  4:577830/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                    GCGAAAAAGCGGCTCGCCTTGAACGACATAAGAGCGCAGCCGTTCAACCTGCAGCCAAAGATAAAGATGCGATTGCTGCCGCTCTGGCGCGGGTGAAAGAGAAACAGGCCCAGGCTACACAGCCTATTGTGATTAAAGCGGGCGAACGCCCGGATAACAGTGCAATTATTGCAGCACGGGAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGC  >  NC_000913/1708324‑1708615

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑