Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,541,768	T→A	10.9%	I327F (ATC→TTC)	dnaB ←	replicative DNA helicase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,541,768	0	T	A	10.9%	183.5 / 9.0	64	I327F (ATC→TTC)	dnaB	replicative DNA helicase
Reads supporting (aligned to +/- strand):  ref base T (33/24);  new base A (3/4);  total (36/28)
Fisher's exact test for biased strand distribution p-value = 6.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

GCTGCGAAAGGGCAATGACCGGGCAGTTGAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAGTCGTCGTCGTCCAGCT  >  NC_002947/5541621‑5541911                                                                                                                                                    |                                                                                                                                                gCTGCGAAAGGGCAATGACCGGGCAGTTGAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGTAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGCttt                                                                                                                                               <  2:87859/150‑1 (MQ=255) gCTGCGAAAGGGCAATGACCGGGCAGTTGAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGAttt                                                                                                                                               >  2:743752/1‑150 (MQ=255)  cTGCGAAAGGGAAATGACCGGGCAGTTGAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTc                                                                                                                                              <  1:743752/150‑1 (MQ=255)   tGCGAAAGGGCAATGACCGGGCAGTTGAATTCTTTGGCCAGGGCCTTGAGGTAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCg                                                                                                                                             >  1:1531694/1‑150 (MQ=255)   tGCGAAAGGGCAATGACCGGGCAGTTGAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCg                                                                                                                                             >  2:302045/1‑150 (MQ=255)         aGGGCAATGACCGGGCAGTTGAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGAAATCGCCATga                                                                                                                                       <  2:980329/150‑2 (MQ=255)                   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)                                           gggCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTc                                                                                                                                              <  1:1524685/109‑1 (MQ=255)                                           gggCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTc                                                                                                                                              >  2:1524685/1‑109 (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                                                                                                                             >  1:208996/1‑103 (MQ=255)                                                                  tcGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCa                                                                                                                             <  2:208996/103‑1 (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:77156/1‑150 (MQ=255)                                                                                                      cccGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTc                                          >  1:1196365/1‑150 (MQ=255)                                                                                                       ccGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTAtcgt                                                                          >  2:689436/1‑117 (MQ=255)                                                                                                       ccGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTAtcgt                                                                          <  1:689436/117‑1 (MQ=255)                                                                                                             tCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGtt                                           >  2:1067171/1‑142 (MQ=255)                                                                                                             tCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGtt                                           <  1:1067171/142‑1 (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                                                                                                 <  2:1499619/82‑1 (MQ=255)                                                                                                                   tCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGa                                                                                                 >  1:1499619/1‑82 (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                                               >  2:1075048/1‑114 (MQ=255)                                                                                                                                     ccATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCa                                               <  1:1075048/114‑1 (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)                                                                                                                                                    |                                                                                                                                                GCTGCGAAAGGGCAATGACCGGGCAGTTGAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAGTCGTCGTCGTCCAGCT  >  NC_002947/5541621‑5541911

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