Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,233,710	A→G	6.1%	G336G (GGT→GGC)	nhaB ←	sodium:proton antiporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,233,710	0	A	G	6.1%	201.1 / 3.5	65	G336G (GGT→GGC)	nhaB	sodium:proton antiporter
Reads supporting (aligned to +/- strand):  ref base A (31/30);  new base G (2/2);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.11e-01

TCCGATGCCTGCAACACAAACTGAATAATTGGCGAAAACAGTTGTTGGTCGATAATCACCGCGACTACCGAGAAAAACACCGTCAACAGTGCGGTGAATGGCAGAGATTCGGTGAAGGCTTTACCGATAGCATGCTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACAATCAGACGGATTTTATCCTGACGGGTACGCTGGTG  >  NC_000913/1233564‑1233853                                                                                                                                                   |                                                                                                                                                tCCGATGCCTGCAACACAAACTGAATAATTGGCGAAAACAGTTGTTGGTCGATAATCACCGCGACTACCGAGAAAAACACCGTCAACAGTGCGGTGAATGGCAGAGATTCGGTGAAGGCTTTACCGATAGCATGCTCATCGGTGacac                                                                                                                                                >  5:291948/1‑148 (MQ=255)      tGCCTGCAACACAAACTGAATAATTGGCGAAAACAGTTGTTGGTCGATAATCACCGCGACTACCGAGAAAAACACCGTCAACAGTGCGGTGAATGGCAGAGATTCGGTGAAGGCTTTACCGATAGCATGCTCATCGGTGACACCGGt                                                                                                                                            >  6:67905/1‑147 (MQ=255)         cTGCAACACAAACTGAATAATTGGCGAAAACAGTTGTTGGTCGATAATCACCGCGACTACCGAGAAAAACACCGTCAACAGTGCGGTGAATGGCAGAGATTCGGTGAAGGCTTTACCGATAGCATGCTCATCGGTGACACCGGTCAat                                                                                                                                        <  8:125194/148‑1 (MQ=255)           gCAACACAAACTGAATAATTGGCGAAAACAGTTGTTGGTCGATAATCCCCGCGACTACCGAGAAAAACACCGTCAACAGTGCGGTGAATGGCAGAGATTCGGTGAAGGCTTTACCGATAGCATGCTCATCGGTGACACCGGTCAatga                                                                                                                                      >  1:288257/1‑148 (MQ=255)            cAACACAAACTGAATAATTGGCGAAAACAGTTGTTGGTCGATAATCACCGCGACTACCGAGAAAAACACCGTCAACAGTGCGGTGAATGGCAGAGATTCGGTGAAGGCTTTACCGATAGCATGCTCATCGGTGACACCGGTCAatgat                                                                                                                                     <  5:224643/148‑1 (MQ=255)             aaCACAAACTGAATAATTGGCGAAAACAGTTGTTGGTCGATAATCACCGCGACTACCGAGAAAAACACCGTCAACAGTGCGGTGAATGGCAGAGATTCGGTGAAGGCTTTACCGATAGCATGCTCATCGGTGACACCGGTCAatgatg                                                                                                                                    <  6:230549/148‑1 (MQ=255)                    aCTGAATAATTGGCGAAAACAGTTGTTGGTCTATAATCACCGCGACTACCGAGAAAAACACCGTCAACAGTGAGGTGAATTGCAGATATTCGGTGAAGGCTTTACCGATAGCATGCTCATCGGTGACACCGGTCAATGATGTTGc                                                                                                                                >  3:261668/1‑145 (MQ=255)                       gAATAATTGGCGAAAACAGTTGTTGGTCGATAATCACCGCGACTACCGAGAAAAACACCGTCAACAGTGCGGTGAATGGCAGAGATTCGGTGAAGGCTTTACCGATAGCATGCTCATCGGTGACACCGGTCAATGATGATGCCAGaa                                                                                                                           >  7:7594/1‑147 (MQ=255)                         ataatTGGCGAAAACAGTTGTTGGTCGATAATCACCGCGACTACCGAGAAAAACACCGTCAACAGTGCGGTGAATGGCAGAGATTCGGTGAAGGCTTTACCGATAGCATGCTCATCGGTGACACCGGTCaa                                                                                                                                         <  1:298816/131‑1 (MQ=255)                         ataatTGGCGAAAACAGTTGTTGGTCGATAATCACCGCGACTACCGAGAAAAACACCGTCAACAGTGCGGTGAATGGCAGAGATTCGGTGAAGGCTTTACCGATAGCATGCGCATCGGTGACACCGGTCaa                                                                                                                                         >  2:298816/1‑131 (MQ=255)                              tGGCGAAAACAGTTGTTGGTCGATAATCACCGCGACTACCGAGAAAAACACCGTCAACAGTGCGGTGAATGGCAGAGATTCGGTGAAGGCTTTACCGATAGCATGCTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAAtgact                                                                                                                   >  7:305983/1‑148 (MQ=255)                              tGGCGAAAACAGTTGTTGGTCGATAATCACCGCGACTACCGAGAAAAACACCGTCAACAGTGCGGTGAATGGCAGAGATTCGGTGAAGGCTTTACCGATAGCATGCTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAAtgact                                                                                                                   >  3:270658/1‑148 (MQ=255)                              tGGCGAAAACAGTTGTTGGTCGATAATCACCGCGACTACCGAGAAAAACACCGTCAACAGTGCGGTGAATGGCAGAGATTCGGTGAAGGCTTTACCGATAGCATGCTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAAtgact                                                                                                                   >  7:300416/1‑148 (MQ=255)                                 cGAAAACAGTTGTTGGTCGATAATCACCGCGACTACCGAGAAAAACACCGTCAACAGTGCGGTGAATGGCAGAGATTCGGTGAAGGCTTTACCGATAGCATGCTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAAtgactga                                                                                                                 >  5:244564/1‑147 (MQ=255)                                              tgGTCGATAATCACCGCGACTACCGAGAAAAACACCGTCAACAGTGCGGTGAATGGCAGAGATTCGGTGAAGGCTTTACCGATAGCATGCTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGc                                                                                                   >  8:184336/1‑148 (MQ=255)                                                gTCGATAATCACCGCGACTACCGAGAAAAACACCGTCAACAGTGCGGTGAATGGCAGAGATTCGGTGAAGGCTTTACCGATAGCATGCTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGcc                                                                                                  >  5:242369/1‑147 (MQ=255)                                                gTCGATAATCACCGCGACTACCGAGAAAAACACCGTCAACAGTGCGGTGAATGGCAGAGATTCGGTGAAGGCTTTACCGAAAGCATGCTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCa                                                                                                 <  4:250983/148‑1 (MQ=255)                                                 tCGATAATCACCGCGACTACCGAGAAAAACACCGTCAACAGTGCGGTGAATGGCAGAGATTCGGTGAAGGCTTTACCGATAGCATGCTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCaa                                                                                                >  2:146031/1‑148 (MQ=255)                                                     tAATCACCGCGACTACCGAGAAAAACACCGTCAACAGTGCGGTGAATGGCAGAGATTCGGTGAAGGCTTTACCGATAGCATGCTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTc                                                                                            >  5:294558/1‑148 (MQ=255)                                                      aaTCACCGCGACTACCGAGAAAAACACCGTCAACAGTGCGGTGAATGGCAGAGATTCGGTGAAGGCTTTACCGATAGCATGCTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTcc                                                                                           <  4:296774/148‑1 (MQ=255)                                                       aTCACCGCGACTACCGAGAAAAACACCGTCAACAGTGCGGTGAATGGCAGAGATTCGGTGAAGGCTTTACCGATAGCATGCTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCg                                                                                          <  5:206695/148‑1 (MQ=255)                                                         cACCGCGACTACCGAGAAAAACACCGTCAACAGTGCGGTGAATGGCAGAGATTCGGTGAAGGCTTTACCGATAGCATGCTCATCGGTGACGCCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGc                                                                                         >  3:284368/1‑147 (MQ=255)                                                         cACCGCGACTACCGAGAAAAACACCGTCAACAGTGCGGTGAATGGCAGAGATTCGGTGAAGGCTTTACCGATAGCATGCTCATCGGTGACGCCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGc                                                                                         >  3:289396/1‑147 (MQ=255)                                                         cACCGCGACTACCGAGAAAAACACCGTCAACAGTGCGGTGAATGGCAGAGATTCGGTGAAGGCTTTACCGATAGCATGCTCATCGGTGACGCCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGc                                                                                         <  4:284368/147‑1 (MQ=255)                                                         cACCGCGACTACCGAGAAAAACACCGTCAACAGTGCGGTGAATGGCAGAGATTCGGTGAAGGCTTTACCGATAGCATGCTCATCGGTGACGCCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGc                                                                                         <  4:289396/147‑1 (MQ=255)                                                            cgcgACTACCGAGAAAAACACCGTCAACAGTGCGGTGAATGGCAGAGATTCGGTGAAGGCTTTACCGATAGCATGCTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGc                                                                                         >  7:414163/1‑144 (MQ=255)                                                            cgcgACTACCGAGAAAAACACCGTCAACAGTGCGGTGAATGGCAGAGATTCGGTGAAGGCTTTACCGATAGCATGCTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGa                                                                                     <  5:67905/148‑1 (MQ=255)                                                                 cTACCGAGAAAAACACCGTCAACAGTGCGGTGAATGGCAGAGATTCGGTGAAGGCTTTACCGATAGCATGCTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAAtg                                                                                                                      >  8:235155/1‑110 (MQ=255)                                                                 cTACCGAGAAAAACACCGTCAACAGTGCGGTGAATGGCAGAGATTCGGTGAAGGCTTTACCGATAGCATGCTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAAtg                                                                                                                      <  7:235155/110‑1 (MQ=255)                                                                      gagaAAAACACCGTCAACAGTGCGGTGAATGGCAGAGATTCGGTGAAGGCTTTACCGATAGCATGCTCATCGGTGCCACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCACCGc                                                                             >  6:22703/1‑146 (MQ=255)                                                                        gaAAAACACCGTCAACAGTGCGGTGAATGGCAGAGATTCGGTGAAGGCTTTACCGATAGCATGCTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCa                                                                                                 >  8:146265/1‑124 (MQ=255)                                                                        gaAAAACACCGTCAACAGTGCGGTGAATGGCAGAGATTCGGTGAAGGCTTTACCGATAGCATGCTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCa                                                                                                 <  7:146265/124‑1 (MQ=255)                                                                         aaaaaCACCGTCAACAGTGCGGTGAATGGCAGAGATTCGGTGAAGGCTTTACCGATAGCATGCTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAgcg                                                                        >  5:166987/1‑148 (MQ=255)                                                                          aaaaCACCGTCAACAGTGCGGTGAATGGCAGAGATTCGGTGAAGGCTTTACCGATAGCATGCTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAgcgc                                                                       <  5:310181/148‑1 (MQ=255)                                                                            aaCACCGTCAACAGTGCGGTGAATGGCAGAGATTCGGTGAAGGCTTTACCGATAGCATGCTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAg                                                                                      >  8:96543/1‑131 (MQ=255)                                                                            aaCACCGTCAACAGTGCGGTGAATGGCAGAGATTCGGTGAAGGCTTTACCGATAGCATGCTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAg                                                                                      <  7:96543/131‑1 (MQ=255)                                                                             acacCGTCAACAGTGCGGTGAATGGCAGAGATTCGGTGAAGGCTTTACCGATAGCATGCTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGt                                                                    <  3:234165/148‑1 (MQ=255)                                                                              cacCGTCAACAGTGCGGTGAATGGCAGAGATTCGGTGAAGGCTTTACCGATAGCATGCTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTc                                                                   <  8:414163/148‑1 (MQ=255)                                                                               acCGTCAACAGTGCGGTGAATGGCAGAGATTCGGTGAAGGCTTTACCGATAGCATGCTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCa                                                                  <  6:242369/148‑1 (MQ=255)                                                                                ccGTCAACAGTGCGGTGAATGGCAGAGATTCGGTGAAGGCTTTACCGATAGCATGCTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGCCTGCCAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCAc                                                                 >  7:209714/1‑148 (MQ=255)                                                                                ccGTCAACAGTGCGGTGAATGGCAGAGATTCGGTGAAGGCTTTACCGATAGCATGCTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCAc                                                                 <  1:146031/148‑1 (MQ=255)                                                                                    cAACAGTGCGGTGAATGGCAGAGATTCGGTGAAGGCTTTACCGATAGCATGCTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCAccagc                                                             >  3:126122/1‑148 (MQ=255)                                                                                    cAACAGTGCGGTGAATGGCAGAGATTCGGTGAAGGCTTTACCGATAGCATGCTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCAccagc                                                             >  6:22180/1‑148 (MQ=255)                                                                                               tGAATTGCAGAGATTCGGTGAAGGCTTTACCGATAGCATGCTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCat                                                   >  3:167909/1‑146 (MQ=255)                                                                                                  aTGGCAGAGATTCGGTGAAGGCTTTACCGATAGCATGCTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCaataat                                               <  8:305983/148‑1 (MQ=255)                                                                                                  aTGGCAGAGATTCGGTGAAGGCTTTACCGATAGCATGCTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGAAAAACCAATCAAGCCAAATTCCGCCAGATTCAACGCAAGCGCAGTCACCAGCCAGACGACaataat                                               <  8:300416/148‑1 (MQ=255)                                                                                                     gCAGAGATTCGGTGAAGGCTTTACCGATAGCATGCTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGc                                            <  8:963/148‑1 (MQ=255)                                                                                                     gCAGAGATTCGGTGAAGGCTTTACCGATAGCATGCTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGc                                            <  6:166987/148‑1 (MQ=255)                                                                                                      cAGAGATTCGGTGAAGGCTTTACCGATAGCATGCTCATCGGTGACACCGATCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGcc                                           <  6:291948/148‑1 (MQ=255)                                                                                                         agaTTCGGTGAAGGCTTTACCGATAGCATGCTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCAc                                                                 <  1:274733/123‑1 (MQ=255)                                                                                                         agaTTCGGTGAAGGCTTTACCGATAGCATGCTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCAc                                                                 >  2:274733/1‑123 (MQ=255)                                                                                                         agaTTCGGTGAAGGCTTTACCGATAGCATGCTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTgg                                        <  8:21937/148‑1 (MQ=255)                                                                                                            ttCGGTGAAGGCTTTACCGATAGCATGCTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGAc                                      >  4:37323/1‑147 (MQ=255)                                                                                                             tCGGTGAAGGCTTTACCGATAGCATGCTCATCGGTGACACCGGTCAATGATGTTGCCAGATTATTCCGTG‑CAAACCAATCCACCCAACGTCCGCCAGCTGCAAAGCCAGCGCAGTAAACCGCCAGCCCCCAATAATCCCCGGGACaa                                    >  7:195176/1‑147 (MQ=255)                                                                                                             tCGGTGAAGGCTTTACCGATAGCATGCTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGcc                                                     <  1:341251/131‑1 (MQ=255)                                                                                                             tCGGTGAAGGCTTTACCGATAGCATGCTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAACCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGcc                                                     >  2:341251/1‑131 (MQ=255)                                                                                                               ggTGAAGGCTTTACCGATAGCATGCTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACAAt                                   <  3:37323/147‑1 (MQ=255)                                                                                                                       cTTTACCGATAGCATGCTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACAATCAGACGGAt                          <  7:184336/148‑1 (MQ=255)                                                                                                                          tACCGATAGCATGCTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACAATCAGACGGAttt                        >  8:396322/1‑147 (MQ=255)                                                                                                                            ccGATAGCATGCTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGcc                                                                                        <  6:154764/81‑1 (MQ=255)                                                                                                                            ccGATAGCATGCTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGcc                                                                                        >  5:154764/1‑81 (MQ=255)                                                                                                                              gATAGCATGCTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACAATCAGACGGATTTTAt                     >  7:76817/1‑146 (MQ=255)                                                                                                                              gATAGCATGCTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACAATCAGACGGATTTTAt                     <  8:76817/146‑1 (MQ=255)                                                                                                                                tAGCATGCTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACAATCAGACGGATTTTATCCt                  >  6:233027/1‑147 (MQ=255)                                                                                                                                          aTCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACAATCAGACGGATTTTATCCTGACGGGTACg        <  4:126122/147‑1 (MQ=255)                                                                                                                                               tGACACCGGTCAATGATGTTGCCAGAATTATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACAATCAGACGGATTTTATCCTGACGGGTACGCtggtg  <  8:7594/148‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                TCCGATGCCTGCAACACAAACTGAATAATTGGCGAAAACAGTTGTTGGTCGATAATCACCGCGACTACCGAGAAAAACACCGTCAACAGTGCGGTGAATGGCAGAGATTCGGTGAAGGCTTTACCGATAGCATGCTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACAATCAGACGGATTTTATCCTGACGGGTACGCTGGTG  >  NC_000913/1233564‑1233853

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