Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,437,506	A→G	5.5%	I18T (ATC→ACC)	PP_4775 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,437,506	0	A	G	5.5%	216.9 / 2.0	73	I18T (ATC→ACC)	PP_4775	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (34/35);  new base G (2/2);  total (36/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.00e-01

TCCTGGTAATACCCCTGGCGCGGCGGTTGGGTCTGGCGCACTTCGTCGCGCGAGTTGAGCGAGCGCCCCGGTCCACCCTCACCTTGCACCTGCTGGCCACGTTCACCGTGCCCCTCGCCGGGCCCTTCAGCCAAGGCCCCAAGGCTGATGCTCAAACCCAGTAAACCAACACCTGCCAACTGCCAGATGCGCGACGTCATGGTGTTCCTCTCGGTAGAAAAAACTGCAACCAACGCCAGTCTACCGTGCCCACCACACAGGCCAATTAAAATTCCCGACATGAAAAAAGGGA  >  NC_002947/5437359‑5437650                                                                                                                                                    |                                                                                                                                                 tCCTTGTAATTCCCCTGGCGCGGCGGTTGGGTCTGGCGCACTTCGTCGCGCGAGTTGAGCGAGCGCCCCGGTCCACCCTCACCTTGCACCTGCTGGCCACGTTCACCGTGCCCCTCGCCGGGCCCTTCAGCCAAGGCCCCAAGGCTGAt                                                                                                                                                 <  8:131366/149‑1 (MQ=255) tCCTGGTAATACCCCTGGCGCGGCGGTTGGGTCTGGCGCACTTCGTCGCGCGAGTTGAGCGAGCGCCCCGGTCCACCCTCACCTTGCACCTGCTGGCCACGTTCACCGTGCCCCTCGCCGGGCCCTTCAGCCAAGGCCCCAAGGCTGAt                                                                                                                                                 <  3:131633/149‑1 (MQ=255)    tGGTAATACCCCTGGCGCGGCGGTTGGGTCTGGCGCACTTCGTCGCGCGAGTTGAGCGAGCGCCCCGGTCCACCCTCACCTTGCACCTGCTGGCCACGTTCACCGTGCCCCTCGCCGGGCCCTTCAGCCAAGGCCCCAAGGCTGATGCt                                                                                                                                              >  2:345187/1‑149 (MQ=255)      gTAATACCCCTGGCGCGGCGGTTGGGTCTGGCGCACTTCGTCGCGCGAGTTGAGCGAGCGCCCCGGTCCACCCTCACCTTGCACCTGCTGGCCACGTTCACCGTGCCCCTCGCCGGGCCCTTCAGCCAAGGCCCCAAGGCTGATGCt                                                                                                                                              <  2:63348/147‑1 (MQ=255)      gTAATACCCCTGGCGCGGCGGTTGGGTCTGGCGCACTTCGTCGCGCGAGTTGAGCGAGCGCCCCGGTCCACCCTCACCTTGCACCTGCTGGCCACGTTCACCGTGCCCCTCGCCGGGCCCTTCAGCCAAGGCCCCAAGGCTGATGCTc                                                                                                                                             <  1:7519/148‑1 (MQ=255)         aTACCCCTGGCGCGGCGGTTGGGTCTGGCGCACTTCGTCGCGCGAGTTGAGCGAGCGCCCCGGTCCACCCTCACCTTGCACCTGCTGGCCACGTTCACCGTGCCCCTCGCCGGGCCCTTCAGCCAAGGCCCCAAGGCTGATGCTCAAAc                                                                                                                                         <  4:120706/149‑1 (MQ=255)         aTACCCCTGGCGCGGCGGTTGGGTCTGGCGCACTTCGTCGCGCGAGTTGAGCGAGCGCCCCGGTCCACCCTCACCTTGCACCTGCTGGCCACGTTCACCGTGCCCCTCGCCGGGCCCTTCAGCCAAGGCCCCAAGGCTGATGCTCAAAc                                                                                                                                         >  1:228397/1‑149 (MQ=255)            ccccTGGCGCGGCGGTTGGGTCTGGCGCACTTCGTCGCGCGAGTTGAGCGAGCGCCCCGGTCCACCCTCACCTTGCACCTGCTGGCCACGTTCACCGTGCCCCTCGCCGGGCCCTTCAGCCAAGGCCCCAAGGCTGATGCTCaaa                                                                                                                                          <  1:221215/145‑1 (MQ=255)            ccccTGGCGCGGCGGTTGGGTCTGGCGCACTTCGTCGCGCGAGTTGAGCGAGCGCCCCGGTCCACCCTCACCTTGCACCTGCTGGCCACGTTCACCGTGCCCCTCGCCGGGCCCTTCAGCCAAGGCCCCAAGGCTGATGCTCaaa                                                                                                                                          >  2:221215/1‑145 (MQ=255)            ccccTGGCGCGGCGGTTGGGTCTGGCGCACTTCGTCGCGCGAGTTGAGCGAGCGCCCCGGTCCACCCTCACCTTGCACCTGCTGGCCACGTTCACCGTGCCCCTCGCCGGGCCCTTCAGCCAAGGCCCCAAGGCTGATGCTCAAACCCa                                                                                                                                      <  4:56624/149‑1 (MQ=255)              ccTGGCGCGGCGGTTGGGTCTGGCGCACTTCGTCGCGCGAGTTGAGCGAGCGCCCCGGTCCACCCTCACCTTGCACCTGCTGGCCACGTTCACCGTGCCCCTCGCCGGGCCCTTCAGCCAAGGCCCCAAGGCTGATGCTCAAACCCAg                                                                                                                                     >  1:306229/1‑148 (MQ=255)                tGGCGCGGCGGTTGGGTCTGGCGCACTTCGTCGCGCGAGTTGAGCGAGCTCCCCGGTCCACCCTCACCTTGCACCTGCTGGCCACGTTCACCGTGCCCCTCGCCGGGCCCTTCAGCCAAGGCCCCAAGGCTGATGCTCAAACCCAGTa                                                                                                                                   >  1:90908/1‑148 (MQ=255)                             gggTCTGGCGCACTTCGTCGCGCGAGTTGAGCGAGCGCCCCGGTCCACCCTCACCTTGCACCTGCTGGCCACGTTCACCGTGCCCCTCGCCGGGCCCTTCAGCCAAGGCCCCAAGGCTGATGCTCAAACCCAGTAAACCAACACCTGc                                                                                                                      >  8:301256/1‑148 (MQ=255)                             gggTCTGGCGCACTTCGTCGCGCGAGTTGAGCGAGCGCCCCGGTCCACCCTCACCTTGCACCTGCTGGCCACGTTCACCGTGCCCCTCGCCGGGCCCTTCAGCCAAGGCCCCAAGGCTGATGCTCAAACCCAGTAAACCAACACCTGc                                                                                                                      <  7:308193/148‑1 (MQ=255)                               gTCTGGCGCACTTCGTCGCGCGAGTTGAGCGAGCGCCCCGGTCCACCCTCACCTTGCACCTGCTGGCCACGTTCACCGTGCCCCTCGCCGGGCCCTTCAGCCAAGGCCCCAAGGCTGATGCTCAAACCCAGTAAACCAACACCTGCCaa                                                                                                                   >  4:380163/1‑149 (MQ=255)                                    gcgcACTTCGTCGCGCGAGTTGAGCGAGCGCCCCGGTCCACCCTCACCTTGCACCTGCTGGCCACGTTCACCGTGCCCCTCGCCGGGCCCTTCATCCAACGCCCCAACGCTCATGCTCAAACCCAGTAAACCAACACCTGCCAACTGc                                                                                                               >  6:259548/1‑148 (MQ=255)                                    gcgcACTTCGTCGCGCGAGTTGAGCGAGCGCCCCGGTCCACCCTCACCTTGCACCTGCTGGCCACGTTCACCGTGCCCCTCGCCGGGCCCTTCAGCCAAGGCCCCAAGGCTGATGCTCAAACCCAGTAAACCAACACCTGCCAACTGc                                                                                                               >  5:250714/1‑148 (MQ=255)                                     cgcACTTCGTCGCGCGAGTTGAGCGAGCGCCCCGGTCCACCCTCACCTTGCACCTGCTGGCCACGTTCACCGTGCCCCTCGCCGGGCCCTTCAGCCAAGGCCCCAAGGCTGATGCTCAAACCCAGTAAACCAACACCTGCCAACTGCCa                                                                                                             <  2:20338/149‑1 (MQ=255)                                      gcACTTCGTCGCGCGAGTTGAGCGAGCGCCCCGGTCCACCCTCACCTTGCACCTGCTGGCCACGTTCACCGTGCCCCTCGCCGGGCCCTTCAGCCAAGGCCCCAAGGCTGATGCTCAAACCCAGTAAACCAACACCTGCCAACTGc                                                                                                               >  7:238447/1‑146 (MQ=255)                                             gtcgCGCGAGTTGAGCGAGCGCCCCGGTCCACCCTCACCTTGCACCTGCTGGCCACGTTCACCGTGCCCCTCGCCGGGCCCTTCAGCCAAGGCCCCAAGGCTGGTGCTCAAACCCAGTAAACCAACACCTGCCAACTGCCAGATGcgcg                                                                                                     >  5:144981/1‑149 (MQ=255)                                             gtcgCGCGAGTTGAGCGAGCGCCCCGGTCCACCCTCACCTTGCACCTGCTGGCCACGTTCACCGTGCCCCTCGCCGGGCCCTTCAGCCAAGGCCCCAAGGCTGGTGCTCAAACCCAGTAAACCAACACCTGCCAACTGCCAGATGcgcg                                                                                                     >  5:138806/1‑149 (MQ=255)                                                  gcgAGTTGAGCGAGCGCCCCGGTCCACCCTCACCTTGCACCTGCTGGCCACGTTCACCGTGCCCCTCGCCGGGCCCTTCAGCCAAGGCCCCAAGGCTGATGCTCAAACCCAGTAAACCAACACCTGCCAACTGCCAGATGCGCGACGt                                                                                                 >  8:24417/1‑148 (MQ=255)                                                     aGTTGAGCGAGCGCCCCGGTCCACCCTCACCTTGCACCTGCTGGCCACGTTCACCGTGCCCCTCGCCGGGCCCTTCAGCCAAGGCCCCAAGGCTGATGCTCAGACCCAGTAAACCAACACCTGGCACCTGCCAGATGCGCGACGTCATg                                                                                             >  8:258279/1‑149 (MQ=255)                                                     aGTTGAGCGAGCGCCCCGGTCCACCCTCACCTTGCACCTGCTGGCCACGTTCACCGTGCCCCTCGCCGGGCCCTTCAGCCAAGGCCCCAAGGCTGATGCTCAAACCCAGTAAACCAACACCTGCCAACTGCCAGATGCGCGACGTCATg                                                                                             <  3:184930/149‑1 (MQ=255)                                                           gcgagcgCCCCGGTCCACCCTCACCTTGCACCTGCTGGCCACGTTCACCGTGCCCCTCGCCGGGCCCTTCAGCCAAGGCCCCAAGGCTGATGCTCAAACCCAGTAAACCAACACCTGCCAACTGCCAGATGCGCGACGTCATGGTGTTc                                                                                       >  1:196083/1‑149 (MQ=255)                                                                  ccccGGTCCACCCTCACCTTGCACCTGCTGGCCACGTTCACCGTGCCCCTCGCCGGGCCCTTCAGCCAAGGCCCCAAGGCTGGTGCTCAAACCCAGTAAACCAACACCTGCCAACTGCCAGATGCGCGACGTCATGGTGTTCCTCTCgg                                                                                <  6:138806/149‑1 (MQ=255)                                                                  ccccGGTCCACCCTCACCTTGCACCTGCTGGCCACGTTCACCGTGCCCCTCGCCGGGCCCTTCAGCCAAGGCCCCAAGGCTGGTGCTCAAACCCAGTAAACCAACACCTGCCAACTGCCAGATGCGCGACGTCATGGTGTTCCTCTCgg                                                                                <  6:144981/149‑1 (MQ=255)                                                                   cccGGTCCACCCTCACCTTGCACCTGCTGGCCACGTTCACCGTGCCCCTCGCCGGGCCCTTCAGCCAAGGCCCCAAGGCTGATGCTCAAACCCAGTAAACCAACACCTGCCAACTGCCAGATGCGCGACGTCATGGTGTTCCTCTCgg                                                                                >  8:26165/1‑148 (MQ=255)                                                                   cccGGTCCACCCTCACCTTGCACCTGCTGGCCACGTTCACCGTGCCCCTCGCCGGGCCCTTCAGCCAAGGCCCCAACGCTGATGCTCAAACCCAGGAAACCCTCACCTGCCAACTCCCAGCTGCGCGACGTCAGTGTGTTCgtcggggt                                                                               >  8:127025/1‑144 (MQ=255)                                                                        gccACCCTCACCTTGCACCTGCTGGCCACGGTCACCGTGCCCCTCGCCGGGCCCTTCAGCCCAGGCCCCAAGGCTGATGCTCCAACCCAGGAACCCCCCCCCTGCCAACTGCCACCTGCGCGACGCCATCGTGCTCCTCTCGGCAGaaa                                                                          <  6:95824/148‑1 (MQ=255)                                                                           aCCCTCACCTTGCACCTGCTGGCCACGTTCACCGTGCCCCTCGCCGGGCCCTTCAGCCAAGGCCCCAAGGCTGATGCTCAAACCCGGTAAACCAACACCTGCCAACTGCCAGATGCGCGACGTCATGGTGTTCCTCTCGGTAGaaaaaa                                                                       <  4:202311/149‑1 (MQ=255)                                                                           aCCCTCACCTTGCACCTGCTGGCCACGTTCACCGTGCCCCTCGCCGGGCCCTTCAGCCAAGGCCCCAAGGCTGATGCTCAAACCCAGTAAACCAACACCTGCCAACTGCCAGATGCGCGACGTCATGGTGTTCCTCTCGGTAGaaaaaa                                                                       >  8:238783/1‑149 (MQ=255)                                                                             ccTCACCTTGCACCTGCTGGCCACGTTCACCGTGCCCCTCGCCGGGCCCTTCAGCCAAGGCCCCAAGGCTGATGCTCAAACCCAGTAAACCAACACCTGCCAACTGCCAGATgcg                                                                                                       >  4:292074/1‑115 (MQ=255)                                                                             ccTCACCTTGCACCTGCTGGCCACGTTCACCGTGCCCCTCGCCGGGCCCTTCAGCCAAGGCCCCAAGGCTGATGCTCAAACCCAGTAAACCAACACCTGCCAACTGCCAGATgcg                                                                                                       <  3:292074/115‑1 (MQ=255)                                                                             ccTCACCTTGCACCTGCTGGCCACGTTCACCGTGCCCCTCGCCGGGCCCTTCAGCCAAGGCCCCAAGGCTGATGCTCAAACCCAGTAAACCAACACCTGCCAACTGCCAGATGCGCGACGTCATGGTGTTCCTCTCGGTAGAAAAAACt                                                                     <  4:165083/149‑1 (MQ=255)                                                                               tCACCTTGCACCTGCTGGCCACGTTCACCGTGCCCCTCGCCGGGCCCTTCAGCCAAGGCCCCAAGGCTGATGCTCAAACCCAGTAAACCAACACCTGCCAACTGCCAGATGCGCGACGTCATGGTGTTCCTCTCGGTAGAAAAAACTg                                                                    >  2:141611/1‑148 (MQ=255)                                                                                 aCCTTGCACCTGCTGGCCACGTTCACCGTGCCCCTCGCCGGGCCCTTCAGCCAAGGCCCCAAGGCTGATGCTCAAACCCAGTAAACCAACACCTGCCAACTGCCAGATGCGCGACGTCATGGTGTTCCTCTCGGTAGAAAAAACTGc                                                                   <  8:275453/147‑1 (MQ=255)                                                                                  ccTTGCACGTGCTTGCCACGTTCACCTTGCCTCTCGCCGGGCACTTCAGCCAAGTCCCCAAGGCTGATGCTCAAACCCAGTAAACCAACACCTTCCAACTGCCAGATGCGCGACGTCATGGTGTTCCTCTCGGTAGAAAAAACTGcaac                                                                <  5:259548/149‑1 (MQ=255)                                                                                  ccTTGCACCTGCTGGCCACGTTCACCGTGCCCCTCGCCGGGCCCTTCAGCCAAGGCCCCAAGGCTGATGCTCAAACCCAGTAAACCAACACCTGCCAACTGCCAGATGCGCGACGTCATGGTGTTCCTCTCGGTAGAAAAAACTGcaac                                                                <  2:306229/149‑1 (MQ=255)                                                                                    ttGCACCTGCTGGCCACGTTCACCGTGCCCCTCGCCGGGCCCTTCAGCCAAGGCCCCAAGGCTGATGCTCAAACCCAGTAAACCAACACCTGCCAACTGCCAGATGCGCGACGTCATGGTGTTCCTCTCGGTAGAAAAAACTGcaacc                                                               >  5:167089/1‑148 (MQ=255)                                                                                       cACCTGCTGGCCACGTTCACCGTGCCCCTCGCCGGGCCCTTCAGCCAAGGCCCCAAGGCTGATGCTCAAACCCAGTAAACCAACACCTGCCAACTGCCAGATGCGCGACGTCATGGTGTCCCTCTCGGTAGAAAAAACTGCAACCAACg                                                           >  1:191834/1‑149 (MQ=255)                                                                                        cccTGCTGGCCACGTTCACCGTGCCCCTCGCCGGGCCCTTCAGCCAAGGCCCCAAGGCTGATGCTCAAACCCAGTAAACCAACACCTGCCAACTGCCAGATGCGCGACGTCATGGTGTTCCTCTCGGTAGAAAAAACTGCAACCAACGc                                                          <  1:141611/148‑1 (MQ=255)                                                                                        aCCTGCTGGCCACGTTCACCGTGCCCCTCGCCGGGCCCTTCAGCCAAGGCCCCAAGGCTGATGCTCAAACCCAGTAAACCAACACCTGCCAACTGCCAGATGCGCGCCGTCATGGTGTTCCTCTCGGTAGAAAAAACTGCAACCAACGc                                                          >  4:262832/1‑149 (MQ=255)                                                                                        aCCTGCTGGCCACGTTCACCGTGCCCCTCGCCGGGCCCTTCAGCCAAGGCCCCAAGGCTGATGCTCAAACCCAGTAAACCAACACCTGCCAACTGCCAGATGCGCGACGTCATGGTGTTCCTCTCGGTAGAAAAAACTGCAACCAACGc                                                          <  7:301256/149‑1 (MQ=255)                                                                                        aCCTGCTGGCCACGTTCACCGTGCCCCTCGCCGGGCCCTTCAGCCAAGGCCCCAAGGCTGATGCTCAAACCCAGTAAACCAACACCTGCCAACTGCCAGATGCGCGACGTCATGGTGTTCCTCTCGGTAGAAAAAACTGCAACCAACGc                                                          <  8:238447/149‑1 (MQ=255)                                                                                          ctgctgGCCACGTTCACCGTGCCCCTCGCCGGGCCCTTCAGCCAAGGACCCAAGGCTGATGCTCAAACCCAGTAAACCAACACCTGCCAACTGCCAGATGCGCGACGTCATGGTGTTCCTCTCTGTAGAAAAAACTGCAACCAACGCCa                                                        <  3:322108/149‑1 (MQ=255)                                                                                                gCCACGTTCACCGTGCCCCTCGCCGGGCCCTTCAGCCAAGGCCCCAAGGCTGATGCTCAAACCCAGTAAACCAACACCTGCCAACTGCCAGATGCGCGACGTCATGGTGTTCCTCTCGGTAGAAAAAACTGCAACCAACGCCAGTCTa                                                   <  8:124367/148‑1 (MQ=255)                                                                                                   aCGTTCACCGTGCCCCTCGCCGGGCCCTTCAGCCAAGGCCCCAAGGCTGATGCTCAAAACCAGTAAACCAACACCTGCCAACTGCCAGATGCGCTACGTCATGGTGTTCCTCTCGGTAGAAAAAACTGCAACCAACGCAAGTCTACCGt                                               >  2:117060/1‑149 (MQ=255)                                                                                                        cACCGTGCCCCTCGCCGGGCCCTTCAGCCAAGGCCCCAAGGCTGATGCTCAAACCCAGTAAACCAACACCTGCCAACTGCCAGATGCGCGACGTCATGGTGTTCCTCTCGGTAGAAAAAACTGCAACCAACGCCAGTCTACCGTGccc                                           <  7:197429/148‑1 (MQ=255)                                                                                                        cACCGTGCCCCTCGCCGGGCCCTTCAGCCAAGGCCCCAAGGCTGATGCTCAAACCCAGTAAACCAACACCTGCCAACTGCCAGATGCGCGACGTCATGGTGTTCCTCTCGGTAGAAAAAACTGCAACCAACGCCAGTCTACCGTGccc                                           <  2:90908/148‑1 (MQ=255)                                                                                                           cGTGCCCCTCGCCGGGCCCTTCAGCCAAGTCCCCTAGGCTGATGCTCAAACCCAGTAAACGAACACCGGCCAACTGCCAGATGCGCGACGTCATGGTGTTGCTCTCGGTTGAAAAAACTGCCACCAACGCCAGTCTACCGTGCccacca                                       >  5:367514/1‑149 (MQ=255)                                                                                                           cGTGCCCCTCGCCGGGCCCTTCAGCCAAGGCCCCAAGGCTGATGCTCAAACCCAGTAAACCAACACCTGCCAACTGCCAGATGCGCGACGTCATGGTGTTCCTCTCGGTAGAAAAAACt                                                                     >  2:220455/1‑119 (MQ=255)                                                                                                           cGTGCCCCTCGCCGGGCCCTTCAGCCAAGGCCCCAAGGCTGATGCTCAAACCCAGTAAACCAACACCTGCCAACTGCCAGATGCGCGACGACATGGTGTTCCTCTCGGTAGAAAAAACt                                                                     <  1:220455/119‑1 (MQ=255)                                                                                                               ccccTCGCCGGGCCCTTCAGCCAAGGCCCCAAGGCTGATGCTCAAACCCAGTAAACCAACACCTGCCAACTGCCAGATGCGCGACGTCATGGTGTTCCTCTCGGTAGAAAAAACTGCAACCAACGCCAGTCTACCGTGCCCACcacaca                                   <  6:250714/149‑1 (MQ=255)                                                                                                                cccTCGCCGGGCCCTTCAGCCAAGGCCCCAAGGCTGATGCTCAAACCCAGTAAACCAACACCTGCCAACTGCCAGATGCGCGACGTCATGGTGTTCCTCTCGGTAGAAAAAACTGCAACCAACGCCAGTCTACCGTGCCCACCACACAg                                  <  7:196578/149‑1 (MQ=255)                                                                                                                 ccTCGCCGGGCCCTTCAGCCAAGGCCCCAAGGCTGATGCTCAAACCCAGTAAACCAACACCTGCCAACTGCCAGATGCGCGACGTCATGGTGTTCCTCTCGGTAGAAAAAACTGCAACCAACGCCAGTCTACCGTGCCCACCACACAgg                                 >  6:307258/1‑149 (MQ=255)                                                                                                                   tCGCCGGGCCCTTCAGCCAAGGCCCCAAGGCTGATGCTCAAACCCAGTAAACCAACACCTGCCAACTGCCAGATGCGCGACGTCATGGTGTTCCTCTCGGTAGAAAAAACTGCAACCAACGCCAGTCTACCGTGCCCACCACACAGGc                                >  6:354059/1‑148 (MQ=255)                                                                                                                     gCCGGGCCCTTCAGCCAAGGCCCCAAGGCTGATGCTCAAACCCAGTAAACCAACACCTGCCAACTGCCAGATGCGCGACGTCATGGTGTTCCTCTCGGTAGAAAAAACTGCAACCAACGCCAGTCTACCGTGCCCACCACACAGGcc                               <  1:348665/147‑1 (MQ=255)                                                                                                                          gCCCTTCAGCCAAGGCCCCAAGGCTGATGCTCAAACCCAGTAAACCAACACCTGCCAACTCCCAGATGCGCGCAGTCATGGTGTTCCTATCGGTAGAAAAAACTGCAACCAACGCCAGTCTACCGTGCCCCCCACACAGGCCAATTaa                         >  7:347004/1‑148 (MQ=255)                                                                                                                              ttCAGCCAAGGCCCCAAGGCTGATGCTCAAACCCAGTAAACCAACACCTGCCAACTGCCAGATGCGCGACGTCATGGTGTTACTCTCGGTAGAAAAAACTGCAACCAACCCCAGTCTACCGTGCCCACCACCCAGGCCAATTAAAAtta                    >  3:11309/1‑148 (MQ=255)                                                                                                                              ttCAGCCAAGGCCCCAAGGCTGATGCTCAAACCCAATAAACCAACACCTGCCAACTGCCAGATGCGCGACGTCATGGTGTTCCTCTCGGTAGAAAAAACTGCAACCAACGCCAGTCTACCGTGCCCACCACACAGGCCAATTAAAATTc                    >  6:253452/1‑149 (MQ=255)                                                                                                                                cAGCCAAGGCCCCAAGGCTGATGCTCAAACCCAGTAAACCAACACCTGCCAACTGCCAGATGCGCGACGTCATGGTGTTCCTCTCGGTAGAAAAAACTGCAACCAACGCCAGTCTACCGTGCCCACCACACAGGCCAATTAAAATTc                    <  7:26165/147‑1 (MQ=255)                                                                                                                                  gCCAAGGCCCCAAGGCTGATGCTCAAACCCAGTAAACCAACACCTGCCAACTGCCAGATGCGCGACGTCATGGTGTTCCTCTCGGTAGAAAAAACTGCAACCAACGCCAGTCTACCGTGCCCACCACACAGGCCAATTAAAATTCCCGa                <  3:332509/149‑1 (MQ=255)                                                                                                                                   ccAAGGCCCCAAGGCTGATGCTCAAACCCAGTAAACCAACACCTGCCAACTGCCAGATGCGCGACGTCATGGTGTTCCTCTCGGTAGAAAAAACTGCAACCAACGCCAGTCTACCGTGCCCACCACACAGGCCAATTAAAATTCCCGa                >  8:201205/1‑148 (MQ=255)                                                                                                                                   ccAAGGCCCCAAGGCTGATGCTCAAACCCAGTAAACCAACACCTGCCAACTGCCAGATGCGCGACGTCATGGTGTTCCTCTCGGTAGAAAAAACTGCAACCAACGCCAGTCTACCGTGCCCACCACACAGGCCAATTAAAATTCCCGAc               <  6:365779/149‑1 (MQ=255)                                                                                                                                    cAAGGCCCCAAGGCTGATGCTCAAACCCAGTAAACCAACACCTGCCAACTGCCAGATGCGCGACGTCATGGTGTTCCTCTCGGTATAAAAAACTGCAACCAACGCCAGTCTACCGTGCCCACCACACAGGCCAATTAAAATTCCCGACa              <  7:238783/149‑1 (MQ=255)                                                                                                                                    cAAGGCCCCAAGGCTGATGCTCAAACCCAGTAAACCAACACCTGCCAACTGCCAGATGCGCGACGTCATGGTGTTCCTCTCGGTAGAAAAAACTGCAACCAACACGAGTCTACCGTGCCCACCACACAGGCCAATTAAAATTCCCGACa              >  1:328549/1‑149 (MQ=255)                                                                                                                                        gCCCCAAGGCTGATGCTCAAACCCAGTAAACCAACACCTGCCAACTGCCAGATGCGCGACGTCATGGTGTTCCTCTCGGTAGAAAAAACTGCAACCAACGCCAGTCTACCGTGCCCACCACACAGGCCAATTAAAATTCCCGACATGa           >  4:337349/1‑148 (MQ=255)                                                                                                                                             aaGGCTGATGCTCAAACCCAGTAAACCAACACCTGCCAACTGCCAGATGCGCGACGTCATGGTGTTCCTCTCGGTAGAAAAAACTGCAACCAACGCCAGTCTACCGTGCCCACCACACAGGCCAATTAAAATTCCCGACATGAAAAAAg     >  2:212640/1‑149 (MQ=255)                                                                                                                                              aGGCTGATGCTCAAACCCAGTAAACCAACACCTGCCAACTGCCAGATGCGCGACGTCATGGTGTTCCTCTCGGTAGAAAAAACTGCAACCAACGCCAGTCTACCGTGCCCACCACACAGGCCAATTAAAATTCCCGACATGAAAAAAgg    <  7:24417/149‑1 (MQ=255)                                                                                                                                                tcTGATGCTCAAACCCAGTAAACCACCACCTGCCAACTGCCAGATGCGCGACGTCATGGTGTTCCTCTCGGTAGAATACACTGCAACCAACGCCATTCTACCGTGCCCACCACACAGGCCAATTAAAATTCCCGACATGAAAAAAGGga  <  8:347004/148‑1 (MQ=255)                                                                                                                                                gCTGATGCTCAAACCCAGTAAACCAACACCTGCCACCTGCCAGATGCGCGACGTCATGGTGTTCCTCTCGGTAGAAAAAACTGCAACCACCGCCAGTCTACCGTGCCCACCACACAGGCCAATTAAAATTGCCGACATGAAAAAAgg    >  4:53618/1‑147 (MQ=255)                                                                                                                                                gCTGATGCTCAAACCCAGTAAACCAACACCTGCCAACTGCCAGATGCGCGACGTCATGGTGTCCCTCTCGGTAGAAAAAACTGCAACCAACGCCAGTCTACCGTGCCCACCACACAGGCCAATTAAAATTCCCGACATGAAAAAAGGga  <  2:191834/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                 TCCTGGTAATACCCCTGGCGCGGCGGTTGGGTCTGGCGCACTTCGTCGCGCGAGTTGAGCGAGCGCCCCGGTCCACCCTCACCTTGCACCTGCTGGCCACGTTCACCGTGCCCCTCGCCGGGCCCTTCAGCCAAGGCCCCAAGGCTGATGCTCAAACCCAGTAAACCAACACCTGCCAACTGCCAGATGCGCGACGTCATGGTGTTCCTCTCGGTAGAAAAAACTGCAACCAACGCCAGTCTACCGTGCCCACCACACAGGCCAATTAAAATTCCCGACATGAAAAAAGGGA  >  NC_002947/5437359‑5437650

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