Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	554,152:1	+C	5.8%	coding (92/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,152	1	.	C	5.8%	276.6 / 10.1	69	coding (92/714 nt)	rpsC	30S ribosomal protein S3
Reads supporting (aligned to +/- strand):  ref base . (38/27);  new base C (2/2);  total (40/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.90e-01

GTGCCAAAGGCCGTGCTGATCGCATCGTCAAGCGGTCTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGAT‑GGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTG‑AAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCGTAAGCCGTATCGATATTCATCGTCCGGCTCAAACTGCACGCATCACCATCCACAC  >  NC_002947/554004‑554296                                                                                                                                                       |                                                                                                                                                  gTGCCAAAGGCCGTGCTGATCGCATCGTCAAGCGGTCTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGAT‑GGTTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTg                                                                                                                                                    <  4:319086/149‑1 (MQ=255)        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:60286/1‑149 (MQ=255)                                                                                           tCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTG‑AAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCg                                                         >  2:109596/1‑149 (MQ=255)                                                                                           tCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTG‑AAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCg                                                         >  3:175026/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:208624/1‑149 (MQ=255)                                                                                                                      ccTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTG‑AAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCGTAAGCCGTATCGATATTCATCGTCCgg                              >  5:321929/1‑149 (MQ=255)                                                                                                                      ccTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTG‑AAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCGTAAGCCGTATCGATATTCATCGTCCgg                              <  2:314333/149‑1 (MQ=255)                                                                                                                        tGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTG‑AAGGATCTGAAAACGCGTGAGTACCTCCAAGACaaa                                                                         >  2:187401/1‑104 (MQ=255)                                                                                                                        tGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTG‑AAGGATCTGAAAACGCGTGAGTACCTCCAAGACaaa                                                                         <  1:187401/104‑1 (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)                                                                                                                                                       |                                                                                                                                                  GTGCCAAAGGCCGTGCTGATCGCATCGTCAAGCGGTCTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGAT‑GGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTG‑AAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCGTAAGCCGTATCGATATTCATCGTCCGGCTCAAACTGCACGCATCACCATCCACAC  >  NC_002947/554004‑554296

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