Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	554,157	Δ1 bp	5.4%	coding (97/714 nt)	rpsC →	30S ribosomal protein S3

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	554,157	0	G	.	5.4%	242.8 / 6.6	74	coding (97/714 nt)	rpsC	30S ribosomal protein S3
Reads supporting (aligned to +/- strand):  ref base G (41/29);  new base . (2/2);  total (43/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.52e-01

AGGCCGTGCTGATCGCATCGTCAAGCGGTCTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGAT‑GGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTG‑AAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCGTAAGCCGTATCGATATTCATCGTCCGGCTCAAACTGCACGCATCACCATCCACACCGCTC  >  NC_002947/554011‑554301                                                                                                                                                     |                                                                                                                                                  aGGCCGTGCTGATCGCATCGTCAAGCGGTCTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGAT‑GGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCa                                                                                                                                                 <  4:179938/149‑1 (MQ=255) aGGCCGTGCTGATCGCATCGTCAAGCGGTCTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGAT‑GGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCa                                                                                                                                                 <  2:225479/149‑1 (MQ=255)    ccGTGCTGATCGCATCGTCAAGCGGTCTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGAT‑GGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCTGTCTG‑GTACGCAGAc                                                                                                                                              >  8:244246/1‑149 (MQ=255)      gTGCTGATCGCATCGTCAAGCGGTCTTGACATATCACTGTCAAGGTTGCGGACAAGTAACTGAGTCGATCAGAT‑GGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTGCGTAC‑CAGACgg                                                                                                                                            <  5:320297/149‑1 (MQ=255)         cTGATCGCATCGTCAAGCGGTCTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGAT‑GGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCTGTCTG‑GTACGCAGACGGTg                                                                                                                                          <  7:244246/148‑1 (MQ=255)         cTGATCGCATCGTCAAGCGGTCTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGAT‑GGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTg                                                                                                                                          <  8:117438/148‑1 (MQ=255)           gATCGCATCGTCAAGCGGTCTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGAT‑GGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTa                                                                                                                                       >  3:265742/1‑149 (MQ=255)            aTCGCATCGTCAAGCGGTCTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGAT‑GGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTAc                                                                                                                                      <  8:254373/149‑1 (MQ=255)            aTCGCATCGTCAAGCAGTCTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGAT‑GGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCACCGTCTGCGTAC‑CAGACggaggt                                                                                                                                        <  2:25510/147‑5 (MQ=255)             tCGCATCGTCAAGCGGTCTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGAT‑GGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACt                                                                                                                                     >  4:306760/1‑149 (MQ=255)               gCATCGTCAAGCGGTCTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGAT‑GGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACtt                                                                                                                                    <  2:253522/148‑1 (MQ=255)               gCATCGTCAAGCGGTCTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGAT‑GGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACtt                                                                                                                                    <  1:292943/148‑1 (MQ=255)               gCATCGTCAAGCGGTCTTGCCATATCACTGTCAAGGTAGCGGACAAGTAACGGAGTCGATCAGAT‑GGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACTCCTCCGTCTG‑GTACGCAGACGGTGCTACtt                                                                                                                                    <  3:245266/148‑1 (MQ=255)                     tCAAGCGGTCTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGAT‑GGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACg                                                                                                                                 >  5:239896/1‑145 (MQ=255)                         gCGGTCTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGAT‑GGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTAc                                                                                                                         <  6:232804/149‑1 (MQ=255)                            gTCTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCt                                                                                                                       >  1:248909/1‑149 (MQ=255)                              cTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGAT‑GGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGtt                                                                                                                    >  4:196711/1‑149 (MQ=255)                               ttGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGAT‑GGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGtt                                                                                                                    >  1:321113/1‑148 (MQ=255)                               ttGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGAT‑GGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGtt                                                                                                                    >  1:316190/1‑148 (MQ=255)                                    atatCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGAT‑GGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTG‑AAGGa                                                                                                             >  6:162441/1‑149 (MQ=255)                                        cACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGAT‑GGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTG‑AAGGATCt                                                                                                          >  8:177680/1‑148 (MQ=255)                                           tGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGAT‑GGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTG‑AAGGATCTGaa                                                                                                       >  6:210731/1‑148 (MQ=255)                                             tCAAGGTTGCGGACAAGTAACGGAGTCGATCAGAT‑GGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTG‑AAGGATCTGaaaa                                                                                                     <  6:330538/148‑1 (MQ=255)                                               aaGGTTGCGGACAAGTAACGGAGTCGATCAGAT‑GGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTG‑AAGGATCTGAAAAcg                                                                                                   <  7:65019/148‑1 (MQ=255)                                                 ggTTGCGGACAAGTAACGGAGTCGATCAGAT‑GGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTG‑AAGGATCTGAAAACGCGt                                                                                                <  4:265742/149‑1 (MQ=255)                                                 ggTTGCGGACAAGTAACGGAGTCGATCAGAT‑GGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTG‑AAGGATCTGAAAACGCGt                                                                                                <  7:177680/149‑1 (MQ=255)                                                   ttGCGGACAAGTAACGGAGTCGATCAGAT‑GGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTG‑AAGGATCTGAAAACGCGTGa                                                                                              <  5:210731/149‑1 (MQ=255)                                                       ggACAAGTAACGGAGTCGATCAGAT‑GGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAACCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTG‑AAGGATCTGAAAACGCGTGAGTAc                                                                                          >  6:206118/1‑149 (MQ=255)                                                             gTAACGGAGTCGATCAGAT‑GGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTG‑AAGGATCTGAAAACGCGTGAGTACCTcc                                                                                      >  4:198293/1‑147 (MQ=255)                                                             gTAACGGAGTCGATCAGAT‑GGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTG‑AAGGATCTGAAAACGCGTGAGTACCTcc                                                                                      >  1:166019/1‑147 (MQ=255)                                                                aCGGAGTCGATCAGAT‑GGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTG‑AAGGATCTGAAAACGCGTGAGTACCTCCAAGa                                                                                  >  1:314333/1‑148 (MQ=255)                                                                    aGTCGATCAGAT‑GGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTG‑AAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAAc                                                                             <  2:246817/149‑1 (MQ=255)                                                                         aTCAGAT‑GGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTG‑AAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTaaaa                                                                        <  2:248909/149‑1 (MQ=255)                                                                              aT‑GGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTG‑AAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAgcgc                                                                   <  5:162441/149‑1 (MQ=255)                                                                               aaggctCAGAAAGAACATCCCACTGGCCTTCGCCTGGGAATCCTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTG‑AAGGATCTGAAAACGCGGGAGTACCTCCAAGACAAACTAAAAAgcgc                                                                   <  4:52534/144‑1 (MQ=255)                                                                                  ggTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTG‑AAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTc                                                                <  2:173819/149‑1 (MQ=255)                                                                                   gTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTG‑AAGGATCTGAAAACGCGTGAGTAc                                                                                          <  6:155267/122‑1 (MQ=255)                                                                                   gTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTG‑AAGGATCTGAAAACGCGTGAGTAc                                                                                          >  5:155267/1‑122 (MQ=255)                                                                                    tCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTG‑AAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCg                                                              >  3:175026/1‑149 (MQ=255)                                                                                    tCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTG‑AAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCg                                                              >  3:60286/1‑149 (MQ=255)                                                                                    tCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTG‑AAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCg                                                              >  2:109596/1‑149 (MQ=255)                                                                                      aGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTG‑AAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCGt                                                             >  4:206781/1‑148 (MQ=255)                                                                                        aaaGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTG‑AAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCGTAAg                                                          >  2:327048/1‑149 (MQ=255)                                                                                          aGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTG‑AAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCGTAAGcc                                                        <  3:306760/149‑1 (MQ=255)                                                                                               aTCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTG‑AAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCGTAAGCCGTATc                                                   >  5:282582/1‑149 (MQ=255)                                                                                                         cATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTG‑AAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCGTAAGCCGTATCGATATTCATc                                         <  4:60286/149‑1 (MQ=255)                                                                                                               ccTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTG‑AAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCGTAAGCCGTATCGATATTCATCGTCCgg                                   >  5:321929/1‑149 (MQ=255)                                                                                                               ccTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTG‑AAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCGTAAGCCGTATCGATATTCATCGTCCgg                                   >  5:208624/1‑149 (MQ=255)                                                                                                               ccTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTG‑AAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCGTAAGCCGTATCGATATTCATCGTCCgg                                   <  2:314333/149‑1 (MQ=255)                                                                                                                 tGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTG‑AAGGATCTGAAAACGCGTGAGTACCTCCAAGACaaa                                                                              <  1:187401/104‑1 (MQ=255)                                                                                                                 tGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTG‑AAGGATCTGAAAACGCGTGAGTACCTCCAAGACaaa                                                                              >  2:187401/1‑104 (MQ=255)                                                                                                                      tctgTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTG‑AAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCGTAAGCCGTATCGATATTCATCGTCCGGCTCaaa                             <  7:81358/145‑1 (MQ=255)                                                                                                                      aTCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTG‑AAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCGTAAGCCGTATCGATATTCATCGTCCGGCTCAAAc                            >  1:197046/1‑149 (MQ=255)                                                                                                                         gTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTG‑AAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCGTAAGCCGTATCGATATTCATCGTCCGGCTCAAACTg                          >  1:288585/1‑148 (MQ=255)                                                                                                                          tCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTG‑AAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCGTAAGCCGTATCGATATTCATCGTCCGGCTCAAACTGCa                        >  4:99875/1‑149 (MQ=255)                                                                                                                          tCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACATACGCAGATTACCTGTTG‑AAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCTTCAGTAACCCGTATCGATATTAATCGTCCGGCTCAAAcaga                         >  8:81358/1‑145 (MQ=255)                                                                                                                           cAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTG‑AAAGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCGTAAGCCGTATCGATATTCATCGTCCGGCTCAAACTGCAc                       >  3:24286/1‑149 (MQ=255)                                                                                                                            aaGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTG‑AAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCGTAAGCCGTATCGATATTCATCGTCCGGCTCAAACTGCACg                      >  2:256004/1‑149 (MQ=255)                                                                                                                                 gCACTCCTCCGTCTG‑GTACGAAGACGGTGCTACTTACGCAGATTACCTGTTG‑AAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCGTAAGCCGTATCGATATTCATCGTCCGGCTCAAACTGCACGCa                    <  4:299090/146‑1 (MQ=255)                                                                                                                                 gCACAACTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTG‑AAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCGTAAGCCGTATCGATATTCATCGTCAGGCTCAAACTGCACGCa                    >  3:299090/1‑146 (MQ=255)                                                                                                                                  cacaCCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCGTAAGCCGTATCGATATTCATCGTCCGGCTCAAACTGCACGCATCa                 >  1:248319/1‑149 (MQ=255)                                                                                                                                  cacaCCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTG‑AAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCGTAAGCCGTATCGATATTCATCGTCCGGCTCAAACTGCACGCATCa                 >  5:155078/1‑148 (MQ=255)                                                                                                                                  cacaCCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTG‑AAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCGTAAGCCGTATCGATATTCATCGTCCGGCTCAAACTGCACGCATCAc                <  2:166019/149‑1 (MQ=255)                                                                                                                                     acCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTG‑AAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCGTAAGCCGTATCGATATTCATCGTCCGGCTCAAACTGCACGCATCACCAt             >  7:276877/1‑149 (MQ=255)                                                                                                                                     acCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTG‑AAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCGTAAGCCGTATCGATATTCATCGTCCGGCTCAAACTGCACGCATCACCAt             >  8:153132/1‑149 (MQ=255)                                                                                                                                      gcaccGTCTGCGTAC‑CAGACGGTGCTACTTACGCAGATTACCTGTTG‑AAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCGTAAGCCGTATCGATATTCATCGTCCGGCTCAAACTGCACGCATCACCATc            >  1:41460/4‑149 (MQ=255)                                                                                                                                       cTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTG‑AAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCGTAAGCCGTATCGATATTCATCGTCCGGCTCAAACTGCACGCATCACCATc            >  6:91063/1‑148 (MQ=255)                                                                                                                                         ccGTCTGCGTAC‑CAGACGGTGCTACTTACGCAGATTACCTGTTG‑AAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCGTAAGCCGTATCGATATTCATCGTCCGGCTCAAACTGCACGCATCACCATCcac         >  8:243050/1‑149 (MQ=255)                                                                                                                                           gTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTG‑AAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCGTAAGCCGTATCGATATTCATCGTCCGGCTCAAACTGCACGCATCACCATCCacac       >  7:27125/1‑149 (MQ=255)                                                                                                                                                 gTACGCAGACGGTGCTACTTACGCAGATTACCTGTTG‑AAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCGTAAGCCGTATCGATATTCATCGTCCGGCTCAAACTGCACGCATCACCATCCACACCGCTc  <  1:291106/149‑1 (MQ=255)                                                                                                                                                 gTACGCAGACGGTGCTACTTACGCAGATTACCTGTTG‑AAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCGTAAGCCGTATCGATATTCATCGTCCGGCTCAAACTGCACGCATCACCATCCACACCGCTc  >  2:136516/1‑149 (MQ=255)                                                                                                                                                 gTACGCAGACGGTGCTACTTACGCAGATTACCTGTTG‑AAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCGTAAGCCGTATCGATATTCATCGTCCGGCTCAAACTGCACGCATCACCATCCACACCGCTc  >  2:291106/1‑149 (MQ=255)                                                                                                                                                    cGCAGACGGTGCTACTTACGCAGATTACCTGTTG‑AAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCGTAAGCCGTATCGATATCCATCGTCCGGCTCAAACTGCACGCATCa                 <  7:291750/131‑1 (MQ=255)                                                                                                                                                    cGCAGACGGTGCTACTTACGCAGATTACCTGTTG‑AAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCGTAAGCCGTATCGATATCCATCGTCCGGCTCAAACTGCACGCATCa                 >  8:291750/1‑131 (MQ=255)                                                                                                                                                     |                                                                                                                                                  AGGCCGTGCTGATCGCATCGTCAAGCGGTCTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGAT‑GGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTG‑AAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCGTAAGCCGTATCGATATTCATCGTCCGGCTCAAACTGCACGCATCACCATCCACACCGCTC  >  NC_002947/554011‑554301

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