Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,679,748	T→G	6.0%	D70A (GAC→GCC)	rnhA ←	ribonuclease H

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,679,748	0	T	G	6.0%	227.4 / 4.4	78	D70A (GAC→GCC)	rnhA	ribonuclease H
Reads supporting (aligned to +/- strand):  ref base T (35/38);  new base G (3/2);  total (38/40)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.30e-01

TCCAGGTTACCTTGTGGCGGTTGACCTGCTCATCGAGCTGCTGCCACAAGTCAGCGTTCTTCACCGGCTCCTTGGCGGCGGTCTTCCAGCCACGCTTCTTCCAGTTGACCATCCATTCGTTGATGCCCTTCATCACGTACTGCGAGTCGGTGGTCAGCACCACCTCGCATTCACGCTTGAGCGACATCAGGCCCTGGATCGCAGCCATCAGCTCCATGCGGTTGTTGGTGGTCTCACGTTCGCCGCCCCACAGTTCCTTCTCGACGCCCTTGTAGATCATCAGCACGCCCCAG  >  NC_002947/4679602‑4679894                                                                                                                                                   |                                                                                                                                                   tCCAGGTTACCTTGTGGCGGTTGACCTGCTCAACGAGCTGCTGCCACAAGTCAGCGTTCTTCACCGGCTCCATGGCGGCGGTCTTCCAGCCCCGCTTCTTCCAGTTGACCATGCATTCCTTGATGCCCTTCATCACGTACTGCGAGTCg                                                                                                                                                  >  7:355758/1‑149 (MQ=255)  ccAGGTTACCTTGTGGCGGTTGACCTGCTCATCGAGCTGCTGCCACAAGTCAGCGTTCTTCACCGGCTCCTTGGCGGCGGTCTTCCAGCCACGCTTCTTCCAGTTGACCATCCATTCGTTGATGCCCTTCATCACGTACTGCGAGTCgg                                                                                                                                                 >  5:51336/1‑149 (MQ=255)      gTTACCTTGTGGCGGTTGACCTGCTCATCGAGCTGCTGCCACAAGTCAGCGTTCTTCACCGGCTCCTTGGCGGCGGTCTTCCAGCCACGCTTCTTCCAGTTGACCATCCATTCGTTGATGCCCTTCATCACGTACTGCGAGTCggtgtt                                                                                                                                             <  2:249696/149‑3 (MQ=255)          ccTTGTGGCGGTTGACCTGCTCATCGAGCTGCTGCCACAAGTCAGCGTTCTTCACCGGCTCCTTGGCGGCGGTCTTCCAGCCACGCTTCTTCCAGTTGACCATCCATTCGTTGATGCCCTTCATCACGTACTGCGAGTCGGTGGTCAg                                                                                                                                          <  4:319430/148‑1 (MQ=255)          ccTTGTGGCGGTTGACCTGCTCATCGAGCTGCTGCCACAAGTCAGCGTTCTTCACCGGCTCCTTGGCGGCGGTCTTCCAGCCACGCTTCTTCCAGTTGACCATCCATTCGTTGATGCCCTTCATCACGTACTGCGAGTCGGTGGTCAGc                                                                                                                                         <  4:90236/149‑1 (MQ=255)           cTTGTGGCGGTTGACCTGCTCATCGAGCTGCTGCCACAAGTCAGCGTTCTTCACCGGCTCCTTGGCGGCGGTCTTCCAGCCACGCTTCTTCCAGTTGACCATCCATTCGTTGATGCCCTTCATCACGTACTGCGAGTCGGTGGTCAGca                                                                                                                                        >  5:253578/1‑149 (MQ=255)           cTTGTGGCGGTTGACCTGCTCATCGAGCTGCTGCCACAAGTCAGCGTTCTTCACCGGCTCCTTGGCGGCGGTCTTCCAGCCACGCTTCTTCCAGTTGACCATCCATTCGTTGATGCCCTTCATCACGTACTGCGAGTCGGTGGTCAGca                                                                                                                                        >  8:149733/1‑149 (MQ=255)                   ggTTGACCTGCTCATCGAGCTGCTGACACAAGTCAGCGTTCTTCACCTGCTCCTTGGCGGCGGTCTTCCAGCCACGCTTCTTCCAGTTGACCATCCATTCGGTGATGCCCTTCCTCACGTACTGCGAGTCGGAGTTCAGCccccccgc                                                                                                                                 >  4:199926/1‑146 (MQ=255)                     ttGACCTGCTCATCGAGCTGCTGCCACAAGTCAGCGTTCTTCACCGGCTCCTTGGCGGCGGTCTTCCAGCCACGCTTCTTCCAGTTGACCATCCATTCGTTGATGCCCTTCATCACGTACTGCGAGTCGGTGGTCAGCACCACCTCGCa                                                                                                                              <  8:53486/149‑1 (MQ=255)                          cTGCTCATCGAGCTGCTGCCACAAGTCAGCGTTCTTCACCGGCTCCTTGGCGGCGGTCTTCCAGCCACGCTTCTTCCAGTTGACCATCCATTCGTTGATGCCCTTCATCACGTACTGCGAGTCGGTGGTCAGCACCACCTCGCATTCAc                                                                                                                         >  2:208282/1‑149 (MQ=255)                          cTGCGCATCGAGCTGCTGCCACAAGTCAGCGTTCTTCACCGGCTCCTTCGCGGCGGTCTTCCAGCCACGCTTCTTCCAGTTGACCATCCATTCGTTGATGCCCTTCATCACGTACTGCGAGTCGGTCGTCAGCACCACCTCGCATTCAc                                                                                                                         <  8:74315/149‑1 (MQ=255)                                aTCGAGCTGCTGCCACAAGTCAGCGTTCTTCACCGGCTCCTTGGCGGCGGTCTTCCAGCCACGCTTCTTCCAGTTGACCATCCATTCGTTGATGCCCTTCATCACGTACTGCGAGTCGGTGGTCAGCACCACCTCGCATTCACGCtt                                                                                                                     >  6:69099/1‑147 (MQ=255)                                   gAGCTGCTGCCACAAGTCAGCGTTCTTCACCGGCTCCTTGGCGGCGGTCTTCCAGCCACGCTTCTTCCAGTTGACCATCCATTCGTTGATGCCCTTCATCACGTACTGCGAGTCGGTGGTCAGCACCACCTCGCATTCACGCTTGAGc                                                                                                                 <  3:167613/148‑1 (MQ=255)                                     gctgctGCCACAAGTCAGCGTTCTTCACCGGCTCCTTGGCGGCGGTCTTCCAGCCACGCTTCTTCCAGTTGACCATCCATTCGTTGATGCCCTTCATCACGTACTGCGAGTCGGTGGTCAGCACCACCTCGCATTCACGCTTGAGCGa                                                                                                               >  4:11818/1‑148 (MQ=255)                                         ctgcCACAAGTCAGCGTTCTTCACCGGCTCCTTGGCGGCGGTCTTCCAGCCACGCTTCTTCCAGTTGACCATCCATTCGTTGATGCCCTTCATCACGTACTGCGAGTCGGTGGTCAGCACCACCTCGCATTCACGCTTGAGCGACATCa                                                                                                          <  6:253578/149‑1 (MQ=255)                                            ccACAAGTCAGCGTTCTTCACCGGCTCCTTGGCGGCGGTCTTCCAGCCACGCTTCTTCCAGTTGACCATCCATTCGTTGATGCCCTTCATCACGTACTGCGAGTCGGTGGTCAGCACCACCTCGCATTCACGCTTGAGCGACATCAg                                                                                                         <  3:357472/147‑1 (MQ=255)                                            ccACAAGTCAGCGTTCTTCACCGGCTCCTTGGCGGCGGTCTTCCAGCCACGCTTCTTCCAGTTGACCATCCATTCGTTGATGCCCTTCATCACGTACTGCGAGGTGGTGCTGACCACCGACTCGCag                                                                                                                             <  8:96831/127‑2 (MQ=255)                                            ccACAAGTCAGCGTTCTTCACCGGCTCCTTGGCGGCGGTCTTCCAGCCACGCTTCTTCCAGTTGACCATCCATTCGTTGATGCCCTTCATCACGTACTGCGAGGTGGTGCTGACCACCGACTCGCag                                                                                                                             >  7:96831/1‑126 (MQ=255)                                              acaAGTCAGCGTTCTTCACCGGCTCCTTGGCGGCGGTCTTCCAGCCACGCTTCTTCCAGTTGACCATCCATTCGTTGATGCCCTTCATCACGTACTGCGAGTCGGTGGTCAg                                                                                                                                          >  3:266352/1‑112 (MQ=255)                                              acaAGTCAGCGTTCTTCACCGGCTCCTTGGCGGCGGTCTTCCAGCCACGCTTCTTCCAGTTGACCATCCATTCGTTGATGCCCTTCATCACGTACTGCGAGTCGGTGGTCAg                                                                                                                                          <  4:266352/112‑1 (MQ=255)                                              acaAGTCAGCGTTCTTCACCGGCTCCTTGGCGGCGGTCTTCCAGCCACGCTTCTTCCAGTTGACCATCCATTCGTTGATGCCCTTCATCACGTACTGCGAGTCGGTGGTCAGCACCACCTCGCATTCACGCTTGAGCGACATCAGGccc                                                                                                     <  5:69099/149‑1 (MQ=255)                                               caAGTCAGCGTTCTTCACCGGCTCCTTGGCGCCGGTCTTCCAGCCACGCTTCTTCCAGTTGACCATCCATTCGTTGATGCCCTTCATCACGTACTGCGAGTCGGTGGTCAGCACCACCTCGCATTCACGCTTGATCGACATCAGGCCCt                                                                                                    >  5:298624/1‑149 (MQ=255)                                                  gTCAGCGTTCTTCACCGGCTCCTTGGCGGCGGTCTTCCAGCCACGCTTCTTCCAGTTGACCATCCATTCGTTGATGCCCTTCATCACGTACTGCGAGTCGGTGGTCAGCACCACCTCGCATTCACGCTTGAGCGACATCAGGCCCTGGa                                                                                                 <  7:11180/149‑1 (MQ=255)                                                  gTCAGCGTTCTTCACCGGCTCCTTGGCGGCGGTCTTCCAGCCACGCTTCTTCCAGTTGACCATCCATTCGTTGATGCCCTTCATCACGTACTGCGAGTCGGTGGTCAGCACCACCTCGCATTCACGCTTGAGCGACATCAGGCCCTGGa                                                                                                 <  1:208282/149‑1 (MQ=255)                                                    cAGCGTTCTTCACCGGCTCCTTGGCGGCGGTCTTCCAGCCACGCTTCTTCCAGTTGACCATCCATTCGTTGATGCCCTTCATCACGTACTGCGAGTCGGTGGTCAGCACCACCTCGCATTCACGCTTGAGCGACATCAGGCCCTGGAt                                                                                                <  7:108279/148‑1 (MQ=255)                                                    cAGCGTTCTTCACCGGCTCCTTGGCGGCGGTCTTCCAGCCACGCTTCTTCCAGTTGACCATCCATTCGTTGATGCCCTTCATCACGTACTGCGAGTCGGTGGTCAGCACCACCTCGCATTCACGCTTGAGCGACATCAGGCCCTGGATc                                                                                               >  7:281109/1‑149 (MQ=255)                                                        gTTCTTCACCGGCTCCTTGGCGGCGGTCTTCCAGCCACGCTTCTTCCAGTTGACCATCCATTCGTTGATGCCCTTCATCACGTACTGCGAGTCGGTGGTCAGCACCACCTCGCATTCACGCTTGATCGACATCAGGCCCTGGATCGCAg                                                                                           <  6:374654/149‑1 (MQ=255)                                                        gTTCTTCACCGGCTCCTTGGCGGCGGTCTTCCAGCCACGCTTCTTCCAGTTGACCATCCATTCGTTGATGCCCTTCATCACGTACTGCGAGTCGGTGGTCAGCACCACCTCGCATTCACGCTTGAGCGACATCAGGCCCTGGATCg                                                                                              <  2:37706/146‑1 (MQ=255)                                                        gTTCTTCACCGGCTCCTTGGCGGCGGTCTTCCAGCCACGCTTCTTCCAGTTGACCATCCATTCGTTGATGACCTTCATCACGTACTGCGAGTCGGTGGTCAGCACCACCTCGCATTCACGCTTGAGCGACATCAGGCCCTGGATCg                                                                                              >  1:37706/1‑146 (MQ=255)                                                           cttcTCCGGCTCCTTGGCGGGGGTCTTCCAGCCACGCTTCTTCCAGTTGACCATCCATTCGTTGATGCCCATCATCACGTACTGCGAGTCGGTGGTCAGCACCACCTCGCATTCACGCTTGAGCGACATCAGGCCCTGGATCGCAGcc                                                                                         <  8:355758/148‑1 (MQ=255)                                                            ttcACCGGCTCCTTGGCGGCGGTCTTCCAGCCACGCTTCTTCCAGTTGACCATCCATTCGTTGATGCCCTTCATCACGTACTGCGAGTCGGTGGTCAGCACCACCTCGCATTCACGCTTGAGCGACATCAGGCCCTGGATCGCAGCCAt                                                                                       >  2:377093/1‑149 (MQ=255)                                                            ttcACCGGCTCCTTGGCGGCGGTCTTCCAGCCACGCTTCTTCCAGTTGACCATCCATTCGTTGATGCCCTTCATCACGTACTGCGAGTCGGTGGTCAGCACCACCTCGCATTCACGCTTGAGCGACATCAGGCCCTGGATCGCAGCCAt                                                                                       >  2:102851/1‑149 (MQ=255)                                                             tcACCGGCTCCTTGGCGGCGGTCTTCCAGCCACGCTTCTTCCAGTTGACCATCCATTCGTTGATGCCCTTCATCACGTACTGCGAGTCGGTGGTCAGCACCACCTCGCATTCACGCTTGAGCGACTTCAGGCCCTGGATCGCAGCCAt                                                                                       >  6:291220/1‑148 (MQ=255)                                                              cACCGGCTCCTTGGCGGCGGTCTTCCAGCCACGCTTCTTCCAGTTGACCATCCATTCGTTGATGCCCTTCATCACGTACTGCGAGTCGGTGGTCAGCACCACCTCGCATTCACGCTTGAGCGACATCAGGCCCTGGATCGCAGCCAt                                                                                       >  1:201697/1‑147 (MQ=255)                                                               aCCGGCTCCTTGGCGGCGGTCTTCCAGCCACGCTTCTTCCAGTTGACCATCCATTCGTTGATGCCCTTCATCACGTACTGCGAGTCGGTGGTCAGCACCACCTCGCATTCACGCTTGAGCGACATCAGGCCCTGGATCGCAGCCATCa                                                                                     <  5:132638/148‑1 (MQ=255)                                                               aCCGGCTCCTTGGCGGCGGTCTTCCAGCCACGCTTCTTCCAGTTGACCATCCATTCGTTGATGCCCTTCATCACGTACTGCGAGTCGGTGGTCAGCACCACCTCGCATTCACGCTTGAGCGACATCAGGCCCTGGATCGCAGCCATCAg                                                                                    <  6:333657/149‑1 (MQ=255)                                                                ccGGCTCCTTGGCGGGGTGCTTCCAGCCACGCTTCTTCCAGTGGACCATACATTCGTTGATGCCCTTCATCACGTACTGCGAGTCGGTGGTCAGCACCACCTCTCATTCACGCTTGAGCGACATCAGGCCCTGGATCGCAGCCATCAGc                                                                                   <  3:199926/149‑1 (MQ=255)                                                                   gCTCCTTGGCGGCGGTCTTCCAGCCACGCTTCTTCCAGTTGACCATCCATTCGTTGATGCCCTTCATCACGTACTGCGAGTCGGTGGTCAGCACCACCTCGCATTCACGCTTGAGCGACATCAGGTCCTGGATCGCAGCCATCAGCTcc                                                                                <  7:91487/149‑1 (MQ=255)                                                                   gCTCCTTGGCGGCGGTCTTCCAGCCACGCTTCTTCCAGTTGACCATCCATTCGTTGATGCCCTTCATCACGTACTGCGAGTCGGTGGTCAGCACCACCTCGCATTCACGCTTGAGCGACATCAGGCCCTGGATCGCAGCCATCAGCTcc                                                                                <  1:102851/149‑1 (MQ=255)                                                                             ggcggTCTTCCAGCCACGCTTCTTCCAGTTGACCATCCATTCGTTGATGCCCTTCATCACGTACTGCGAGTCGGTGGTCAGCACCACCTCGCATTCACGCTTGAGCGACATCAGGCCCTGGATCGCAGCCATCAGCTCCATGCGgttgt                                                                      >  5:90280/1‑149 (MQ=255)                                                                              gcggTCTTCCAGCCACGCTTCTTCCAGTTGACCATCCATTCGTTGATGCCCTTCATCACGTACTGCGAGTCGGTGGTCAGCACCACCTCGCATTCACGCTTGAGCGACATCAGGCCCTGGATCGCAGCCATCAGCTCCATGCGgttgtt                                                                     <  7:92404/149‑1 (MQ=255)                                                                                ggTCTTCCAGCCACGCTTCTTCCAGTTGACCATCCATTCGTTGATGCCCTTCATCACGTACTGCGAGTCGGTGGTCAGCACCACCTCGCATTCACGCTTGAGCGACATCAGGCCCTGGATCGCAGCCATCAGCTCCATGCGGTTGTtgg                                                                   >  7:279313/1‑149 (MQ=255)                                                                                 gTCTTCCAGCCACGCTTCTTCCAGTTGACCATCCATTCGTTGATGCCCTTCATCACGTACTGCGAGTCGGTGGTCAGCACCACCTCGCATTCACGCTTGAGCGACATCAGGCCCTGGATCGCAGCCATCAGCTCCATGCGGTTGTtggt                                                                  <  2:201697/149‑1 (MQ=255)                                                                                  tCTTCCAGCCACGCTTCTTCCAGTTGACCATCCATTCGTTGATGCCCTTCATCACGTACTGCGAGTCGGTGGTCAGCACCACCTCGCATTCACGCTTGAGCGACATCAGGCCCTGGATCGCAGCCATCAGCTCCATGCGGTTGTtggtg                                                                 <  1:377093/149‑1 (MQ=255)                                                                                   cTTCCAGCCACGCTTCTTCCAGTTGACCATCCATTCGTTGATGCCCTTCATCACGTACTGCGAGTCGGTGGTCAGCACCACCTCGCATTCACGCTTGAGCGACATCAGGCCCTGGATCGCAGCCATCAGCTCCAt                                                                              >  4:322913/1‑135 (MQ=255)                                                                                   cTTCCAGCCACGCTTCTTCCAGTTGACCATCCATTCGTTGATGCCCTTCATCACGTACTGCGAGTCGGTGGTCAGCACCACCTCGCATTCACGCTTGAGCGACATCAGGCCCTGGATCGCAGCCATCAGCTCCAt                                                                              <  3:322913/135‑1 (MQ=255)                                                                                      ccAGCCACGCTTCTTCCAGTTGACCATCCATTCGTTGATGCCCTTCATCACGTACTGCGAGTCGGTGGTCAGCACCACCTCGCATTCACGCTTGAGCGACATCAGGCCCTGGATCGCAGCCATCAGCTCCATGCGGTTGTTggtggt                                                               <  7:117707/147‑1 (MQ=255)                                                                                           cACGCTTCTTCCAGTTGACCATACATTCTTTGAAGAACTTCATAACGTACTGCGAGTCGGTGTTAAGCACCAACTAGCATTCACGCTTGAGCGACATCACGAACTAGATCGAAGCCATAATCAACATGCCGTTTTTGGTGGTCTCatag                                                        >  8:299384/1‑146 (MQ=255)                                                                                             cGCTTCTTCCAGTTGACCATCCATTCGTTGATGCCCTTCATCACGTACTGCGAGTCGGTGGTCAGCACCACCTCGCATTCACGCTTGAGCGACATCAGGCCCTGGATCGCAGCCATCAGCTCCATGCGGTTGTTGGTGGTCTCACGTTc                                                      <  8:281109/149‑1 (MQ=255)                                                                                               cttcttCCAGTTGACCATCCATTCGTTGATGCCCTTCATCACGTACTGCGAGTCGGTGGTCAGCACCACCTCGCATTCACGCTTGAGCGACATCAGGCCCTGGATCGCAGCCATCAGCTCCATGCGGTTGTTGGTGGTCTCACGTTcgc                                                    <  6:51336/149‑1 (MQ=255)                                                                                                     ccAGTTGACCATCCCTTCGGTGATGCCCATCAGCACGTACTGCGAGGCGGCGGTCAGCACCACCTCGCATTCACGCGTGAGCGACATCAGGCCCTGGATCGCAGCCATCAGCTCCATGCGGTTGTTGGTGGTCTCACGTTCGCCGccc                                               <  4:227029/148‑1 (MQ=255)                                                                                                     ccAGTTGACCATCCATTCCTTGATGCCCTTCTTCACGTTCTGCGAGTCGGAGGGCAGCACCAGCTCGCATTCACGCTTGAGCGACATCAGTCCCTGGATCTCAGCCATCAGCTGCATGCGGTAGTTTGTGGTCTCACGTTCgccgcc                                                <  6:298624/147‑1 (MQ=255)                                                                                                      cAGTTGACCATCCATTCGTTGATGCCCTTCATCACGTACTGCGAGTCGGTGGTCAGCACCACCTCGCATTCACGCTTGAGCGACATCAGGCCCTGGATCGCAGCCATCAGCTCCATGCGGATGTTGGTGGTCTCACGTTCGCCGCCCca                                             >  1:65457/1‑149 (MQ=255)                                                                                                         ttGACCATCCATTCGTTGATGCCCTTCATCACGTACTGCGAGTCGGTGGTCAGCACCACCTCGCATTCACGCTTGAGGGACATCAGGCCCTGGATCGCAGCCATCCGGTCCATGCGGTTGTTGGTGGTCTCACGTTCGCCGCCCCACAg                                          >  6:303240/1‑149 (MQ=255)                                                                                                           gACCATCCATTCGTTGATGCCCTTCATCACGTACTGCGAGTCGGTGGTCAGCACCACCTCGCATTCACGCTTGAGCGACATCAGGCCCTGGATCGCAGCCATCAGCTCCATGCGGTTGTTGGTGGTCTCACGTTCGCCGCCCCACAGt                                         <  5:291220/148‑1 (MQ=255)                                                                                                           gACCATCCATTCGTTGATGCCCTTCATCACGTACTGCGAGTCGGTGGTCAGCACCACCTCGCATTCACGCTTGAGCGACATCAGGCCCTGGATCGCAGCCATCAGCTCCATGCGGTTGTTGGTGGTCTCACGTTCGCCGCCCCACAGt                                         <  5:182108/148‑1 (MQ=255)                                                                                                           gACCATCAATTCGTTGATGCCCTTCATAACGTACTGCGAGTCGGTGGTAAGAACCACCTCGCATTAACGCTTGAGCGACATCAGGCCCTGGTTCGAAGCCATCAGCTCCATGCGGTTGATGGTGGTCTCACGTTCGCCTCCCCACAGt                                         >  5:59865/1‑148 (MQ=255)                                                                                                              catccatTCGTTGATGCCCTTCATCACGTACTGCGAGTCGGTGGTCAGCACCACCTCGCATTCACGCTTGAGCGACATCAGGCCCTGGATCGCAGCCATCAGCTCCATGCGGTTGTTGGTGGTCTCACGTTCGCCGCCCCACAGTTCCt                                     >  4:239569/1‑149 (MQ=255)                                                                                                              catccatTCGTTGATGCCCTTCATCACGTACTGCGAGTCGGTGGTCAGCACCACCTCGCATTCACGCTTGAGCGACATCAGGCCCTGGATCGCAGCCATCAGCTCCATGCGGTTGTTGGTGGTCTCACGTTCGCCGCCCCACAGTTCCt                                     >  4:244418/1‑149 (MQ=255)                                                                                                                   atTCGTTGATGCCCTTCATCACGTACTGCGAGTCGGTGGTCAGCACCACCTCGCATTCACGCTTGAGCGACATCAGGCCCTGGATCGCAGCCATCAGCTCCATGCGGTTGTTGGTGGTCTCACGTTCGCCGCCCCACAGTTCCTtctc                                 >  8:170832/1‑148 (MQ=255)                                                                                                                           aTGCCCTTCATCACGTACTGCGAGTCGGTGGTCAGCACAACCTCGCATTCACGCTTGAGCGACATCAGGACCTGCATCGCAGCCATCAGCTCCATGCGGTTGTTGGTGGTCTCACGTTCGCCGCCCCCCAGTTCCTTCTCGACGCCCtt                        >  3:24408/1‑149 (MQ=255)                                                                                                                            tGCCCTTCATCACGTACTGCGAGTCGGTGGTCAGCACCACCTCGCATTCACGCTTGAGCGACATCAGGCCCTGGATCGCAGCCATCAGCTCCATGCGGTTGTTGGTGGTCTCACGTTCGCCGCCCCACAGTTCCTTCTCGACGCCCtt                        <  7:170832/148‑1 (MQ=255)                                                                                                                              cccTTCATCACGTACTGCGAGTCGGTGGTCAGCACCACCTCGCATTCACGCTTGAGCGACATCAGGCCCTGGATCGCAGCCATCAGCTCCATGCGGTTGTTGGTGGTCTCACGTTCGCCGCCCCACAGTTCCTTCTCGACGCCCTTGt                      <  7:149733/148‑1 (MQ=255)                                                                                                                               ccTTCATCACGTACTGCGAGTCGGTGGTCAGCACCACCTCGCATTCACGCTTGAGCGACATCAGGCCCTGGATCGCAGCCATCAGCTCCATGCGGTTGTTGGTGGTCTCACGTTCGCCGCCCCACAGTTCCTTCTCGACGCCCTTg                       >  7:83384/1‑146 (MQ=255)                                                                                                                               ccTTCATCACGTACTGCGAGTCGGTGGTCAGCACCACCTCGCATTCACGCTTGAGCGACATCAGGCCCTGGATCGCAGCCATCAGCTCCATGCGGTTGTTGGTGGTCTCACGTTCGCCGCCCCACAGTTCCTTCTCGACGCCCTTg                       >  4:61814/1‑146 (MQ=255)                                                                                                                               ccTTCATCACGTACTGCGAGTCGGGTGTCAGCACCACCTCGCATTCACGCTTGAGCGACATCAGGCCCTGGATCGCAGCCATCAGCTCCATGCGGTTGTTGGTGGTCTCACGTTGGCCGCCCCACAGTTCCTTCTCGACGCCCTTGTAg                    <  6:59865/149‑1 (MQ=255)                                                                                                                                   catcaCGTACTGCGAGTCGGTGGTCAGCACCACCTCGCATTCACGCTTGAGCGACATCAGGCCCTGGATCGCAGCCATCAGCTCCATGCGGTTGTTGGTGGTCTCACGTTCGCCGCCCCACAGTTCCTTCTCGACGCCCTTGTAGatca                >  1:258099/1‑149 (MQ=255)                                                                                                                                    atcaCGTACTGCGAGTCGGTGGTCAGCACCACCTCGCATTCACGCTTGAGCGACATCAGGCCCTGGATCGCAGCCATCAGCTCCATGCGGTTGTTGGTGGTCTCACGTTCGCCGCCCCACAGTTCCTTCTCGACGCCCTTGTAGatcat               <  8:279313/149‑1 (MQ=255)                                                                                                                                     tcaCGTACTGCGAGTCGGTGGTCAGCACCACCTCGCATTCACGCTTGAGCGACATCAGGCCCTGGATCGCAGCCATCAGCTCCATGCGGTTGTTGGTGGTCTCACGTTCGCCGCCCCACAGTTCCTTCTCGACGCCCTTGTAGatcatc              <  3:61814/149‑1 (MQ=255)                                                                                                                                      caCGTACTGCGAGTCGGTGGTCAGCACCACCTCGCATTCACGCTTGAGCGACATCAGGCCCTGGATCGCAGCCATCAGCTCCATGCGGTTGTTGGTGGTCTCACGTTCGCCGCCCCACAGTTCCTTCTCGACGCCCTTGTAGatcat               >  4:8181/1‑147 (MQ=255)                                                                                                                                        cGTACTGCGAGTCGGTGGTCAGCACCACCTCGCATTCACGCTTGAGCGACATCAGGCCCTGGATCGCAGCCATCAGCTCCATGCGGTTGTTGGTGGTCTCACGTTCGCCGCCCCACAGTTCCTTCTCGACGCCCTTGTAGATCATCAGc           <  6:127226/149‑1 (MQ=255)                                                                                                                                         tgaatGCGAGGTGGTGCTGACCACCGACTCGCATTCACGCTTGAGCGACATCAGGCCCTGGATCGCAGCCATCAGCTCCATGCGGTTGTTGGTGGTCTCACGTTCGCCGCCCCACAGTTCCTTCTCGACGCCCTTGTAGATCATCAGCa          >  8:30949/5‑149 (MQ=255)                                                                                                                                         gTACTGCGAGTCGGTGGTCAGCACCACCTCGCATTCACGCTTGAGCGACATCAGGCCCTGGATCGCAGCCATCAGCTCCATGCGGTTGTTGGTGGTCTCACGTTCGCCGCCCCACAGTTCCTtct                                  <  6:41264/125‑1 (MQ=255)                                                                                                                                         gTACTGCGAGTCGGTGGTCAGCACCACCTCGCATTCACGCTTGAGCGACATCAGGCCCTGGATCGCAGCCATCAGCTCCATGCGGTTGTTGGTGGTCTCACGTTCGCCGCCCCACAGTTCCTtct                                  >  5:41264/1‑125 (MQ=255)                                                                                                                                           aCTGCGAGTCGGTGGTCAGCACCACCTCGCATTCACGCTTGAGCGACATCAGGCCCTGGATCGCAGCCATCAGCTCCATGCGGTTGTTGGTGGTCTCACGTTCGCCGCCCCACAGTTCCTTCTCGACGCCCTTGTAGATCATCAGCACg        >  1:250575/1‑149 (MQ=255)                                                                                                                                              gCGAGGTGGTGCTGACCACCGACTCGCATTCACGCTTGAGCGACATCAGGCCCTGGATCGCAGCCATCAGCTCCATGCGGTTGTTGGTGGTCTCACGTTCGCCGCCCCACAGTTCCTTCTCGACGCCCTTGTAGATCATCAGCACGccc     >  3:332596/1‑149 (MQ=255)                                                                                                                                               cGAGTCGGTGGTCAGCACCACCTCGCATTCACGCTTGAGCGACATCAGGCCCTGGATCGCAGCCATCAGCTCCATGCGGTTGTTGGTGGTCTCACGTTCGCCGCCCCACAGTTCCTTCTCGACGCCCTTGTAGATCATCAGCACGccc     >  5:252833/1‑148 (MQ=255)                                                                                                                                                 aGTCGGTGGTCAGCACCACCTCGCATTCACGCTTGAGCGACATCAGGCCCTGGATCGCAGCCATCAGCTCCATGCGGTTGTTGGTGGTCTCACGTTCGCCGCCCCACAGTTCCTTCTCGACGCCCTTGTAGATCATCAGCACGCCCCAg  >  6:116262/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                   TCCAGGTTACCTTGTGGCGGTTGACCTGCTCATCGAGCTGCTGCCACAAGTCAGCGTTCTTCACCGGCTCCTTGGCGGCGGTCTTCCAGCCACGCTTCTTCCAGTTGACCATCCATTCGTTGATGCCCTTCATCACGTACTGCGAGTCGGTGGTCAGCACCACCTCGCATTCACGCTTGAGCGACATCAGGCCCTGGATCGCAGCCATCAGCTCCATGCGGTTGTTGGTGGTCTCACGTTCGCCGCCCCACAGTTCCTTCTCGACGCCCTTGTAGATCATCAGCACGCCCCAG  >  NC_002947/4679602‑4679894

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