Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,409,134	T→A	6.3%	noncoding (6/122 nt)	fnrS →	FNR‑activated anaerobic sRNA antisense regulator down regulates cydDC, metE, sodA, sodB, Hfq‑dependent; mediates negative FNR regulation

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,409,134	0	T	A	6.3%	161.9 / 2.6	64	noncoding (6/122 nt)	fnrS	FNR‑activated anaerobic sRNA antisense regulator down regulates cydDC, metE, sodA, sodB, Hfq‑dependent; mediates negative FNR regulation
Reads supporting (aligned to +/- strand):  ref base T (31/29);  new base A (2/2);  total (33/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 = 2.72e-01

TTATTGGTGGTGTTGCTTTATGGTTGCATCGTAGTAAATGGTTGTAACAAAAGCAATTTTTCCGGCTGTCTGTATACAAAAACGCCGCAAAGTTTGAGCGAAGTCAATAAACTCTCTACCCATTCAGGGCAATATCTCTCTTGCAGGTGAATGCAACGTCAAGCGATGGGCGTTGCGCTCCATATTGTCTTACTTCCTTTTTTGAATTACTGCATAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGCTTTTTTTTGCCTGTTATTTATCAGCGTCTACCCT  >  NC_000913/1408987‑1409276                                                                                                                                                    |                                                                                                                                               ttATTGGTGGTGTTGCTTTATGGTTGCATCGTAGTAAATGGTTGTAACAAAAGCAATTTTTACGGCTGTCTGTATACAAAAACTCCGCAAAGTTTGAGCGAAGTCAATAAACTCTCTACCCATTCAGGGCAATATCTCTATTGCAGCTg                                                                                                                                               >  1:70974/1‑149 (MQ=255)          gtgtTGCTTTATGGTTGCATCGTAGTAAATGGTTGTAACAAAAGCAATTTTTCCGGCTGTCTGTATACAAAAACGCCGCAAAGTTTGAGCGAAGTCAATAAACTCTCTACCCATTCAGGGCAATATCTCTCTTGCAGGTGAATGCAACg                                                                                                                                      <  1:1127141/149‑1 (MQ=255)          gtgtTGCTTTATGGTTGCATCGTAGTAAATGGTTGTAACAAAAGCAATTTTTCCGGCTGTCTGTATACAAAAACGCCGCAAAGTTTGAGCGAAGTCAATAAACTCTCTACCCATTCAGGGCAATATCTCTCTTGCAGGTGAATGCAACg                                                                                                                                      >  2:606292/1‑149 (MQ=255)             tttcTTTATGGTTGCATCGTAGTAAATGGTTGTAACAAAAGAAATTTTTCCGGCTGTCTTTATACAAAAACGCCGCAAAGTTTGAGCGAAGTCAATAAACTCTCTACCCATTCAGGGCAATATCTCTCTTGCAGGTGAATGCAAAGTct                                                                                                                                   <  2:579195/146‑2 (MQ=255)             ttGCTTTATGGTTGCATCGTAGTAAATGGTTGTAACAAAAGCAATTTTTCCGGCTGTCTGTATACAAAAACGCCGCAAAGTTTGAGCGAAGTCAATAAACTCTCTACCCATTCAGGGCAATATCTCTCTTGCAGGTGAATGCAACGTCa                                                                                                                                   >  1:388738/1‑149 (MQ=255)               gCTTTATGGTTGCATCGTAGTAAATGGTTGTAACAAAAGCAATTTTTCCGGCTGTCTGTATACAAAAACGCCGCAAAGTTTGAGCGAAGTCAATAAACTCTCTACCCATTCAGGGCAATATCTCTCTTGCAGGTGAATGCAACGTCAAg                                                                                                                                 >  2:249613/1‑149 (MQ=255)                 tttATGGTTGCATCGTAGTAAATGGTTGTAACAAAAGCAATTTTTCCGGCTGTCTGTATACAAAAACGCCGCAAAGTTTGAGCGAAGTCAATAAACTCTCTACCCATTCAGGGCAATATCTCTCTTGCAGGTGAATGCAACGTCAAGCg                                                                                                                               >  1:968258/1‑149 (MQ=255)                   tATGGTTGCATCGTAGTAAATGGTTGTAACAAAAGCAATTTTTCCGGCTGTCTGTATACAAAAACGCCGCAAAGTTTGAGCGAAGTCAATAAACTCTCTACCCATTCAGGGCAATATCTCTCTTGCAGGTGAATGCAACGTCAAGCGAt                                                                                                                             <  1:981341/149‑1 (MQ=255)                   tATGGTTGCATCGTAGTAAATGGTTGTAACAAAAGCAATTTTTCCGGCTGTCTGTATACAAAAACGCCGCAAAGTTTGAGCGAAGTCAATAAACTCTCTACCCATTCAGGGCAATATCTCTCTTGCAGGTGAATGCAACGTCAAGCGAt                                                                                                                             <  1:761350/149‑1 (MQ=255)                        ttGCATCGTAGTAAATGGTTGTAACAAAAGCAATTTTTCCGGCTGTCTGTATACAAAAACGCCGCAAAGTTTGAGCGAAGTCAATAAACTATCTACCCATTCAGGGCAATATCTCTCTTGCAGGTGAATGCAACGTCAAGCGATGGGCg                                                                                                                        >  1:247193/1‑149 (MQ=255)                                tagtaAATGGTTGTAACAAAAGCAATTTTTCCGGCTGTCTGTATACAAAAACGCCGCAAAGTTTGAGCGAAGTCAATAAACTCTCTACCCATTCAGGGCAATATCTCTCTTGCAGGTGAATGCAACGTCAAGCGATGGGCGTTGCGCTc                                                                                                                >  2:976002/1‑149 (MQ=255)                                  gtaAATGGTTGTAACAAAAGCAATTTTTCCGGCTGTCTGTATACAAAAACGCCGCAAAGTTTGAGCGAAGTCAATAAACTCTCTACCCATTCAGGGCAATATCTCTCTTGCAGGTGAATGCAACGTCAAGCGATGGGCTTTGCGCTcaa                                                                                                              >  2:664382/1‑147 (MQ=255)                                    aaaTGGTTGTAACAAAAGCAATTTTTCCGGCTGTCTGTATACAAAAACGCCGGAAAGTTTGAGCGAAGTCAATAAAGTCTCTACCCATTCAGGGCAATATGTCTCTTGCAGGTGAATGCAACGTCAAGCGATGGGCGAGGAGCGCCata                                                                                                            >  2:526448/1‑149 (MQ=255)                                     aaTGGTTGTAACAAAAGCAATTTTTCCGGCTGTCTGTATACAAAACCGCCGCAAAGTTTGAGCGAAGTCAATAAACTCTCTACCCATTCAGGGCAATATCTCTCTTGCAGGTGAATGCAACGTCAAGCGATGGGCGTTGCGCTCCatat                                                                                                           >  1:592002/1‑149 (MQ=255)                                         gTTGTAACAAAAGCAATTTTTCCGGCTGTCTGTATACAAAAACGCCGCAAAGTTTGAGCGAAGTCAATAAACTCTCTACCCATTCAGGGCAATATCTCTCTTGCAGGTGAATGCAACGTCAAGCGATGGGCGTTGCGCTCCATATTGTc                                                                                                       <  2:259818/149‑1 (MQ=255)                                                cAAAAGCAATTTTTCCGGCTGTCTGTATACAAAAACGCCGCAAAGTTAGAGCGATGTCAATTAACTCTCTACCCAGTGACCGCCATAGCTCTCTAGCTGGTGAATGGAACGTCAAGCGATGGCCGTGGCGCTCCAGATTGTcttccttc                                                                                                >  2:921364/1‑149 (MQ=255)                                                  aaaGCAATTTTTCCGGCTGTCTGTATACAAAAACGCCGCAAAGGGTGAGCGAAGTCAATAAACTCTCTACCCAATCAGGGCAATATCTCTCTTGAAGGTGAATGCAACGTCAAGCGATGGGCGTTGAGCTCCATAACGTATTACTTcat                                                                                              >  1:680958/1‑147 (MQ=255)                                                   aaGCAATTTTTCCGGCTGTCTGTATACAAAAACGCCGCAAAGTTTGAGCGAAGTCAATAAACTCTCTACCCATTCAGGGCAATATCTCTCTTGCAGGTGAATGCAACGTCAAGCGATGGGCGTTGCGCTCCATATTGTCTTACTTCCtt                                                                                             >  2:554644/1‑149 (MQ=255)                                                      cAATTTTTCCGGCCGTCTGTATACAAACACGCCGCAAAGTTTGAGCGAAGTCAATAAACTCTCTACCCATTCAGGGCAATATCTCTCTTGCAGGTGAATGCAACGTCAAGCGATGGGCGTTGCGCTCCATATTGTCTTACTTCCttttt                                                                                          <  2:372191/149‑1 (MQ=255)                                                       aaTTTTTCCGGCTGTCTGTATACAAAAACGCCGCAAAGTTTGAGCGAAGTCAATAAACTCTCTACCCATTCAGGGCAATATCTCTCTTGCAGGTGATTGCAACGTCAAGCGATGGGCGTTGCGCTCCATATTGTCTTACTTCCtttttt                                                                                         <  1:664382/149‑1 (MQ=255)                                                       aaTTTTTCCGGCTGTCTGTATACAAAAACGCCGCAAAGTTTGAGCGAAGTCAATAAACTCTCTACCCATTCAGGGCAATATCTCTCTTGCAGGTGAATGCAACGTCAAGCGATGGGCGTTGCGCTCCATATTGTCTTACTTCCtttttt                                                                                         <  2:981970/149‑1 (MQ=255)                                                       aaTTTTTCCGGCTGTCTGTATACAAAAACGCCGCAAAGTTTGAGCGAAGTCAATAAACTCTCTACCCATTCAGGGCAATATCTCTCTTGCAGGTGAATGCAACGTCAAGCGATGGGCGTTGCGCTCCATATTGTCTTACTTCCtttttt                                                                                         <  2:500122/149‑1 (MQ=255)                                                        aTTTTTCCGGCTGTCTGTATACAAAAACGCCGCAAAGTTTGAGCGAAGTCAATAAACTCTCTACCCATTCAGGGCACTATCTCTCTTGCAGGTGAATGCAACGTCAAGCGATGGGCGTTGCGCTCCATATTGGCTTACTTCCTTTTTTg                                                                                        <  2:700954/149‑1 (MQ=255)                                                        aTTTTTCCGGCTGTCTGAATACAAAAACGCAGCAAAGTTTGAGCGAATACAATAAACTCTCTACCCATTCAGGGAAATATATCTCTTGCAGCTGAATGCAACGTCCAGCAATGGTCGTTGCTCTACATATTTTATTACTTACttatttt                                                                                        >  2:415322/1‑148 (MQ=255)                                                          ttttCCGGCTGTCTGTATACAAAAACGCCGCAAAGTTTGAGCGAAGTCAATAAACTCTCTACCCATTCAGGGCAATATCTCTCTTGCAGGTGAATGCAACGTCAAGCGATGGGCGTTGCGCTCCATATTGTCTTACTTCCTTTTTTGaa                                                                                      <  1:387419/149‑1 (MQ=255)                                                            ttCCGGCTGTCTGTATACAAAAACGCCGCAAAGTTTGAGCGAAGTCAATAAACTCTCTACCCATTCAGGGCAATATCTCTCTTGCAGGTGAATGCAACGTCAAGCGATGGGCGTTGCGCTCCATATTGTCTTACTTCCTTTTTTGAAtt                                                                                    <  2:143156/149‑1 (MQ=255)                                                               cGGCTGTCTGTATACAAAAACGCCGCAAAGTTTGAGCGAAGTCAATAAACTCTCTACCCATTCAGGGCAATATCTCTCTTGCAGGTGAATGCAACGTCAAGCGATGGGCGTTGCGCTCCATATTGTCTTACTTCCTTTTTTGAATTACt                                                                                 >  1:586997/1‑149 (MQ=255)                                                               cGGCTGTCTGTATACAAAAACGCCGCAAAGTTTGAGCGAAGTCAATAAACTCTCTACCCATTCAGGGCAATATCTCTCTTGCAGGTGAATGCAACGTCAAGCGATGGGCGTTGCGCTCCATATTGTCTTACTTCCTTTTTTGAATTACt                                                                                 <  2:388738/149‑1 (MQ=255)                                                                    gtcttgataCAAAAACGCCGCACAGTTTGAGGGAAGTCAATAAACTCTCTACCCATTCAGGGCAATATCTCTCTTGCAGGTGAATGCAACGTCAAGCGATGGGCGTTGCGCTCCATATTGTCTTACTTCCTTTTTTGAATTACTGCTTa                                                                            <  2:844766/143‑1 (MQ=255)                                                                        gtATACAAAAACGCCGCAAAGTTTGAGCGAAGTCAATAAACTCTCTACCCATTCAGGGCAATATCTCTCTTGCAGGAGAATGCAACGCCCATCGCTTGACGTTgca                                                                                                                   >  1:173661/1‑105 (MQ=255)                                                                        gtATACAAAAACGCCGCAAAGTTTGAGCGAAGTCAATAAACTCTCTACCCATTCAGGGCAATATCTCTCTTGCAGGAGAATGCAACGCCCATCGCTTGACGTTgca                                                                                                                   <  2:173664/106‑2 (MQ=255)                                                                        gtATACAAAAACGCCGCAAAGTTTGAGCGAAGTCAATAAACTCTCTACCCATTCAGGGCAATATCTCTCTTGCAGGAGAATGCAACGCCCATCGCTTGACGTTgca                                                                                                                   <  2:173661/106‑2 (MQ=255)                                                                        gtATACAAAAACGCCGCAAAGTTTGAGCGAAGTAAATAAACTCTCTACCCATTCAGGGCAATATCTCTCTTGCAGGAGAATGCAACGCCCATCGCTTGACGTTgca                                                                                                                   >  1:173664/1‑105 (MQ=255)                                                                          ataCAAAAACGCCGCAAAGTTTGAGCGAAGTCAATAAACTCTCTACCCATTCAGGGCAATATCTCTCTTGCAGGTGAATGCAACGTCAAGCGATGGGCGTTGCGCTCCATATTGTCTTACTTCCTTTTTTGAATTACTGCATAGCACaa                                                                      >  1:597714/1‑149 (MQ=255)                                                                          ataCAAAAACGCCGCAAAGTTTGAGCGAAGTCAATAAACTCTCTACCCATTCAGGGCAATATCTCTCTTGCAGGTGAATGCAACGACAAGCGATGGGCGTTGCGCTCCATATTGTCTTACTTACTTTTTTGAATTACTGCATAGCACaa                                                                      >  1:993226/1‑149 (MQ=255)                                                                            aCAAAAACGCCGCAAAGTTTGAGCGAAGTCAATAAACTCTCTACCCATTCAGGGCAATATCTCTCTTGCAGGTGAATGCAACGTCAAGCGATGGGCGTTGCGCTCCATATTGTCTTACTTCCTTTTTTGAATTACTGCATAGCACAAtt                                                                    <  1:1047293/149‑1 (MQ=255)                                                                              aaaaaCGCCGCAAAGTTTGAGCGAAGTCAATAAACTCTCTACCCATTCAGGGCAATATCTCTCTTGCAGGTGAATGCAACGTCAAGCGATGGGCGTTGCGCTCCATATTGTCTTACTTCCTTTTTTGAATTACTGCATAGCACAATTGa                                                                  >  1:412277/1‑149 (MQ=255)                                                                               aaaaCGCCGCAAAGTTTGAGCGAAGTCAATAAACTCTCTACCCATTCAGGGCAATATCTCTCTTGCAGGTGAATGCAACGTCAAGCGATGGGCGTTGCGCTCCATATTGTCTTACTTCCTTTTTTGAATTACTGCATAGCACAATTGAt                                                                 <  1:254832/149‑1 (MQ=255)                                                                                 aaCGCCGCAAAGTTTGAGCGAAGTCAATAAACTCTCTACCCATTCAGGGCAATATCTCTCTTGCAGGTGAATGCAACGTCAAGCGATGGGCGTTGCGCTCCATATTGTCTTACTTCCTTTTTTGAATTACTGCATAGCACAATTGATTc                                                               >  2:1121967/1‑149 (MQ=255)                                                                                    gccgcAAAGTTTGAGCGAAGTCAATAAACTCTCTACCCATTCAGGGCAATATCTCTCTTGCAGGTGAATGCAACGTCAAGCGATGGGCGTTGCGCTCCATATTGTCTTACTTCCTTTTTTGAATTACTGCATAGCACAATTGATTCGTa                                                            <  1:976002/149‑1 (MQ=255)                                                                                     ccgaAAAGTTTGAGCGAAGTCAATAAACTCTCTACCCATTCAGGGCAATATCTCTCTTTCATGTGAATGCAACGTCAATCGATTGGCGTTGCGATCCATATTGTCTTACTTCCTTTTTTGAATTACTTCATCGCACAATTGATTCGTac                                                           <  2:135053/149‑1 (MQ=255)                                                                                       gcAAAGTTTGAGCGAAGTCAATAAACTCTCTACCCATTCAGGGCAATATCTCTCTTGCAGGTGAATGCAACGTCAAGCGATGGGCGTTGCGCTCCATATTGTCTTACTTCCTTTTTTGAATTACTGCATAGCACAATTGATTCGTacga                                                         <  2:730526/149‑1 (MQ=255)                                                                                        cAAAGTTTGAGCGAAGTCAATAAACTCTCTACCCATTCAGGGCAATATCTCTCTTGCAGGTGAATGCAACGTCAAGCGATGGGCGTTGCGCTCCATATTGTCTTACTTCCTTTTTTGAATTACTGCATAGCACAATTGATTCGTacgac                                                        >  1:995708/1‑149 (MQ=255)                                                                                             tttGAGCGAAGTCAATAAACTCTCTACCCATTCAGGGCAATATCTCTCTTGCAGGTGAATGCAACGTCAAGCGATGGGCGTTGCGCTCCATATTGTCTTACTTCCTTTTTTGAATTACTGCATAGCACAATTGATTCGTACGACGCCgc                                                   >  1:816482/1‑148 (MQ=255)                                                                                                   cGAAGTCAATAAACTCTCTACCCATTCAGGGCAATATCTCTCTTGCAGGTGAATGCAACGTCAAGCGATGGGCGTTGCGCTCCATATTGTCTTACTTCCTTTTTTGAATTACTGCATAGCACAATTGATTCGTACGACGCCGACTTtga                                             <  2:968258/149‑1 (MQ=255)                                                                                                    gAAGTCAATAAACTCTCTACCCATTCAGGGCAATATCTCTCTTGCAGGTGAATGCAACGTCAAGCGATGGGCGTTGCGCTCCATATTGTCTTACTTCCTTTTTTGAATTACTGCCTCGCACAATTGATTCGTACGACGCCGACTTtgat                                            <  2:336505/149‑1 (MQ=255)                                                                                                       gTCAATAAACTCTCTACCCATTCAGGGCAATATCTCTCTTGCAGGTGAATGCAACGTCAAGCGATGGGCGTTGCGCTCCATATTGTCTTACTTCCTTTTTTGAATTACTGCATAGCACAATTGATTCGTACGACGCCGACTTTGATGAg                                         >  1:776794/1‑149 (MQ=255)                                                                                                        tCAATAAACTCTCTACCCATTCAGGGCAATATCTCTCTTGCAGGTGAATGCAACGTCAAGCGATGGGCGTTGCGCTCCATATTGTCTTACTTCCTTTTTTGAATTACTGCATAGCACAATTGATTCGTACGACGCCGACTTTGATGAGt                                        >  1:611767/1‑149 (MQ=255)                                                                                                          aaTAAACTCTCTACCCATTCAGGGCAATATCTCTCTTGCAGGTGAATGCAACGTCAAGCGATGGGCGTTGCGCTCCATATTGTCTTACTTCCTTTTTTGAATTACTGCATAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCg                                      >  1:297431/1‑149 (MQ=255)                                                                                                              aaCTCTCTACCCATTCAGGGCAATATCTCTCTTGCAGGTGAATGCAACGTCAAGCGATGGGCGTTGCGCTCCATTTTGTCTTACTTCCTTTTTTGAATTACTGCATAGCACAATTGATTCGTACGACGCCGACTTTGAGGAGTCGGCtt                                  >  2:185588/1‑149 (MQ=255)                                                                                                                  ctgtACCCATTCAGGGAAATATCTCTCTTGGAGGTTAATGCAACGTGAAGCGATGGGCGTTGCGCTTGATATTGTCTTACTTCCTTTTTTGAATTTCTGCATAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGCtttttt                              <  2:58940/146‑1 (MQ=255)                                                                                                                  ctctACCCATTCAGGGCAATATCTCTCTTGCAGGTGAATGCAACGTCAAGCGATGGGCGTTGCGCTCCATATTGTCTTACTTCCTTTTTTGAATTACTGCATAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGCtttttt                              <  1:368916/149‑1 (MQ=255)                                                                                                                   tctACCCATTCAGGGCAATATCTCTCTTGCAGGTGAATGCAACGTCAAGCGATGGGCGTTGCGCTCCATATTGTCTTACTTCCTTTTTTGAATTACTGCATAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGCttttttt                             <  1:554644/149‑1 (MQ=255)                                                                                                                    ctACCCATTCATGGCAATATCTCTCTTGCCGGTGACAGCAACGTCAAGCGATGGGCGTTGCGCTGCATATTGTCTTACTTCCTTTTTTGAATTACTGCATAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGCtttttttt                            <  1:689623/149‑1 (MQ=255)                                                                                                                      aCCCATTCAGGGCAATATCTCTCTTGCAGGTGAATGCAACGTCAAGCGATGGGCGTTGCGCTCCATATTGTCTTACTTCCTTTTTTGAATTACTGCATAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGCTTTTTTTTGc                          >  2:237654/1‑149 (MQ=255)                                                                                                                      aCCCATTCAGGGCAATATCTCTCTTGCAGGTGAATGCAACGGCAAGCGATGGGCGTTGCGCTCCATATTGTCTTACTTCCTTTTTTGAATTACTGCATAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGCTTTTTTTTGc                          <  2:237524/149‑1 (MQ=255)                                                                                                                        ccATTCAGGGCAATATCTCTCTTGCAGGTGAATGCAACGTCAAGCGATGGGCGTTGCGCTCCATATTGTCTTACTTCCTTTTTTGAATTACTGCATAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGCTTTTTTTTGCCt                        >  1:64784/1‑149 (MQ=255)                                                                                                                             cAGGGCAATATCTCTCTTGCAGGTGAATGCAACGTCAAGCGTTGGTCGTTGCGCTCCATATTGTATTACTTCCTTTTTTGAATTACTGCATAGCACAATTGATTCGTACGTCGCCGACTTTGATGAGTCGGCTTTTTTTTGCCTGttat                   <  2:72698/149‑1 (MQ=255)                                                                                                                              aGGTCAATATCTCTCTTGCAGGTGAATGCAACGTCAAGCGATGGGCGTTGCGCTCCATATTGTCTTACTTCCTTTTTTGAATTACTGCATAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGC‑TTTTTTTGCCTGttattt                 <  1:526448/149‑1 (MQ=255)                                                                                                                                 gCAATATCTCTCTTGCAGGTGAATGCAACGTCAAGCGATGGGCGTTGCGCTCCATATTGTCTTACTTCCTTTTTTGAATTACTGCATAGCACAATTGATTCGTACGACTCCGACTTTGATGACTCGGCTGTTTTTTGCCAGttatttat               >  2:896070/1‑149 (MQ=255)                                                                                                                                       tctctcTTGCAGGTGAATGCAACGTCTAGCGATGGGCGTTGCGCTCCATATTGTCTTACTTGCTTGTTTGAATTACTGGATAGCCCAATTGATTCGTACGACGCCGAGTGTGATGAGTCGGCTTTTTTTGGCCTGTTATTTATGAGCgg         >  2:138748/1‑148 (MQ=255)                                                                                                                                            ctTGCAGGTGAATGCAACGTCAAGCGATGGGCGTTGCGCTCCATATTGTCTTACTTCCTTTTTTGAATTACTGCATAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGCTTTTTTTTGCCTGTTATTTATCAGCGTCTAcc    >  1:1119677/1‑149 (MQ=255)                                                                                                                                             ttGCAGGTGAATGCAACGTCAAGCGATGGGCGTTGCGCTCCATATTGTCTTACTTCCTTTTTTGAATTACTGCATAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGCTTTTTTTTGCCTGTTATTTATCAGCGTCTAccc   >  2:536331/1‑149 (MQ=255)                                                                                                                                              tGCAGGTGAATGCTACGTCAAGCGATGGGCGTTGCGCTCCATATTGTCTTACTTACTTTTTTGAATTACTGCATAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGCTTTTTTTTGCCTGTTATTTATCAGCGTCTACCCt  <  1:273807/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                               TTATTGGTGGTGTTGCTTTATGGTTGCATCGTAGTAAATGGTTGTAACAAAAGCAATTTTTCCGGCTGTCTGTATACAAAAACGCCGCAAAGTTTGAGCGAAGTCAATAAACTCTCTACCCATTCAGGGCAATATCTCTCTTGCAGGTGAATGCAACGTCAAGCGATGGGCGTTGCGCTCCATATTGTCTTACTTCCTTTTTTGAATTACTGCATAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGCTTTTTTTTGCCTGTTATTTATCAGCGTCTACCCT  >  NC_000913/1408987‑1409276

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑