Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,107,295	G→A	6.7%	R289C (CGC→TGC)	ydiJ ←	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,107,295	0	G	A	6.7%	172.9 / 2.8	60	R289C (CGC→TGC)	ydiJ	putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (31/25);  new base A (2/2);  total (33/27)
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.10e-01

TCCGCCGGGAAGTACTCGGCAACGCTGTGCCAGACGATGTCTTTCATCGCCAGCATCAGCACCTTGGAGTCGACCGTCTCGATCGACAGCGGCTTGTGCGCCATCAGCGCGTTGGCATCACGCAGCGCATCCATGAAGCTGGTGTAGCGCACATTGACCAGCACCGCATACTTGGGAATCGGCAGCACGTTGAGCTTGGCTTCGACCACGTAG‑‑CCCAGCGA‑GCCCTCGGCGCCGCACAGTACGCTGTTGAGGTTGAAGCGGCCCTGCTCGTCGCGCAGGTGCGCCAGGTCGTAGCCGG  >  NC_002947/5107147‑5107444                                                                                                                                                     |                                                                                                                                                         tCCGCCGGGAAGTACTCGGCAACGCTGTGCCAGACGATGTCTTTCATCGCCAGCATCAGCACCTTGGAGTCGACCGTCTCGATCGACAGCGGCTTGTGCGCCATCAGCGCGTTGGCATCACGCAGCGCATCCATGAAGCTGGTGTAgcgc                                                                                                                                                         <  1:1220859/150‑1 (MQ=255)  ccgccgGGAAGTACTCGGCAACGCTGTGCCAGACGATGTCTTTCATCGCCAGCATCAGCACCTTGGAGTCGACCGTCTCGATCGACAGCGGCTTGTGCGCCATCAGCGCGTTGGCATCACGCAGCGCATCCATGAAGCTGGTGTAGCGca                                                                                                                                                        <  2:1631438/150‑1 (MQ=255)   cgccgGGAAGTACTCGGCAACGCTGTGCCAGACGATGTCTTTCATCGCCAGCATCAGCACCTTGGAGTCGACCGTCTCGATCGACAGCGGCTTGTGCGCCATCAGCGCGTTGGCATCACGCAGCGCATCCATGAAGCTGGTGTAGCGcac                                                                                                                                                       >  1:163001/1‑150 (MQ=255)        ggAGGTGCTCAGCAACGCTTTGCCAGACGAGGTGGTTCATCGCCAGCATCAGCACCTTGGAGTCGACCGTCTCGATCGACAGCGGCTTGTGCGCCATCAGCGCGTTGGCATCACGCAGCGCATCCATGAAGCTGGTGTAGCGCACATTGa                                                                                                                                                  >  2:518195/1‑150 (MQ=255)         gAAGTACTCGGCAACGCTGTGCCAGACGATGTCTTTCATCGCCAGCATCAGCACCTTGGAGTCGACCGTCTCGATCGACAGCGGCTTGTGCGCCATCAGCGCGTTGGCATCACGCAGCGCATCCATGAAGCTGGTGTAGCGCACATTGAc                                                                                                                                                 <  2:1285055/150‑1 (MQ=255)           aGTACTCGGCAACGCTGTGCCAGACGATGTCTTTCATCGCCAGCATCAGCACCTTGGAGTCGACCGTCTCGATCGACAGCGGCTTGTGCGCCATCAGCGCGTTGGCATCACGCAGCGCATCCATGAAGCTGGTGTAGCGCACATTGACCa                                                                                                                                               <  2:1682646/150‑1 (MQ=255)             tACTCGGCAACGCTGTGCCAGACGATGTCTTTCATCGCCAGCATCAGCACCTTGGAGTCGACCGTCTCGATCGACAGCGGCTTGTGCGCCATCAGCGCGTTGGCATCACGCAGCGCATCCATGAAGCTGGTGTAGCGCACATTGa                                                                                                                                                  <  1:1127386/145‑1 (MQ=255)             tACTCGGCAACGCTGTGCCAGACGATGTCTTTCATCGCCAGCATCAGCACCTTGGAGTCGACCGTCTCGATCGACAGCGGCTTGTGCGCCATCAGCGCGTTGGCATCACGCAGCGCATCCATGAAGCTGGTGTAGCGCACATTGa                                                                                                                                                  >  2:1127386/1‑145 (MQ=255)             tACTCGGCAACGCTGTGCCAGACGATGTCTTTCATCGCCAGCATCAGCACCTTGGAGTCGACCGTCTCGATCGACAGCGGCTTGTGCGCCATCAGCGCGTTGGCATCACGCAGCGCATCCATGAAGCTGGTGTAGCGCACATTGACCAGc                                                                                                                                             >  1:843977/1‑150 (MQ=255)                tCGGCAACGCTGTGCCAGACGATGTCTTTCATCGCCAGCATCAGCACCTTGGAGTCGACCGTCTCGATCGACAGCGGCTTGTGCGCCATCAGCGCGTTGGCATCACGCAGCGCATCCATGAAGCTGGTGTAGCGCACATTGACCAGCAcc                                                                                                                                          >  2:976326/1‑150 (MQ=255)                         cTGTGCCAGACGATGTCTTTCATCGCCAGCATCAGCACCTTGGAGTCGACCGTCTCGATCGACAGCGGCTTGTGCGCCATCAGCGCGTTGGCATCACGCAGCGCATCCATGAAGCTGGTGTAGCGCACATTGACCAGCACCGCATACTTg                                                                                                                                 <  1:369763/150‑1 (MQ=255)                                  aCGATGTCTTTCATCGCCAGCATCAGCACCTTGGAGTCGACCGTCTCGATCGACAGCGGCTTGTGCGCCATCAGCGCGTTGGCATCACGCAGCGCATCCATGAAGCTGGTGTAGCGCACATTGACCAGCACCGCATACTTGGGAATCCgc                                                                                                                        <  1:1935208/150‑1 (MQ=255)                                        tCTTTCATCGCCAGCATCAGCACCTTGGAGTCGACCGTCTCGATCGACAGCGGCTTGTGCGCCATCAGCGCGTTGGCATCACGCAGCGCATCCATGAAGCTGGTGTAGCGCACATTGACCAGCACCGCATACTTgg                                                                                                                                >  2:1091358/1‑136 (MQ=255)                                        tCTTTCATCGCCAGCATCAGCACCTTGGAGTCGACCGTCTCGATCGACAGCGGCTTGTGCGCCATCAGCGCGTTGGCATCACGCAGCGCATCCATGAAGCTGGTGTAGCGCACATTGACCAGCACCGCATACTTgg                                                                                                                                <  1:1091358/136‑1 (MQ=255)                                          tttCATCGCCAGCATCAGCACCTTGGAGTCGACCGTCTCGATCGACAGCGGCTTGTGCGCCATCAGCGCGTTGGCATCACGCAGCGCATCCATGAAGCTGGTGTAGCGCACATTGACCAGCACCGCATACTTGGGAATCGGCAGCACGtt                                                                                                                >  2:576977/1‑150 (MQ=255)                                          tttCATCGCCAGCATCAGCACCTTGGAGTCGACCGTCTCGATCGACAGCGGCTTGTGCGCCATCAGCGCGTTGGCATCACGCAGCGCATCCATGAAGCTGGTGTAGCGCACATTGACCAGCACCGCATACTTGGGAATCGGCAGCACGtt                                                                                                                >  2:1216044/1‑150 (MQ=255)                                            tCATCGCCAGCATCAGCACCTTGGAGTCGACCGTCTCGATCGACAGCGGCTTGTGCGCCATCAGCGCGTTGGCATCACGCAGCGCATCCATGAAGCTGGTGTAGCGCACATTGACCAGCACCGCATACTTGGGAATCGGCAGCACGTTGa                                                                                                              <  1:576977/150‑1 (MQ=255)                                             cATCGCCAGCATCAGCACCTTGGAGTCGACCGTCTCGATCGACAGCGGCTTGTGCGCCATCAGCGCGTTGGCATCACGCAGCGCATCCATGAAGCTGGTGTAGCACACATTGACCAGCACCGCATACTTGGGAATCGGCAGCACGtt                                                                                                                >  2:1783781/1‑147 (MQ=255)                                             cATCGCCAGCATCAGCACCTTGGAGTCGACCGTCTCGATCGACAGCGGCTTGTGCGCCATCAGCGCGTTGGCATCACGCAGCGCATCCATGAAGCTGGTGTAGCACACATTGACCAGCACCGCATACTTGGGAATCGGCAGCACGtt                                                                                                                >  2:1784147/1‑147 (MQ=255)                                             cATCGCCAGCATCAGCACCTTGGAGTCGACCGTCTCGATCGACAGCGGCTTGTGCGCCATCAGCGCGTTGGCATCACGCAGCGCATCCATGAAGCTGGTGTAGCACACATTGACCAGCACCGCATACTTGGGAATCGGCAGCACGtt                                                                                                                <  1:1783781/147‑1 (MQ=255)                                             cATCGCCAGCATCAGCACCTTGGAGTCGACCGTCTCGATCGACAGCGGCTTGTGCGCCATCAGCGCGTTGGCATCACGCAGCGCATCCATGAAGCTGGTGTAGCACACATTGACCAGCACCGCATACTTGGGAATCGGCAGCACGtt                                                                                                                <  1:1784147/147‑1 (MQ=255)                                              aTCGCCAGCATCAGCACCTTGGAGTCGACCGTCTCGATCGACAGCGGCTTGTGCGCCATCAGCGCGTTGGCATCACGCAGCGCATCCATGAAGCTGGTGTAGCGCACATTGACCAGCACCGCATACTTGGGAATCGGCAGCACGTTGAGc                                                                                                            >  2:1635745/1‑150 (MQ=255)                                                cGCCAGCATCAGCACCTTGGAGTCGACCGTCTCGATCGACAGCGGCTTGTGCGCCATCAGCGCGTTGGCATCACGCAGCGCATCCATGAAGCTGGTGTAGCGCACATTGACCAGCACCGCATACTTGGGAATCGGCAGCACGTTGAGCtt                                                                                                          <  1:518195/150‑1 (MQ=255)                                                  ccAGCATCAGCACTTTGGAGTCGACCGTCTCGATCGACAGCGGCTTGTGCGCCATCAGCGCGTTGGCATCACGCAGCGCATCCATGAAGCTGGTGTAGCGCACATTGACCAGCACCGCATACTTGGGAATCGGCAGCACGTTGAGCTTgg                                                                                                        <  2:936410/150‑1 (MQ=255)                                                  ccAGCATCAGCACCTTGGAGTCGACCGTCTCGATCGACAGCGGCTTGTGCGCCATCAGCGCGTTGGCATCACGCAGCGCATCCATGAAGCTGGTGTAGCGCACATTGACCAGCACCGCATACTTGGGAATc                                                                                                                           >  1:399967/1‑131 (MQ=255)                                                  ccAGCATCAGCACCTTGGAGTCGACCGTCTCGATCGACAGCGGCTTGTGCGCCATCAGCGCGTTGGCATCACGCAGCGCATCCATGAAGCTGGTGTAGCGCACATTGACCAGCACCGCATACTTGGGAATc                                                                                                                           <  2:399967/131‑1 (MQ=255)                                                  ccAGCATCAGCACCTTGGAGTCGACCGTCTCGATCGACAGCGGCTTGTGCGCCATCAGCGCGTTGGCATCACGCAGCGCATCCATGAAGCTGGTGTAGCGCACATTGACCAGCACCGCATACTTGGGAATCGGCAGCACGTTGAGCTTgg                                                                                                        >  2:1899350/1‑150 (MQ=255)                                                      cATCAGCACCTTGGAGTCGACCGTCTCGATCGACAGCGGCTTGTGCGCCATCAGCGCGTTGGCATCACGCAGCGCATCCATGAAGCTGGTGTAGCGCACATTGACCAGCACCGCATACTTGGGAATCGGCAGCACGTTGAGCTTGGCTTc                                                                                                    <  2:1639010/150‑1 (MQ=255)                                                        tCAGCACCTTGGAGTCGACCGTCTCGATCGACAGCGGCTTGTGCGCCATCAGCGCGTTGGCATCACGCAGCGCATCCATGAAGCTGGTGTAGCGCACATTGACCAGCACCGCATACTTGGGAATCGGCAGCACGTTGAGCTTGGCt                                                                                                      >  2:66702/1‑146 (MQ=255)                                                        tCAGCACCTTGGAGTCGACCGTCTCGATCGACAGCGGCTTGTGCGCCATCAGCGCGTTGGCATCACGCAGCGCATCCATGAAGCTGGTGTAGCGCACATTGACCAGCACCGCATACTTGGGAATCGGCAGCACGTTGAGCTTGGCt                                                                                                      <  1:66702/146‑1 (MQ=255)                                                           gCACCTTGGAGTCGACCGTCTCGATCGACAGCGGCTTGTGCGCCATCAGCGCGTTGGTATCACGCAGCTCATCCATGTTGCTGGTGTAGCGAACATTTGGCAGCACCGCATGCGTGGGAATCGGCAGCACGTTGAGCTTGGCTTCGACCa                                                                                               <  1:1493197/150‑1 (MQ=255)                                                             aCCTTGGAGTCGACCGTCTCGATCGACAGCGGCTTGTGCGCCATCAGCGCGTTGGCATCACGCAGCGCATCCATGAAGCTGGTGTAGCGCACATTGACCa                                                                                                                                               >  2:349764/1‑100 (MQ=255)                                                             aCCTTGGAGTCGACCGTCTCGATCGACAGCGGCTTGTGCGCCATCAGCGCGTTGGCATCACGCAGCGCATCCATGAAGCTGGTGTAGCGCACATTGACCa                                                                                                                                               <  1:349764/100‑1 (MQ=255)                                                                ttGGAGTCGACCGTCTCGATCGACAGCGGCTTGTGCGCCATCAGCGCGTTGGCATCACGCAGCGCATCCATGAAGCTGGTGTAGCGCACATTGACCAGCACCGCATACTTGGGAATCGGCAGCACGTTGAGCTTGGCTTCGACCACGTAg                                                                                          >  2:1647786/1‑150 (MQ=255)                                                                    aGTCGACCGTCTCGATCGACAGCGGCTTGTGCGCCATCAGCGCGTTGGCATCACGCAGCGCATCCATGAAGCTGGTGTAGCGCACATTGACCAGCACCGCATACTTGGGAATCGGCAGCACGTTGAGCTTGGCTTCGACCACGTAG‑‑CCCa                                                                                    <  1:1899350/150‑1 (MQ=255)                                                                    aGTCGACCGTCTCGATCGACAGCGGCTTGTGCGCCATCAGCGCGTTGGCATCACGCAGCGCATCCATGAAGCTGGTGTAGCGCACATTGACCAGCACCGCATACTTGGGAAGCGGCAGCACGTTGAGCTTGGCTTCGACCACGTAG‑‑CCCa                                                                                    <  2:163001/150‑1 (MQ=255)                                                                     gTCGACCGTCTCGATCGACAGCGGCTTGTGCGCCATCAGCGCGTTGGCATCACGCAGCGCATCCATGAAGCTGGTGTAGCGCACATTGACCAGCACCGCATACTTGGGAATCGGCAGCACGTTGAGCTTGGCTTCGACCACGTAG‑‑CCCAg                                                                                   >  1:1732582/1‑150 (MQ=255)                                                                        gACCGTCTCGATCGACAGCGGCTTGTGCGCCATCAGCGCGTTGGCATCACGCAGCGCATCCATGAAGCTGGTGTAGCGCACATTGACCAGCACCGCATACTTGGGAATCGGCAGCACGTTGAGCTTGGCTTCGACCACGTAG‑‑CCCAGCGa                                                                                >  1:1346453/1‑150 (MQ=255)                                                                        gACCGTCTCGATCGACAGCGGCTTGTGCGCCATCAGCGCGTTGGCATCACGCAGCGCATCCATGAAGCTGGTGTAGCGCACATTGACCAGCACCGCATACTTGGGAATCGGCAGCACGTTGAGCTTGGCTTCGACCACGTAG‑‑CCCAGCGa                                                                                <  2:843977/150‑1 (MQ=255)                                                                              ctcGATCGACAGCGGCTTGTGCGCCATCAGCCCGTTGGCATGACGCACGGCATCCAGTAAGCTGGTTTAGCGGACATTGACCAGCACCGCATACTTGGGAATCGGCAGCACGTTGAGCTTGGCTTCGACCACGTAG‑‑CCCAGCGA‑GCCCTc                                                                         >  1:1008805/1‑150 (MQ=255)                                                                                 gaccgccAGCTGCCTGTGCGCCATCAGCGCGTTGGCATCACGCAGCGCATCCATGAAGCTGGTGTAGCGCACATTGACCAGCACCGCATACTTGGGAATCGGCAGCACGTTGAGCTTGGCTTCGACCGCGTAG‑‑CCCAGCGA‑GCCCTCGgc                                                                      >  2:483596/7‑150 (MQ=255)                                                                                           ggCTTGTGCGCCATCAGCGCGTTGGCATCACGCAGCGCATCCATGAAGCTGGTGTAGCGCACATTGACCAGCACCGCATACTTGGGAATCGGCAGCACGTTGAGCTTGGCTTCGACCACGTAG‑‑CCCAGCGA‑GCCCTCGGCGCCGCACAGt                                                            >  1:1775533/1‑150 (MQ=255)                                                                                                          aGCGCGTTGGCATCACGCAGCGCATCCATGAAGCTGGTGTAGCGCACATTGACCAGCACCGCATACTTGGGAATCGGCAGCACGTTGAGCTTGGCTTCGACCACGTAG‑‑cc                                                                                      >  2:769051/1‑110 (MQ=255)                                                                                                          aGCGCGTTGGCATCACGCAGCGCATCCATGAAGCTGGTGTAGCGCACATTGACCAGCACCGCATACTTGGGAATCGGCAGCACGTTGAGCTTGGCTTCGACCACGTAG‑‑cc                                                                                      <  1:769051/110‑1 (MQ=255)                                                                                                          aGCGCGTTGGCATCACGCAGCGCATCCATGAAGCTGGTGTAGCGCACATTGACCAGCACCGCATACTTGGGAATCGGCAGCACGTTGAGCTTGGCTTCGACCACGTAG‑‑CCCAGCGA‑GCCCTCGGCGCCGCACAGTACGCTGTTGAGGTTg                                             >  2:857753/1‑150 (MQ=255)                                                                                                            cgcgTTGGCATCACGCAGCGCATCCATGAAGCTGGTGTAGCGCACATTGACCAGCACCGCATACTTGGGAATCGGCAGCACGTTGAGCTTGGCTTCGACCACGTAG‑‑CCCAGCGA‑GCCCTCGGCGCCGCACAGTACGCTGTTGAGGTTGaa                                           <  2:1008805/150‑1 (MQ=255)                                                                                                                    cATCACGCAGCGCATCCATGAAGCTGGTGTAGCGCACATTGACCAGCACCGCATACTTGGGAATCGGCAGCACGTTGAGCTTGGCTTCGACCACGTAG‑‑CCCAGCGA‑GCCCTCGGCGCCGCACAGTACGCTGTTGAGGTTGAAGCGGCCCt                                   <  1:1635745/150‑1 (MQ=255)                                                                                                                     atctcacAGCGCATCCATGAAGCTGGTGTAGCGCACATTGACCAGCACCGCATACTTGGGAATCGGCAGCACGTTGAGCTTGGCTTCGACCACGTAG‑‑CCCAGCGA‑GCCCTCGGCGCCGCACAGTACGCTGTTGAGGTTGAAGCGGCCCTg                                  >  2:485100/7‑150 (MQ=255)                                                                                                                        aCGCAGCGCATCCATGAAGCTGGTGTAGAGCCCATTGACAAGCACCGCATACTTGGGAATCGGCAGCACGTTGAGCTCATGATTGAGAACAAAC‑‑TCCAGCGA‑GCCCTCGGCGCCGCACAGTACGCTGTTGAGGTTGAAGCGGCCCTGCtc                               >  1:925544/1‑150 (MQ=255)                                                                                                                          gCAGCGCATCCATGAAGCCGGTGTAGCGCACATTGACCAGCACCGCATACTTGGGAATCGGCAGCACGTTGAGCTTGGCTTCGACCACGTAG‑‑CCCAGCGA‑GCCCTCGGCGCCGCACAGTACGCTGTTGAGGTTGAAGCGGCCCTGCtcgt                             >  1:1057534/1‑150 (MQ=255)                                                                                                                           cAGCGCATCCATGAAGCTGGTGTAGCGCACATTGACCAGCACCGCATACTTGGGAATCGGCAGCACGTTGAGCTTGGCTTCGACCACGTAG‑‑CCCAGCGA‑GCCCTCGGCGCCGCACAGTACGCTGTTGAGGTTGAAGc                                         <  2:1835156/137‑1 (MQ=255)                                                                                                                           cAGCGCATCCATGAAGCTGGTGTAGCGCACATTGACCAGCACCGCATACTTGGGAATCGGCAGCACGTTGAGCTTGGCTTCGACCACGTAG‑‑CCCAGCGA‑GCCCTCGGCGCCGCACAGTACGCTGTTGAGGTTGAAGc                                         >  1:1835156/1‑137 (MQ=255)                                                                                                                           cAGCGCATCCATGAAGCTGGTGTAGCGCACATTGACCAGCACCGCATACTTGGGAATCGGCAGCACGTTGAGCTTGGCTTCGACCACGTAG‑‑CCCAGCGA‑GCCCTCGGCGCCGCACAGTACGCTGTTGAGGTTGAAGCGGCCCTGCtcgtc                            >  2:1636869/1‑150 (MQ=255)                                                                                                                                     aTGAAGCTGGTGTAGCGCACATTGACCAGCACCGCATACTTGGGAATCGGCAGCACGTTGAGCTTGGCTTCGACCACGTAG‑‑CCCAGCGA‑GCCCTCGGCGCCGCACAGTACGCTGTTGAGGTTGAAGCGGCCCTGCTCGTCGCGCAGGTgc                  >  2:832614/1‑150 (MQ=255)                                                                                                                                     aTGAAGCTGGTGTAGCGCACATTGACCAGCACCGCAGACTTGGGAATCGGCAGCACGTTGAGCTTGGCTTCGAACACGGAG‑‑CTCAGCGA‑GCCCTCGGCCCCGCACGGTACTGTGTTGCGGGGGAAGCGGCCCTGCTCGTCGCGCAgctac                  >  2:40023/1‑146 (MQ=255)                                                                                                                                           cTGGTGTAGCGCACATTGACCAGCACCGCATACTTGGGAATCGGCAGCACGTTGAGCTTGGCTTCGACCACGTAG‑‑CCCAGCGA‑GCCCTCGGCGCCGCACAGTACGCTGTTGAGGTTGAAGCGGCCCTGCTCGTCGCGCAGGTGCGCCAgg            <  1:1216044/150‑1 (MQ=255)                                                                                                                                            tGGCGTAGCGCACATTGACCAGCACCGCACACTTGGGAATCGGCAGCACGTTGAGCTTGGCTTCGACCACGTAG‑‑CCCAGCGA‑GCCCTCGGCGCCGCACAGTACGCTGTTGAGGTTGAAGCGGCCCTGCTCGTCGCGCAGGTGCGCCAGGt           <  1:1857815/150‑1 (MQ=255)                                                                                                                                              gtgtAGCGCACATTGACCAGCACCGCATACTTGGGAATCGGCAGCACGTTGAGCTTGGCTTCGACCACGTAG‑‑CCCAGCGA‑GCCCTCGGCGCCGCACAGTACGCTGTTGAGGTTGAAGCGGCCCTGCTCGTCGCGCAGGTGCGCCAGGTCg         >  1:1421708/1‑150 (MQ=255)                                                                                                                                                gtAGCGCACATTGACCAGCACCGCATACCTGGAAATCGGCAGCACGTTGAGCTTGGCTTCGACGACGTAGATCACAGCCATGGCCACCGCGCCGC‑C‑G‑ATGCTGTTGAGGTTGAAGCGGCCCTGCTCGTCGCGCAGGTGCGCCAGGTCGTa       >  1:1407211/1‑150 (MQ=255)                                                                                                                                                 tAGCGCACATTGACCAGCACCGCATACTTGGGAATCGGCAGCACGTTGAGCTTGGCTTCGACCACGTAG‑‑CCCAGCGA‑GCCCTCGGCGCCGCACAGTACGCTGTTGAGGTTGAAGCGGCCCTGCTCGTAGCGCAGGTGCGCCAGGTCGTAg      >  1:62023/1‑150 (MQ=255)                                                                                                                                                  aGCGCACATTGACCAGCACCGCATACTTGGGAATCGGCAGCACGTTGAGCTTGGCTTCGACCACGTAG‑‑CCCAGCGA‑GCCCTCGGCGCCGCACAGTACGCTGTTGAGGTTGAAGCGGCCCTGCTCGTCGCGCAGGTGCGCCAGGTCGTAGc     <  1:832614/150‑1 (MQ=255)                                                                                                                                                   gcgcACATTGACCAGCACCGCATACTTGGGAATCGGCAGCACGTTGAGCTTGGCTTCGACCACGTAG‑‑CCCAGCGA‑GCCCTCGGCGCCGCACAGTACGCTGTTGAGGTTGAAGCGGCCCTGCTCGTCGCGCAGGTGCGCCAGGTCGTAGcc    >  2:1204134/1‑150 (MQ=255)                                                                                                                                                     gcACATTGACCAGCACCGCATACTTGGGAATCGGCAGCACGTTGAGCTTGGCTTCGACCACGTAG‑‑CCCAGCGA‑GCCCTCGGCGCCGCACAGTACGCTGTTGAGGTTGAAGCGGCCCTGCTCGTCGCGCAGGTGCGCCAGGTCGTAGCCgg  >  1:988593/1‑150 (MQ=255)                                                                                                                                                     |                                                                                                                                                         TCCGCCGGGAAGTACTCGGCAACGCTGTGCCAGACGATGTCTTTCATCGCCAGCATCAGCACCTTGGAGTCGACCGTCTCGATCGACAGCGGCTTGTGCGCCATCAGCGCGTTGGCATCACGCAGCGCATCCATGAAGCTGGTGTAGCGCACATTGACCAGCACCGCATACTTGGGAATCGGCAGCACGTTGAGCTTGGCTTCGACCACGTAG‑‑CCCAGCGA‑GCCCTCGGCGCCGCACAGTACGCTGTTGAGGTTGAAGCGGCCCTGCTCGTCGCGCAGGTGCGCCAGGTCGTAGCCGG  >  NC_002947/5107147‑5107444

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