Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,835,659	G→A	6.5%	G49S (GGC→AGC)	ynjA →	carboxymuconolactone decarboxylase family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,835,659	0	G	A	6.5%	192.7 / 2.2	61	G49S (GGC→AGC)	ynjA	carboxymuconolactone decarboxylase family protein
Reads supporting (aligned to +/- strand):  ref base G (29/28);  new base A (2/2);  total (31/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 = 1.00e+00

TGATGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCGCTTGTCAGCGCGCGCATTGCGCAAATGTGCCTGTGTGAGTTTTGTGTGGATATCACCAGTATGAAAGTCGCCGAGCGCACCGGCAGCAGCGAT  >  NC_000913/1835513‑1835799                                                                                                                                                   |                                                                                                                                             tGATGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTtggc                                                                                                                                            >  2:60980/1‑149 (MQ=255)       gTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCTTTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTGGATGTTTGTTAgctggctg                                                                                                                                      >  7:257145/1‑149 (MQ=255)       gTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTAgctggctg                                                                                                                                      >  7:263083/1‑149 (MQ=255)       gTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGCCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAATCCCCATTCGCTGGTGGGTGCGGATCCCGTTTAGCTTTTATCTTGTGTCGATGTTTGTTGgctggctg                                                                                                                                      >  6:317081/1‑149 (MQ=255)         ttACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGgctggctg                                                                                                                                      >  3:14953/1‑147 (MQ=255)          tACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGa                                                                                                                                    <  8:26100/148‑1 (MQ=255)                cgcTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCaaa                                                                                                                             <  4:14953/149‑1 (MQ=255)                 gcTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGAGCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAAc                                                                                                                            >  4:149006/1‑149 (MQ=255)                                tttAATTTTACTTCCACAGGCGTGGATGCATCGTCGAAATTAAGGCTAGGTGATAAGCCCCATTCTCTTGTGGGGGCGGATCCAGTTTATCTTTTATCTGGTGTCGATGTTTGTTGTCTGGCTGGAGAGAAAACTCTCAACGCTAGATc                                                                                                             >  6:270292/1‑149 (MQ=255)                                tttAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATc                                                                                                             >  1:152274/1‑149 (MQ=255)                                tttAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATc                                                                                                             >  3:81304/1‑149 (MQ=255)                                tttAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATc                                                                                                             >  2:244620/1‑149 (MQ=255)                                tttAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCAATTCGCTTGTGCGGGCGGATCCCGTTTATCGTTTATCTGGTCTCGATGTTTGTTGGCTGGCTGGAGCTCAAACGCTCACAGCTAGATc                                                                                                             >  5:142348/1‑149 (MQ=255)                                    aTTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCgg                                                                                                          <  1:60980/148‑1 (MQ=255)                                    aTTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCgg                                                                                                          >  6:370341/1‑148 (MQ=255)                                     ttttACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCgg                                                                                                          >  1:341470/1‑147 (MQ=255)                                              cacaGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATc                                                                                                <  2:152274/148‑1 (MQ=255)                                               acaGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCGc                                                                                              >  1:38508/1‑149 (MQ=255)                                               acaGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCGc                                                                                              <  5:370341/149‑1 (MQ=255)                                                          tGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCGCTTGTCa                                                                                        >  8:30054/1‑144 (MQ=255)                                                          tGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCGCTTGTCa                                                                                        <  7:30054/144‑1 (MQ=255)                                                                gtcgCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCGTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCGCTTGTCAGCGCGCGCAtt                                                                             >  2:344470/1‑149 (MQ=255)                                                                 tcgCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTTGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCGCTTGTCAGCGCGCGCATTg                                                                            <  8:327991/149‑1 (MQ=255)                                                                 tcgCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCGCTTGTCAGCGCGCGCATTg                                                                            <  2:341470/149‑1 (MQ=255)                                                                   gCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGGTTATCTTTTATCTGGGGGCGATGTTTGTTGGCGGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCGCTTGTCAGCGCGCGCATTgc                                                                           <  4:32464/148‑1 (MQ=255)                                                                       ttACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGGCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCGCTTGTCAGCGCGCGCATTGCGCaaa                                                                      <  2:196264/149‑1 (MQ=255)                                                                          cGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCGCTTGTCAGCGCGCGCATTGCGCAAAtgt                                                                   <  2:38508/149‑1 (MQ=255)                                                                            gCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCGCTTGTCAGCGCGCGCATTGCGCAAATGTGc                                                                 >  4:149053/1‑149 (MQ=255)                                                                              gAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCGCTTGTCAGCGCGCGCATTGCGCAAATGTGCCt                                                               >  3:121567/1‑149 (MQ=255)                                                                                    cTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCGCTTGTCAGCGCGCGCATTGCGCAAATGTGCCtgtgt                                                           >  1:176157/1‑147 (MQ=255)                                                                                    cTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCGCTTGTCAGCGCGCGCATTGCGCAAATGTGCCTGTGTGa                                                         <  5:165287/149‑1 (MQ=255)                                                                                       aGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCGGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCGCTTGTCAGCGCGCGCAGTGCGCAAATGTGCCTGTGTGAg                                                        <  5:146962/147‑1 (MQ=255)                                                                                        gCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCGCtt                                                                                            >  6:303592/1‑110 (MQ=255)                                                                                         ccccATGCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGTTGTCGATGGTTGTTGGCTGGCTGGAGCGCAAACGCTGACCGCGCGATCCGGTAGTAAGATCGCTTGTCAGCGCGCGCATTGCGCAAATGTGCCTGTGGGAGgttt                                                    <  8:151501/149‑1 (MQ=255)                                                                                          cccTTTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCTATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTTGTACGATCGCTTGTCAGCGCGCGCATTGCGCAAATGTGCCTGTGTGAGtttt                                                    <  6:142348/148‑1 (MQ=255)                                                                                           ccATTCGCTGGTGGGGGCGGATCCCGGGTATCTTTTATCTGGGGTCGATGTGTGTTGGCTGGCTGGAGCGCAAACTCTCACCGCTCGATCCGGTAGTACGATCGCAGGTCAGCGCGCGCATTGAGCAAATGTGCCTGTGTGAGtt                                                      >  3:318671/1‑145 (MQ=255)                                                                                           cacttCTCTGTTGGTGGCGGATCCCGTTTATCTTTTCTCTGGTGTCGATGTTTGTTGGCTGGCTTGAGCGCAAACGCTCACCGCTCGATCCGGTAGTCCGATCGCTTGTCAGCGCGCGCATTGCGCAAATGTGCCTGTGTGAGtt                                                      <  4:318671/142‑1 (MQ=255)                                                                                             aTTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCGCTTGTCAGCGCGCGCATTGCGCAAATGTGCCTGTGTGAGTTTtgtgt                                                <  8:173638/149‑1 (MQ=255)                                                                                              ttCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCGCTTGTCAGCGCGCGCATTGCGCAAATGTGCCTGTGTGAGTTTtgtgtg                                               <  7:393404/149‑1 (MQ=255)                                                                                                cGCTGGGGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGGAGTACGATCGCTTGTCAGCGCGCGCATTGCGCAAATGTGCCTGTGTGAGTTTTGTGTgg                                              <  4:308522/148‑1 (MQ=255)                                                                                                   tggtggGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCGCTTGTCAGCGCGCGCATTGCGCAAATGTGCCTGTGTGAGTTTTGTGTGGatat                                          >  5:18074/1‑149 (MQ=255)                                                                                                       gggggCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTAGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCGCTTGTCAGCGCGCGCATTGCGCAAATGTGCCTGTGTGAGTTTTGTGTGGATATCAcc                                      <  8:263083/149‑1 (MQ=255)                                                                                                        ggggCTGATCCCGATTATCTTTAATCTGGTGTCGATGTTTGTTAGCTGGCTGGACCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCGCTTGTCAGCGCGCGCATTGCGCAAATGTGCCTGTGTGAGTTTTGTGTGGATATCAcc                                      <  8:257145/148‑1 (MQ=255)                                                                                                        ggggCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCGCTTGTCAGCGCGCGCATTGCGCAAATGTGCCTGTGTGAGTTTTGTGTGGATATCAc                                       >  6:344332/1‑147 (MQ=255)                                                                                                            cGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCGCTTGTCAGCGCGCGCATTGCGCAAATGTGCCTGTGTGAGTTTTGTGTGGATATCACCAGTAt                                 <  4:121567/149‑1 (MQ=255)                                                                                                                 cccGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCGCTTGTCAGCGCGCGCATTGCGCAAATGTGCCTGTGTGAGTTTtg                                                   <  7:79830/126‑1 (MQ=255)                                                                                                                   cGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCGCTTGTCAGCGCGCGCATTGCGCAAATGTGCCTGTGTGAGTTTTGTGTGGATATCACCAGTATGAAAGTc                          >  7:353530/1‑149 (MQ=255)                                                                                                                    gTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCGCTTGTCAGCGCGCGCATTGCGCAAATGTGCCTGTGTGAGtttt                                                    >  6:169143/1‑122 (MQ=255)                                                                                                                    gTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCGCTTGTCAGCGCGCGCATTGCGCAAATGTGCCTGTGTGAGtttt                                                    <  5:169143/122‑1 (MQ=255)                                                                                                                    gTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCGCTTGTCAGCGCGCGCATTGCGCAAATGTGCCTGTGTGAGTTTTGTGTGGATATCACCAGTATGAAAGTCg                         >  4:279882/1‑149 (MQ=255)                                                                                                                       tatcttttatctGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCGCTTGTCAGCGCGCGCATTGCGCAAATGTGCCTGTGTGAGTTTTGTGTGGATATCACCAGTATGa                               <  1:171705/140‑1 (MQ=255)                                                                                                                       tatcttttatctGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCGCTTGTCAGCGCGCGCATTGCGCAAATGTGCCTGTGTGAGTTTTGTGTGGATATCACCAGTATGa                               >  2:171705/1‑140 (MQ=255)                                                                                                                        atcttttatctGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCGCTTGTCAGCGCGCGCATTGCGCAAATGTGCCTGTGTGAGTTTTGTGTGGATATCACCAGTATGAAAGTCGcc                       >  2:220113/1‑147 (MQ=255)                                                                                                                           ttttatctGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCGCTTGTCAGCGCGCGCATTGCGCAAATGTGCCTGTGTGAGTTTTGTGTGGATATCACCAGTATGAAAGTCGCCGAgcg                  >  1:2392/1‑149 (MQ=255)                                                                                                                            tttatctGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCGCTTGTCAGCGCGCGCATTGCGCAAATGTGCCTGTGTGAGTTTTGTGTGGATATCACCAg                                    >  2:221300/1‑130 (MQ=255)                                                                                                                            tttatctGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCGCTTGTCAGCGCGCGCATTGCGCAAATGTGCCTGTGTGAGTTTTGTGTGGATATCACCAg                                    <  1:221300/130‑1 (MQ=255)                                                                                                                             ttatctGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCGCTTGTCAGCGCGCGCATTGCGCAAATGTGCCTGTGTGAGTTTTGTGTGGATATCACCAGTATGAAAGTCGCCGAgc                   <  5:344332/146‑1 (MQ=255)                                                                                                                                  tGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCGCTTGTCAGCGCGCGCATTGCGCAAATGTGCCTGTGTGAGTTTTGTGTGGATATCACCAGTATGAAAGTCGCCGAGCGCACCgg            >  3:152220/1‑148 (MQ=255)                                                                                                                                    gtgtCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCGCTTGTCAGCGCGCGCATTGCGCAAATGTGCCTGTGTGAGTTTTGTGTGGATATCACCAGTATGAAAGTCGCCGAGCGCACCgg            <  4:81304/146‑1 (MQ=255)                                                                                                                                           tttttgttGGATGGCTGGAGAGCAAAATCTAACAGCTCGATCCGGTAGTACTATAGCTTTTCAGCGCGCGCAATGCGCAAATGTGCCTGTGTGAGTTTTGTGTTGATATCACCAGTATGAAAGTCGCCGAGCGCACCGGCAGCAGCGAt  <  5:270292/147‑1 (MQ=255)                                                                                                                                            gtttgttGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCGCTTGTCAGCGCGCGCATTGCGCAAATGTGCCTGTGTGAGTTTTGTGTGGATATCa                                        <  7:318662/110‑1 (MQ=255)                                                                                                                                            gtttgttGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCGCTTGTCAGCGCGCGCATTGCGCAAATGTGCCTGTGTGAGTTTTGTGTGGATATCa                                        >  8:318662/1‑110 (MQ=255)                                                                                                                                                   |                                                                                                                                             TGATGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCGCTTGTCAGCGCGCGCATTGCGCAAATGTGCCTGTGTGAGTTTTGTGTGGATATCACCAGTATGAAAGTCGCCGAGCGCACCGGCAGCAGCGAT  >  NC_000913/1835513‑1835799

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