Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,541,778	C→A	8.2%	E323D (GAG→GAT)	dnaB ←	replicative DNA helicase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,541,778	0	C	A	8.2%	191.6 / 3.1	65	E323D (GAG→GAT)	dnaB	replicative DNA helicase
Reads supporting (aligned to +/- strand):  ref base C (31/27);  new base A (3/3);  total (35/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

CCGGGCAGTTGAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAGTCGTCGTCGTCCAGCTGGCCGGAACGGACCT  >  NC_002947/5541639‑5541926                                                                                                                                            |                                                                                                                                                     ccGGGCAGTTGAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCa                                                                                                                                            <  2:1631829/150‑1 (MQ=255)       aGTTGAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAggcggc                                                                                                                                      <  2:1631689/150‑1 (MQ=255)         ttGAATTCTTTGGCCAGGGTCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGcggcgg                                                                                                                                    <  2:224992/150‑1 (MQ=255)         ttGAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCATGGGCCAGGcggcgg                                                                                                                                    >  1:555836/1‑150 (MQ=255)         ttGAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGcggcgg                                                                                                                                    >  2:1272846/1‑150 (MQ=255)            aaTTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTa                                                                                                                                 >  2:907810/1‑150 (MQ=255)                 tttGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCACGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGGTTTCGCCATGCTCACGGGCCACGCGGCGGGTAcgcgc                                                                                                                            >  2:777900/1‑150 (MQ=255)                      ccAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCATGGGCCAGGCGGCGGGTACGCGCGCGCa                                                                                                                       <  2:555836/150‑1 (MQ=255)                         gggCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGAAATCGCCATGATCatg                                                                                                                                                 <  2:1519051/121‑3 (MQ=255)                         gggCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTAACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGAAATCGCCATGATCatg                                                                                                                                                 >  1:1519051/1‑119 (MQ=255)                              ttGAGCGAGCGGGAAATATCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGAGCATGGCGAAATCGCCATGATCATggtcga                                                                                                                                            <  2:858927/121‑5 (MQ=38)                              ttGAGCGAGCGGGAAATATCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGAGCATGGCGAAATCGCCATGATCATggtcga                                                                                                                                            >  1:858927/1‑117 (MQ=38)                                        gggAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCt                                                                                                          <  2:962217/145‑1 (MQ=255)                                        gggAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCt                                                                                                          >  1:962217/1‑145 (MQ=255)                                                tcGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCCACCCGGTTTGCCTATGCTCACTTTCCAGGCGGCGGGCACGCGCGCGCATTTCCGAGTGGCTGATACCCGCCGTa                                                                                             >  1:849653/1‑150 (MQ=255)                                                tcGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCa                                                                                                                                            <  2:208996/103‑1 (MQ=255)                                                tcGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCa                                                                                                                                            >  1:208996/1‑103 (MQ=255)                                                tcGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTa                                                                                             >  2:1515779/1‑150 (MQ=255)                                                tcGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTa                                                                                             >  2:1552485/1‑150 (MQ=255)                                                   gAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTAtcg                                                                                          >  2:1537935/1‑150 (MQ=255)                                                         tCGTTGGTACGGTCGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATg                                                                                    <  2:1531694/150‑1 (MQ=255)                                                           tttGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGaa                                                                                  <  2:1164303/149‑1 (MQ=255)                                                            ttGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAAc                                                                                 >  1:271788/1‑150 (MQ=255)                                                               gTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAAt                                                                              >  2:1474167/1‑150 (MQ=255)                                                                tACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTTCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAAtt                                                                             <  1:428678/150‑1 (MQ=255)                                                                tACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTAATGCAGGTGTACTGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAAtt                                                                             >  2:875525/1‑150 (MQ=255)                                                                   ggTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTc                                                                          >  2:1085820/1‑150 (MQ=255)                                                                          aCCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGTGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGACGTCTCGTCGATGAACAATTTTCGGTCGtt                                                                   >  2:380006/1‑150 (MQ=255)                                                                          aCCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGtt                                                                   <  1:907810/150‑1 (MQ=255)                                                                          aCCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGtt                                                                   <  2:271788/150‑1 (MQ=255)                                                                                    cccGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTc                                                         >  1:1196365/1‑150 (MQ=255)                                                                                    cccGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTc                                                         >  1:77156/1‑150 (MQ=255)                                                                                     ccGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTAtcgt                                                                                         <  1:689436/117‑1 (MQ=255)                                                                                     ccGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTAtcgt                                                                                         >  2:689436/1‑117 (MQ=255)                                                                                           tCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGtt                                                          <  1:1067171/142‑1 (MQ=255)                                                                                           tCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGtt                                                          >  2:1067171/1‑142 (MQ=255)                                                                                            cTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGa                                                 <  1:1537935/150‑1 (MQ=255)                                                                                             tGCATGAGCTGGAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGAGCGCGCATTTGCGAGGGGCTGATACCCGCCGTATCGTCGATGAAGAAGTTTCGGTTGTTGAGCAGGTTCACTGACGAg                                                >  1:582784/1‑150 (MQ=255)                                                                                             tGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAg                                                <  2:1196365/150‑1 (MQ=255)                                                                                               cATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGt                                              >  1:1560021/1‑150 (MQ=255)                                                                                                 tCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGa                                                                                                                >  1:1499619/1‑82 (MQ=255)                                                                                                 tCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGa                                                                                                                <  2:1499619/82‑1 (MQ=255)                                                                                                  cAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCgcagtaag                                           <  2:684120/150‑8 (MQ=255)                                                                                                         aGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAg                                           >  2:991581/1‑143 (MQ=255)                                                                                                         aGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAg                                           <  1:991581/143‑1 (MQ=255)                                                                                                                  aCCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAgtcgtcg                           >  2:149056/1‑150 (MQ=255)                                                                                                                   ccATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCa                                                              <  1:1075048/114‑1 (MQ=255)                                                                                                                   ccATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCa                                                              >  2:1075048/1‑114 (MQ=255)                                                                                                                    cATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGGGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAgtcgtcgtc                         >  1:857314/1‑150 (MQ=255)                                                                                                                    cATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAgtcgtcgtc                         >  2:371249/1‑150 (MQ=255)                                                                                                                    cATGATCATGGCGATTTCACCATGCTCACGGGCCAGGCGGCGGGCACGCGCGCGCATCACTGAGCCGCTGAGAGCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAgtcgtcgtc                         <  1:777900/150‑1 (MQ=255)                                                                                                                    cATGATCATGGCGATGTCGCCATGCTCACGGGCCAGGGGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAgtcgtcgtc                         <  1:1474167/150‑1 (MQ=255)                                                                                                                     gtgagcATGGCGAAATCGCCATGATCATGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGgc                                         >  2:1637053/6‑133 (MQ=255)                                                                                                                     gtgagcATGGCGAAATCGCCATGATCATGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGgc                                         <  1:1637053/128‑1 (MQ=255)                                                                                                                     aTGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAgtcgtcgtcg                        >  2:479123/1‑150 (MQ=255)                                                                                                                       gATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAgtcgtcgtcgtc                      >  1:494537/1‑150 (MQ=255)                                                                                                                          cATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAGTCGTCGTCGTCCAg                   >  2:737970/1‑150 (MQ=255)                                                                                                                            tGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAGTCGTCGTCGTCCAGCt                 <  1:302045/150‑1 (MQ=255)                                                                                                                                   ttCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAGTCGTCGTCGTCCAGCTGGCCGGa          <  2:1560021/150‑1 (MQ=255)                                                                                                                                       ccATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAGTCGTCGTCGTCCAGCTGGCCGGaa         <  2:1133986/147‑1 (MQ=255)                                                                                                                                       ccATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAGTCGTCGTCGTCCAGCTGGCCGGaa         >  1:1133986/1‑147 (MQ=255)                                                                                                                                        cATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGAt                                                                                     >  2:804866/1‑70 (MQ=255)                                                                                                                                        cATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGAt                                                                                     <  1:804866/70‑1 (MQ=255)                                                                                                                                        cATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAGTCGTCGTCGTCCAGCTGGCCGGAACGGa     <  1:1272846/150‑1 (MQ=255)                                                                                                                                           gCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGGGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAGTCGTCGTCGTCCAGCTGGCCGGAACGGACCt  <  2:857314/150‑1 (MQ=255)                                                                                                                                            |                                                                                                                                                     CCGGGCAGTTGAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAGTCGTCGTCGTCCAGCTGGCCGGAACGGACCT  >  NC_002947/5541639‑5541926

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: ‑