Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,767,352	T→G	5.9%	M1174L (ATG→CTG)	PP_5060 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,767,352	0	T	G	5.9%	215.8 / 3.7	68	M1174L (ATG→CTG)	PP_5060	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (30/33);  new base G (2/2);  total (33/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

AGGCAAGTCACAGAACATCAGCAGATGATGCGTTTGTGCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCG‑ATGTCGGTGATGGCCAACAGGCGATCCAGCTGGCTGGAGACATTGATGATC  >  NC_002947/5767212‑5767500                                                                                                                                             |                                                                                                                                                      aGGCAAGTCACAGAACATCAGCAGATGATGCGTTTGTGCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCATGGGCCTggt                                                                                                                                               <  1:153526/149‑2 (MQ=255)   gCAAGTCACAGAACATCAGCAGATGATGCGTTTGTGCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGaa                                                                                                                                              >  5:156711/1‑148 (MQ=255)     aaGTCACAGAACATCAGCAGATGATGCGTTTGTGCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGaaaaa                                                                                                                                           >  6:1082/1‑149 (MQ=255)      aGTCACAGAACATCAGCAGATGATGCGTTTGTGCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAAc                                                                                                                                          <  3:245494/149‑1 (MQ=255)          acaGAACATCAGCAGATGATGCGTTTGTGCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCgg                                                                                                                                       <  7:245311/148‑1 (MQ=255)          acaGAACATCAGCAGATGATGCGTTTGTGCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCgg                                                                                                                                       <  7:135619/148‑1 (MQ=255)          acaGAACATCAGCAGATGATGCGTTTGTGCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGa                                                                                                                                      <  2:76696/149‑1 (MQ=255)            aGAACATCAGCAGATGATGCGTTTGTGCGGCCAACCGCAGACGCTGACTGGCGAGAGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTTTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCATGGAAAAACCGGATg                                                                                                                                    >  4:145556/1‑149 (MQ=255)             gAACATCAGCAGATGATGCGTTTGTGCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGTAAAACCGGATGt                                                                                                                                   <  4:74902/149‑1 (MQ=255)                     gcagATGATGCGTTTGTGCGGCCACCCGCAGACGCTGACTGGCGAGAGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTg                                                                                                                           >  3:118160/1‑149 (MQ=255)                      cagATGATGCGTTTGTGCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTgg                                                                                                                          >  6:267733/1‑149 (MQ=255)                      cagATGATGCGTTTGTGCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCATGGGCCTggt                                                                                                                                               <  4:145992/128‑2 (MQ=255)                      cagATGATGCGTTTGTGCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCATGGGCCTggt                                                                                                                                               >  3:145992/1‑127 (MQ=255)                      cagAGGATGCGTTTGTGCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATTTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATc                                                                                                                              <  8:150309/145‑1 (MQ=255)                        gatgatGCGTTTGTGCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGgtc                                                                                                                        >  7:37059/1‑149 (MQ=255)                        gatgatGCGTTTGTGCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGgtc                                                                                                                        <  1:138720/149‑1 (MQ=255)                               cGTTTGTGCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTa                                                                                                                  <  4:190813/148‑1 (MQ=255)                                    gtgCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGATGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTc                                                                                                            >  1:152297/1‑149 (MQ=255)                                     tgCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCg                                                                                                                                        >  4:208349/1‑120 (MQ=255)                                     tgCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCg                                                                                                                                        <  3:208349/120‑1 (MQ=255)                                        ggCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTg                                                                                                                                                 <  8:25147/108‑1 (MQ=255)                                        ggCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTg                                                                                                                                                 >  7:25147/1‑108 (MQ=255)                                        ggCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGa                                                                                                        >  4:129275/1‑149 (MQ=255)                                        ggCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGa                                                                                                        >  4:29459/1‑149 (MQ=255)                                           cACCCGCAGACGCTGACTGGCGATCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGTAAAAACCGGAGGTTGATCGTGGGCGTCTTACGGCTCGCGATaa                                                                                                     <  3:145556/149‑1 (MQ=255)                                              ccGCAGATGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGc                                                                                                  <  2:19252/149‑1 (MQ=255)                                                 cAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTcgc                                                                                                          >  3:278522/1‑138 (MQ=255)                                                 cAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTcgc                                                                                                          <  4:278522/138‑1 (MQ=255)                                                 cAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCg                                                                                               <  3:29459/149‑1 (MQ=255)                                                  aGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCga                                                                                              <  5:76554/149‑1 (MQ=255)                                                      gCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACg                                                                                          >  4:70083/1‑149 (MQ=255)                                                               cGAGCGCATGGCGGATGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTg                                                                                 <  2:152297/149‑1 (MQ=255)                                                                gAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTgg                                                                                >  4:98903/1‑149 (MQ=255)                                                                gAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTgg                                                                                >  2:178964/1‑149 (MQ=255)                                                                           ggAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTaa                                                                     <  5:267733/149‑1 (MQ=255)                                                                              aGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGTATGTTGATCTTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAgg                                                                   <  4:118160/148‑1 (MQ=255)                                                                              aGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGc                                                                  <  8:37059/149‑1 (MQ=255)                                                                                         tGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGAACTTGCatgtc                                                       >  3:107965/1‑149 (MQ=255)                                                                                                 aGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCG‑ATGTCgg                                              >  8:119042/1‑149 (MQ=255)                                                                                                  ggCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCG‑ATGTCgg                                              <  3:129275/148‑1 (MQ=255)                                                                                                   gCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCG‑ATGTCGGt                                             >  7:209004/1‑148 (MQ=255)                                                                                                   gCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCG‑ATGTCGGTg                                            >  1:291958/1‑149 (MQ=255)                                                                                                     gTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGAc                                                                         <  6:287022/119‑1 (MQ=255)                                                                                                     gTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGAc                                                                         >  5:287022/1‑119 (MQ=255)                                                                                                     gTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCG‑ATGTCGGt                                             >  2:230603/1‑146 (MQ=255)                                                                                                     gTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCG‑ATGTCGGt                                             <  1:230603/146‑1 (MQ=255)                                                                                                     gTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCG‑ATGTCGGTGAt                                          <  1:178964/149‑1 (MQ=255)                                                                                                      tGCCTTAGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCCTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTTCATGTAG‑ATGTCGGTGATg                                         >  5:14082/1‑149 (MQ=255)                                                                                                         cTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCG‑ATGTCGGTGATGGc                                       >  5:319254/1‑148 (MQ=255)                                                                                                               gCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCG‑ATGTCGGTGATGGCCAACAgg                                >  5:7343/1‑149 (MQ=255)                                                                                                               gCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCG‑ATGTCGGTGATGGCCAACAgg                                >  1:240573/1‑149 (MQ=255)                                                                                                                ccACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCG‑ATGTCGGTGATGGCCAACAgg                                <  6:156711/148‑1 (MQ=255)                                                                                                                ccACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCG‑ATGTCGGTGATGGCCAACAgg                                <  2:240573/148‑1 (MQ=255)                                                                                                                ccACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCG‑ATGTCGGTGATGGCCAACAGGc                               >  1:65664/1‑149 (MQ=255)                                                                                                                ccACTGCAGGGCCAGCTCGACCAGGCCCAAGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCG‑ATGTCGGTGATGGCCAACAGGc                               >  2:195867/1‑149 (MQ=255)                                                                                                                 cACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCG‑ATGTCGGTGATGGCCAACAGGc                               >  6:296402/1‑148 (MQ=255)                                                                                                                  aCTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCG‑ATGTCGGTGATGGCCAACAGGCGa                             <  3:98903/149‑1 (MQ=255)                                                                                                                         ggCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCG‑ATGTCGGTGATGGCCAACAGGCGATCCAgct                      >  5:68834/1‑149 (MQ=255)                                                                                                                          ggcAGCTCGACCACGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTGTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCGTATGTCGGTGATGGCCAACAGGCGATCCAgct                      <  1:195867/147‑1 (MQ=255)                                                                                                                             aGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCG‑ATGTCGGTGATg                                         <  6:204181/126‑1 (MQ=255)                                                                                                                             aGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCG‑ATGTCGGTGATg                                         >  5:204181/1‑126 (MQ=255)                                                                                                                                  gACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCTCTTGCATGTCG‑ATGTCGGTTATGGCCAACAGGCGATCCAGCTGGCTGga               <  1:188286/147‑1 (MQ=255)                                                                                                                                  gACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCG‑ATGTCGGTGATGGCCAACAGGCGATCCAGCTGGCTGga               >  2:188286/1‑147 (MQ=255)                                                                                                                                   tccAGGCCATGGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCG‑ATGTCGGTGATGGCCAACAGGCGATCCAGCTGGCTGgaga             >  1:301079/2‑148 (MQ=255)                                                                                                                                    ccAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCG‑ATGTCGGTGATGGCCAACAGGCGATCCAGCTGGCTGGAGACa           <  7:119042/149‑1 (MQ=255)                                                                                                                                        gCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCG‑ATGTCGGTGATGGCCAACAGGCGATCCAGCTGGCTGGAGACATtga       <  8:209004/149‑1 (MQ=255)                                                                                                                                        gCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCG‑ATGTCGGTGATGGCCAACAGGCGATCCAGCTGGCTGGAGACATtga       <  2:291958/149‑1 (MQ=255)                                                                                                                                        gCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCGCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCG‑ATGTCGGTGATGGCCAACAGGCGATCCAGCTGGCTGGAGACATtga       <  5:1082/149‑1 (MQ=255)                                                                                                                                             tGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCG‑ATGTCGGTGATGGCCAACAGGCGATCCAGCTGGCTGGAGACATTGATGATc  <  2:65664/149‑1 (MQ=255)                                                                                                                                             |                                                                                                                                                      AGGCAAGTCACAGAACATCAGCAGATGATGCGTTTGTGCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCG‑ATGTCGGTGATGGCCAACAGGCGATCCAGCTGGCTGGAGACATTGATGATC  >  NC_002947/5767212‑5767500

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