Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,255,724	A→T	6.8%	Q69L (CAG→CTG)	folM →	bifunctional dihydrofolate reductase/dihydromonapterin reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,255,724	0	A	T	6.8%	173.9 / 4.5	60	Q69L (CAG→CTG)	folM	bifunctional dihydrofolate reductase/dihydromonapterin reductase
Reads supporting (aligned to +/- strand):  ref base A (25/30);  new base T (2/2);  total (27/33)
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.90e-01

CATTCTCGTTACCGGAGCCAGCCAGCGCGTCGGGCTGGCCCTGGCCCTCGAACTGGCACAGGCCGGCCATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTCGACAACCTCGAGACCATGTTCCGCCTGCACGTCGAAGCGCCCTACCAC  >  NC_002947/5255575‑5255872                                                                                                                                                       |                                                                                                                                                     cATTCTCGTTACCGGAGCCAGCCAGCGCGTCGGGCTGGCCCTGGCCCTCGAACTGGCACAGGCCGGCCATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCa                                                                                                                                                      >  1:1034814/1‑150 (MQ=255) cATTCTCGTTACCGGAGCCAGCCAGCGCGTCGGGCTGGCCCTGGCCCTCGAACTGGCACAGGCCGGCCATACGGGGGTCAGTGCCAGCAGCCGCATCATGTTGCAGGCTGCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCa                                                                                                                                                      >  1:1174466/1‑150 (MQ=255)    tctcGTTACCGGAGCCAGCCAGCGCGTCGGGCTGGCCCTGGCCCTCGAACTGGCACAGGCCGGCCATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAggc                                                                                                                                                   <  2:1100183/150‑1 (MQ=255)     ctcGTTACCGGAGCCAGCCAGCGCGTCGGGCTGGCCCTGGCCCTCGAACTGGCACAGGCCGGCCATACGGTGGTCAGTGCCAGCCGCAGCAT‑CAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAggcgg                                                                                                                                                 <  2:573367/150‑1 (MQ=255)      tcGTTACCGGAGCCAGCCAGCGCGTCGGGCTGGCCCTGGCCCTCGAACTGGCACAGGCCGGCCATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAggcgg                                                                                                                                                 <  1:807564/150‑1 (MQ=255)        gTTACCGGAGCCAGCCAGCGCGTCGGGCTGGCCCTGGCCCTCGAACTGGCACAGGCCGGCCATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAggcgggc                                                                                                                                               >  2:1101190/1‑150 (MQ=255)              ggAGCCAGCCAGCGCGTCGGGCTGGCCCTGGCCCTCGAACTGGCACAGGCCGGCCATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAGGCGGGCGACCGc                                                                                                                                         >  2:509476/1‑150 (MQ=255)              ggAGCCAGCCAGCGCGTCGGGCTGGCCCTGGCCCTCGAACTGGCACAGGCCGGCCATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAGGCGGGCGACCGc                                                                                                                                         <  1:53160/150‑1 (MQ=255)                agccagccAGCGCGTCGGGCTGGCCCTGGCCCTCGAACTGGCACAGGCCGGCCATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAGGCGGGCGACCGCCa                                                                                                                                       <  1:718372/150‑1 (MQ=255)                 gccagccaGCGCGTCGGGCTGGCCCTGGCCCTCGAACTGGCACAGGCCGGCCATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCTGGCGGTcgcccgcctg                                                                                                                                      >  2:488084/1‑148 (MQ=255)                 gccagccaGCGCGTCGGGCTGGCCCTGGCCCTCGAACTGGCACAAGCCGGCCATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAGGCGGGCGACCGCCAg                                                                                                                                      <  2:709288/150‑1 (MQ=255)                    agccagcGCGTCGGGCTGGCCCTGGCCCTCGAACTGGCACAGGCCGGCCATACGGTGGTCAGTGCCAACCGCAACATCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAGGCGGGCGACCGCCAGgcg                                                                                                                                   >  1:319970/1‑150 (MQ=255)                       cagcGCGTCGGGCTGGCCCTGGCCCTCGAACTGGCACAGGCCGGCCATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCGG‑GCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTg                                                                                                                                >  1:119789/1‑150 (MQ=255)                       cagcGCGTCGGGCTGGCCCTGGCCCTCGAACTGGCACAGGCCGGCCATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCTGGCGGTCGCCCGCCTGgcacag                                                                                                                                <  1:488084/150‑5 (MQ=255)                        agcGCGTCGGGCTGGCCCTGGCCCTCGAACTGGCACAGGCCGGCCATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTga                                                                                                                               <  2:59406/150‑1 (MQ=255)                         gcgcGTCGGGCTGGCCCTGGCCCTCGAACTGGCACAGGCCGGCCATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTgat                                                                                                                              <  1:954818/150‑1 (MQ=255)                          cgcgTCGGGCTGGCCCTGGCCCTCGAACTGGCACAGGCCGGCCATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGCCCTCGGCGGGCGACCGCCAGGCGCTgatc                                                                                                                             <  1:1069034/150‑1 (MQ=255)                                gggCTGGCCCTGGCCCTCGAACTGGCACAGGCCGGCCATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCAc                                                                                                                       <  1:496429/150‑1 (MQ=255)                                  gCTGGCCCTGGCCCTCGAACTGGCACAGGCCGGCCATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAggcgg                                                                                                                                                 >  1:637877/1‑122 (MQ=255)                                  gCTGGCCCTGGCCCTCGAACTGGCACAGGCCGGCCATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAggcgg                                                                                                                                                 <  2:637877/122‑1 (MQ=255)                                      gccctggccctCGAACTGGCACAGGCCGGCCATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCCACATCGTGCAGTTTCAG‑GACCACCTGTGCCAGGCGCGCGACCGCCAGGCGCTGATCGATCACCTGCAc                                                                                                                 <  2:1043327/150‑1 (MQ=255)                                          tggccctCGAACTGGCACAGGCCGGCCATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATc                                                                                                             <  1:703208/150‑1 (MQ=255)                                             ccctCGAACTGGCACAGGCCGGCCATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAG‑GCCGAACTGTGCCAGGTGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGtt                                                                                                          <  1:788359/150‑1 (MQ=255)                                              cctCGAACTGGCACAGGCCGGCCATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAGACGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTa                                                                                                         >  1:1116888/1‑150 (MQ=255)                                               ctCGAACTGGCACAGGCCGGCCATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTAt                                                                                                        >  2:147657/1‑150 (MQ=255)                                                tCGAACTGGCACAGGCCGGCCATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATg                                                                                                       >  2:700852/1‑150 (MQ=255)                                                          acaGGCCGGCCATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTAcg                                                                                             >  1:1000228/1‑150 (MQ=255)                                                          acaGGCCGGCCATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTAcg                                                                                             >  2:1172947/1‑150 (MQ=255)                                                                cggcgATACGGTGGTCAATGGCAGTGGCAGCAGCCAGCCGCAGGCAGCCCACCCCCACATCGTGCA‑TACGAGCTTCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCcat                                                                                       <  1:677421/150‑1 (MQ=255)                                                                               cAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTc                                                                        <  2:1000228/150‑1 (MQ=255)                                                                                  tGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCt                                                                     >  1:667898/1‑150 (MQ=255)                                                                                    ccAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAGGCTTGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTAt                                                                   <  1:509476/150‑1 (MQ=255)                                                                                    ccAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTAt                                                                   >  2:110928/1‑150 (MQ=255)                                                                                      aGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATgg                                                                 <  2:1174466/150‑1 (MQ=255)                                                                                             gcaTCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTcgacg                                                          >  2:934757/1‑150 (MQ=255)                                                                                                tCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGacgacg                                                       <  2:1034814/150‑1 (MQ=255)                                                                                                   aGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCTGGCGGTcgcccgcctgg                                                                                                                                     >  1:537485/1‑66 (MQ=37)                                                                                                   aGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCTGGCGGTCGCCCGCCTgg                                                                                                                                     <  2:537485/69‑1 (MQ=37)                                                                                                   aGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCCTCCACAACGCCTCGCTATGGCTCGACGACGGcc                                                    >  2:550096/1‑150 (MQ=255)                                                                                                   aGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACCACGCCTCGCTATGGCTCGACGACGGcc                                                    <  2:625613/150‑1 (MQ=255)                                                                                                   aGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGcc                                                    >  1:402922/1‑150 (MQ=255)                                                                                                     ctgcaggcagCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTc                                                  >  1:949188/3‑150 (MQ=255)                                                                                                      cGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCCAGAGTAGATACGATGGGTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTCg                                                 >  1:1114062/1‑150 (MQ=255)                                                                                                             agCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTCGACAACCt                                          >  2:930100/1‑150 (MQ=255)                                                                                                              gCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTCGACAACCTc                                         >  2:224176/1‑150 (MQ=255)                                                                                                               cccACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATg                                                                                                       >  2:217657/1‑87 (MQ=255)                                                                                                               cccACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATg                                                                                                       <  1:217657/87‑1 (MQ=255)                                                                                                                ccACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTCGACAACCTCga                                       <  2:1027003/150‑1 (MQ=255)                                                                                                                   tcccAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAGGCGCGCTACCCCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTCGACAACCTCGAGAc                                    <  2:460651/149‑1 (MQ=255)                                                                                                                    cccAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTCGACAACCTCGAGAcc                                   >  1:78654/1‑150 (MQ=255)                                                                                                                     ccAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTCGACAACCTCGAGACCa                                  >  1:460651/1‑150 (MQ=255)                                                                                                                      cAACATCGTGCAGTTTCAG‑GCCGACCTGTTGCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTCGACAACCTCGAGACCAt                                 >  1:1125507/1‑150 (MQ=255)                                                                                                                      cAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTCGACAACCTCGAGACCAt                                 >  1:461739/1‑150 (MQ=255)                                                                                                                             gTGCAGTTTCAG‑GCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGCCGGACTACGCGCCATTCACGACAACGCCTCGCTATGGCTCGACGACGGCCTCGACAACCTCGAGACCATGTTCCGc                          <  1:147657/150‑1 (MQ=255)                                                                                                                             gTGCAGTTTCAG‑GCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTCGACAACCTCGAGACCATGTTCCGc                          >  1:225754/1‑150 (MQ=255)                                                                                                                             gTGCAGTTTCAG‑GCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTCGACAACCTCGAGACCATGTTCCGc                          <  1:14414/150‑1 (MQ=255)                                                                                                                                          gCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTCGACAACCGCGAGACCATGTTCCGCCTGCACGTCGaa              <  1:550096/150‑1 (MQ=255)                                                                                                                                            cGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTCGACAACCTCGAGACCATGTTCCGCCTGCACGTCGAAgc            <  1:700852/150‑1 (MQ=255)                                                                                                                                              aCCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTCGACAACCTCGAGACCATGTTCCGCCTGCACGTCGAAgcgc          <  2:402922/150‑1 (MQ=255)                                                                                                                                              aCCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTCGACAACCTCGAGACCATGTTCCGCCTGCACGTCGAAgcgc          <  1:930100/150‑1 (MQ=255)                                                                                                                                                    gCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTCGACAACCTCGAGACCATGTTCCGCCTGCACGTCGAAGCGCCCTacc    <  2:78654/150‑1 (MQ=255)                                                                                                                                                      cAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTCGACAACCTCGAGACCATGTTCCGCCTGCACGTCGAAGCGCCCTaccac  <  2:225754/150‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                     CATTCTCGTTACCGGAGCCAGCCAGCGCGTCGGGCTGGCCCTGGCCCTCGAACTGGCACAGGCCGGCCATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTCGACAACCTCGAGACCATGTTCCGCCTGCACGTCGAAGCGCCCTACCAC  >  NC_002947/5255575‑5255872

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