Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,707,029	A→G	100%	V37V (GTA→GTG)	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,707,029	0	A	G	98.9%	300.9 / ‑4.8	91	V37V (GTA→GTG)	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 (0/1);  new base G (41/49);  total (41/50)
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.41e-01

CTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCC  >  NC_000913/1706933‑1707123                                                                                                 |                                                                                               ctctGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGccc                                                                                             >  2:80386/1‑100 (MQ=255) ctctGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGccc                                                                                             >  2:80320/1‑100 (MQ=255)      cATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGgc                                                                                        <  1:47425/100‑1 (MQ=255)       aTTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGgcg                                                                                       <  2:576025/100‑1 (MQ=255)       aTTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCGGCGCCAGGTGCCCCTgg                                                                                         >  2:285616/1‑98 (MQ=255)        ttCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGgcgc                                                                                      <  1:1115385/100‑1 (MQ=255)          cAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAg                                                                                    <  1:196231/100‑1 (MQ=255)          cAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGGGCACCTGGCGCAg                                                                                    <  1:1053160/100‑1 (MQ=255)           aGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGc                                                                                   >  1:272000/1‑100 (MQ=255)           aGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGGGCACCTGGCGCAgg                                                                                   <  1:524484/100‑2 (MQ=255)            gAAAAAATAAAATCTGGGATTTCACCGGCGGCATCCCTCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGGGCACCTGGCGCAgg                                                                                   <  2:69300/99‑2 (MQ=255)            gAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCg                                                                                  >  2:519683/1‑100 (MQ=255)            gAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGAGGAAAACCCAGTCCAACGGTACACCACTGCGCAAGGTGCCCCTGGCGCAGCg                                                                                  >  2:763202/1‑100 (MQ=255)             aaaaaaTAAAATCTGGGATTTCAACGGCGGCATCCATCCCCCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGGGCACCTGGCGCAgggg                                                                                 <  2:2022100/100‑4 (MQ=255)              aaaaaTAAAGCCTGGGAATGCAAAGGCGCGATCCGTTCACTCGAGAGGCGCCCCCCATCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGt                                                                                 <  2:314116/99‑1 (MQ=255)               aaaataaaatCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGttt                                                                               >  1:1945283/1‑100 (MQ=255)                aaataaaatCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGtttt                                                                              <  1:549805/100‑1 (MQ=255)                 aataaaatCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTg                                                                             >  1:54533/1‑100 (MQ=255)                  ataaaatCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGt                                                                            <  1:1948689/100‑1 (MQ=255)                  ataaaatCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGt                                                                            <  1:1965361/100‑1 (MQ=255)                   taaaatCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGtt                                                                           <  1:1614511/100‑1 (MQ=255)                   taaaatCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGtt                                                                           <  1:1379366/100‑1 (MQ=255)                     aaatCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTACAACGGTACACCCCTGCGCCGGGTGCCACTGGCGCAGCGgttttt                                                                            >  2:1098250/1‑95 (MQ=255)                      aatCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTAtt                                                                        <  1:347943/100‑1 (MQ=255)                       atCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTc                                                                       >  1:486898/1‑100 (MQ=255)                        tCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTcc                                                                      <  2:1924403/100‑1 (MQ=255)                        tCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTcc                                                                      <  1:1066313/100‑1 (MQ=255)                        tCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTcc                                                                      <  1:1074514/100‑1 (MQ=255)                         cTGGAATTACAACGACGAAATCCATCCCCCGGAGCTGAAATCCCACTCCACCGGGTCTACCCCGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTcc                                                                      <  2:973614/99‑1 (MQ=255)                          tGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCAc                                                                    <  1:410575/100‑1 (MQ=255)                           gggATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACt                                                                   >  1:1957944/1‑100 (MQ=255)                            ggATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTg                                                                  <  1:52510/100‑1 (MQ=255)                              aTTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGaa                                                                <  1:944934/100‑1 (MQ=255)                                ttCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAAc                                                              <  1:918097/100‑1 (MQ=255)                                ttCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAAc                                                              <  1:1594053/100‑1 (MQ=255)                                 tCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGCTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAAc                                                              >  2:482910/1‑99 (MQ=255)                                  cAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACa                                                             >  2:2517957/1‑99 (MQ=255)                                  cAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAg                                                            <  1:413655/100‑1 (MQ=255)                                  cAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAg                                                            >  1:406824/1‑100 (MQ=255)                                    aaggcggcATCCAGCCACCGGAGATGAAAACACAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCa                                                          >  1:1727563/3‑100 (MQ=255)                                     cggcggCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCat                                                         >  1:1439098/1‑100 (MQ=255)                                     cggcggCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCa                                                          <  2:1412879/99‑1 (MQ=255)                                         ggcATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTg                                                     <  1:1344568/100‑1 (MQ=255)                                           catccatcCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGgc                                                   >  1:559754/1‑100 (MQ=255)                                           catccatcCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGgc                                                   <  1:2036365/100‑1 (MQ=255)                                             tccatccaCCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGgcgc                                                 >  2:1659442/1‑100 (MQ=255)                                               catccaCCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGTTGCCCCTGGCGCACCGTCTTGTTATACCACTGAAACAGCAGATTGgcgc                                                 >  2:1325832/1‑98 (MQ=255)                                                atccaCCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGa                                              <  1:327956/100‑1 (MQ=255)                                                 tccaCCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGTCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGaa                                             <  2:420403/100‑1 (MQ=255)                                                  ccaCCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGaa                                             <  1:734679/99‑1 (MQ=255)                                                    aCCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAg                                            <  1:1614731/98‑1 (MQ=255)                                                     ccGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGt                                          <  1:515653/99‑1 (MQ=255)                                                     ccGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTg                                         >  1:2063743/1‑100 (MQ=255)                                                     ccGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTg                                         >  2:222823/1‑100 (MQ=255)                                                     ccGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTg                                         >  1:2463001/1‑100 (MQ=255)                                                       ggAGATGAAACCCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAg                                       <  1:2289200/100‑1 (MQ=255)                                                        gagaTGAAAACCCAGTCCAACGGTACACCTCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGt                                      <  1:158866/100‑1 (MQ=255)                                                         agaTGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGtt                                     <  2:2308637/100‑1 (MQ=255)                                                            tGAAAACCCAGTCTAACGGTACAGCCATGCGCCTGGTGCCCCTGTCGCAGCGTTTTGTTATTCCACTGAAAGAGCATATTGGCGCTGAAGGTGAGTtg                                    >  2:376494/1‑98 (MQ=255)                                                            tGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTtgtg                                  <  1:1256382/100‑1 (MQ=255)                                                            tGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTtgtg                                  <  1:183699/100‑1 (MQ=255)                                                             gAACACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGATATTCCCCTGAAACAGCATATTGGCGCTGAAGGTGAGTTATGc                                 >  2:2548301/1‑100 (MQ=255)                                                             gAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGc                                 <  1:92005/100‑1 (MQ=255)                                                              aaaaCCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTtgt                                   >  1:1905699/1‑97 (MQ=255)                                                              aaaaCCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCg                                <  2:391837/100‑1 (MQ=255)                                                               aaaCCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGt                               >  2:74520/1‑100 (MQ=255)                                                                aaCCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGtt                              >  2:718565/1‑100 (MQ=255)                                                                 aCCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTa                             >  1:1781735/1‑100 (MQ=255)                                                                 aCCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTa                             >  1:2120682/1‑100 (MQ=255)                                                                    cAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCg                          <  1:1659828/100‑1 (MQ=255)                                                                       tCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGATTTGTTATTCCACTGAAAAAGCATATTGGCGCTGAAGGTGAGGTGTGCGTTAGCGTCg                       >  1:2515908/1‑100 (MQ=255)                                                                        ccAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTCTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCg                       >  2:583994/1‑99 (MQ=255)                                                                         cAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGc                     <  2:4981/100‑1 (MQ=255)                                                                            cGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGAt                  >  1:2146179/1‑100 (MQ=255)                                                                              acaaaacgcTGCGCCAGGGGCACCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATaa                >  1:2363613/9‑100 (MQ=38)                                                                                acacCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAg              <  1:1659454/100‑1 (MQ=255)                                                                                  acCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTa            >  1:1317843/1‑100 (MQ=255)                                                                                   ccccTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTAt           <  2:2018640/100‑1 (MQ=255)                                                                                    cgcTGCGCCAGGGGCACCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTAtt          >  1:1418002/3‑100 (MQ=255)                                                                                    cgcTGCGCCAGGGGCACCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCAGCGATAAAGTAtt          >  2:2521941/3‑100 (MQ=255)                                                                                    cgcTGCGCCAGGGGCACCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGATGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTAtt          >  2:360760/3‑100 (MQ=255)                                                                                    cccTGCGCCAGGTGCCCCTTGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTTAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTa            <  2:1571552/98‑1 (MQ=255)                                                                                     gcTGCGCCAGGGGCACCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTg         >  1:957412/2‑100 (MQ=255)                                                                                     gcTGCGCCAGGGGCACCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTg         >  2:820681/2‑100 (MQ=255)                                                                                     ccTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTg         <  2:2096532/100‑1 (MQ=255)                                                                                      cTGCGCCAGGGGCACCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTgc        >  1:449572/1‑100 (MQ=255)                                                                                       tGCGCCAGGTGCCGCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTgcg       >  1:1721838/1‑100 (MQ=255)                                                                                          gcCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCgg    <  2:682201/100‑1 (MQ=255)                                                                                           cctggTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAACCAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGc   <  1:1325841/97‑1 (MQ=255)                                                                                           ccAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCgg    <  2:1791288/99‑1 (MQ=255)                                                                                            cAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGcc  <  1:718572/100‑1 (MQ=255)                                                                                                 |                                                                                               CTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCC  >  NC_000913/1706933‑1707123

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