Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,381,772	C→T	81.8%	E340K (GAG→AAG)	PP_3857 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,381,772	0	C	T	81.8%	140.6 / 23.0	66	E340K (GAG→AAG)	PP_3857	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (6/6);  new base T (28/26);  total (34/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGAGCCTTGCTGAATACTTCGCCGAAAACGCCTTCTTCAAACAGCGGCAGGCCATTGGCGCCGATCATAAAGCCACGACTGCCCAGTGACAGCGATTGACCCGTGCCTTGACGATGCCACAGCACCGGAGCTTCAGCAATCAGCTCGCCCGACTTCAACTTGGCATGCGGGGTCAGGGCGCCA‑CCAGCACCCGTCGTACGGTCGCTCAAGAAGCGCACGGTACTGCCTGTCGCCGCGGGCGCTGAGTCGTTAAACCCATCGCCACCAGGGAACAGGTCGCGCCAC  >  NC_002947/4381627‑4381911                                                                                                                                                  |                                                                                                                                             ggAGCCTTGCTGAATACTTCGCCGAAAACGCCTTCTTCAAACAGCGGCAGGCCATTGGCGCCGATCATAAAGCCACGACTGCCCAGTGACAGCGATTGACCCGTGCCTTGACGATGCCACAGCACCGGAGCTTCAGCAATCAGCTTGcc                                                                                                                                           >  4:275330/1‑149 (MQ=255)    gCCTTGCTGAATACTTCGCCGAAAACGCCTTCTTCAAACAGCGGCAGGCCATTGGCGCCGATCATAAAGCCACGACTGCCCAGTGACAGCGATTGACCCGTGCCTTGACGATGCCACAGCACCGGAGCTTCAGCAATCAGCTCGCCCGa                                                                                                                                        <  4:74502/149‑1 (MQ=255)      cTTGCTGAATACTTCGCCGAAAACGCCTTCTTCAAACAGCGGCAGGCCATTGGCGCCGATCATAAAGCCACGACTGCCCAGTGACAGCGATTGACCCGTGCCTTGACGATGCCACAGCACCGGAGCTTCAGCAATCAGCTTGCCCGACt                                                                                                                                      >  5:230806/1‑149 (MQ=255)       ttGCTGAATACTTCGCCGAAAACGCCTTCTTCAAACAGCGGCAGGCCATTGGCGCCGATCATAAAGCCACGACTGCCCAGTGACAGCGATTGACCCGTGCCTTGACGATGCCACAGCACCGGAGCTTCAGCAATCAGCTTGCCCGACtt                                                                                                                                     >  7:281503/1‑149 (MQ=255)          cTGAATACTTCGCCGAAAACGCCTTCTTCAAACAGCGGCAGGCCATTGGCGCCGATCATAAAGCCACGACTGCCCAGTGACAGCGATTGACCCGTGCCTTGACGATGCCACAGCACCGGAGCTTCAGCAATCAGCTTGCCCGACTTCaa                                                                                                                                  >  2:340966/1‑149 (MQ=255)                    cGCCGAAAACGCCTTCTTCAAACAGCGGCAGGCCATTGGCGCCGATCATAAAGCCACGACTGCCCAGTGACAGCGATTGACCCGTGCCTTGACGATGCCACAGCACCGGAGCTTCAGCAATCAGCTTGCCCGACTTCAACTTGGCATGc                                                                                                                        >  2:212119/1‑149 (MQ=255)                      ccGAAAACGCCTTCTTCAAACAGCGGCAGGCCATTGGCGCCGATCATAAAGCCACGACTGCCCAGTGACAGCGATTGACCCGTGCCTTGACGATGCCACAGCACCGGAGCTTCAGCAATCAGCTTGCCCGACTTCAACTTGGCATGCgg                                                                                                                      <  6:230806/149‑1 (MQ=255)                       cGAAAACGCCTTCTTCAAACAGCGGCAGGCCATTGGCGCCGATCATAAAGCCACGACTGCCCAGTGACAGCGATTGACCCGTGCCTTGACGATGCCACAGCACCGGAGCTTCAGCAATCAGCTTGCCCGACTTCAACTTGGCATGCggg                                                                                                                     <  1:212119/149‑1 (MQ=255)                              gCCTTCTTCAAACAGCGGCAGGCCATTGGCGCCGATCATAAAGCCACGACTGCCCAGTGACAGCGATTGACCCGTGCCTTGACGATGCCACAGCACCGGAGCTTCAGCAATCAGCTCGCCCGACTTCAACTTGGCATGCGGGGTCAggg                                                                                                              <  2:112722/149‑1 (MQ=255)                                cttcttCAAACAGCGGCAGGCCATTGGCTCCGATCATAAAGCCACGACTGCCCAGTGACAGCGATTGACCCGTGCCTTGACGATGCCACAGCACCGGAGCTTCAGCAATCAGCTTGCCCGACTTCAACTTGGCATGCGGGGTCAGGgcg                                                                                                            >  1:176630/1‑149 (MQ=255)                                     tcAAACAGCGGCAGGCCATTGGCGCCGATCATAAAGCCACGACTGCCCAGTGACAGCGATTGACCCGTGCCTTGACGATGCCACAGCACCGGAGCTTCAGCAATCAGCTTGCCCGACTTCAACTTGGCATGCGGGGTCAGGGCGcca‑c                                                                                                       >  1:193037/1‑148 (MQ=255)                                      cAAACAGCGGCAGGCCATTGGCGCCGATCATAAAGCCACGACTGCACAGTGACAGCGAATGACCCGTGCCTTGACGATGCCACAGCACCGGAGCTTCAGCAATCAGCTTGCCCGACTTCAACTTGCCATGCGCGGTCAGGGCGcca‑cca                                                                                                     >  1:138223/1‑149 (MQ=255)                                       aaaCAGCGGCAGGCCATTGGCGCCGATCATAAAGCCACGACTGCCCAGTGACAGCGGTTGACCCGTGCCTTGACGATGCCACAGCACCGGAGCTTCAGCAATCAGCTTGCCCGACTTCAACTTGGCATGCGGGGTCAGGGCGCCA‑CCAg                                                                                                    <  8:130058/149‑1 (MQ=255)                                              ggCAGGCCATTGGCGCCGATCATAAAGCCACGACTGCCCAGTGACAGCGATTGACCCGTGCCTTGACGATGCCACAGCACCGGAGCTTCAGCAATCAGCTTGCCCGACTTCAACTTGGCATGCGGGGTCAGGGCGCCA‑CCAGCACCcgt                                                                                             <  8:281503/149‑1 (MQ=255)                                                 aGGCCATTGGCGCCGATCATAAAGCCACGACTGCCCAGTGACAGCGATTGACCCGTGCCTTGACGATGCCACAGCACCGGAGCTTCAGCAATCAGCTTGCCCGACTTCAACTTGGCATGCGGGGTCAGGGCGCCA‑CCAGCACCcgtcgt                                                                                          <  8:194401/149‑1 (MQ=255)                                                 aGGCCATTGGCGCCGATCATAAAGCCACGACTGCCCAGTGACAGCGATTGACCCGTGCCTTGACGATGCCACAGCACCGGAGCTTCAGCAATCAGCTTGCCCGACTTCAACTTGGCATGCGGGGTCAGGGCGCCA‑CCAGCACCcgtcgt                                                                                          <  8:189144/149‑1 (MQ=255)                                                  ggCCATTGGCGCCGATCATAAAGCCACGACTGCCCAGTGACAGCGATTGACCCGTGCCTTGACGATGCCACAGCACCGGAGCTTCAGCAATCAGCTTGCCCGACTTCAACTTGGCATGCGGGGTCAGGGCGCCA‑CCAGCACCCGTCGTa                                                                                         <  4:329340/149‑1 (MQ=255)                                                    ccATTGGCGCCGATCATAAAGCCACGACTGCCCAGTGACAGCGATTGACCCGTGCCTTGACGATGCCACAGCACCGGAGCTTCAGCAATCAGCTTGCCCGACTTCAACTTGGCATGCGGGGTCAGGGCGCCA‑CCAGCACCCGTCGTACg                                                                                       <  2:193037/149‑1 (MQ=255)                                                     cATTGGCGCCGATCATAAAGCCACGACTGCCCAGTGACAGCGATTGACCCGTGCCTTGACGATGCCACAGCACCGGAGCTTCAGCAATCAGCTTGCCCGACTTCAACTTGGCATGCGGGGTCAGGGCGCCA‑CCAGCACCCGTCGTACgg                                                                                      <  3:275330/149‑1 (MQ=255)                                                                aTCATAAAGCCACGACTGCCCAGTGACAGCGATTGACCCGTGCCTTGACGATGCCACAGCACCGGAGCTTCAGCAATCAGCTTGCCCGACTTCAACTTGGCATGCGGGGTCAGGGCGCCA‑CCAGCACCCGTCGTACGGTCGCTCaagaa                                                                           <  2:176630/149‑1 (MQ=255)                                                                  cATAAAGCCACGACTGCCCAGTGACAGCGATTGACCCGTGCCTTGACGATGCCACAGCACCGGAGCTTCAGCAATCAGCTTGCCCGACTTCAACTTGGCATGCGGGGTCAGGGCGCCA‑CCAGCACCCGTCGTACGGTCGCTCAAGAAgc                                                                         >  2:292817/1‑149 (MQ=255)                                                                   aTAAAGCCACGACTGCCCAGTGACAGCGATTGACCCGTGCCTTGACGATGCCACAGCACCGGAGCTTCAGCAATCAGCTTGCCCGACTTCAACTTGGCATGCGGGGTCAGGGCGCCA‑CCAGCACCCGTCGTACGGTCGCTCAAGAAgcg                                                                        >  5:197996/1‑149 (MQ=255)                                                                    tAAAGCCACGACTGCCCAGTGACAGCGATTGACCCGTGCCTTGACGATGCCACAGCACCGGAGCTTCAGCAATCAGCTTGCCCGACTTCAACTTGGCATGCGGGGTCAGGGCGCCA‑CCAGCACCCGTCGTACGGTCGCTCAAGAAgcg                                                                        >  8:75026/1‑148 (MQ=255)                                                                    tAAAGCCACGACTGCCCAGTGACAGCGATTGACCCGTGCCTTGACGATGCCACAGCACCGGAGCTTCAGCAATCAGCTCGCCCGACTTCAACTTGGCATGCGGGGTCAGGGCGCCA‑CCAGCACCCGTCGTACGGTCGCTCAAGAAgcg                                                                        >  8:71569/1‑148 (MQ=255)                                                                     aaaGCCACGACTGCCCAGTGACAGCGATTGACCCGTGCCTTGACGATGCCACAGCACCGGAGCTTCAGCAATCAGCTTGCCCGACTTCAACTTGGCATGCGGGGTCAGGGCGCCA‑CCAGCACCCGTCGTACGGTCGCTCaagaag                                                                          <  2:38221/145‑1 (MQ=255)                                                                     aaaGCCACGACTGCCCAGTGACAGCGATTGACCCGTGCCTTGACGATGCCACAGCACCGGAGCTTCAGCAATCAGCTTGCCCGACTTCAACTTGGCATGCGGGGTCAGGGCGCCA‑CCAGCACCCGTCGTACGGTCGCTCaagaag                                                                          >  1:38221/1‑145 (MQ=255)                                                                         ccACGACTGCCCAGTGACAGCGATTGACCCGTGCCTTGACGATGCCACAGCACCGGAGCTTCAGCAATCAGCTCGCCCGAATTCAACTTGGCATGCGCGGTCAGGGCGCCA‑CCAGCACCCCTCGTACGGTCCCTCAAGAAGCGCACgg                                                                   >  8:312938/1‑148 (MQ=255)                                                                          cACGACTGCCCAGTGACAGCGATTGACCCGTGCCTTGACGATGCCACAGCACCGGAGCTTCAGCAATCAGCTTGCCCGACTTCAACTTGGCATGCGGGGTCAGGGCGcca‑cca                                                                                                     <  1:145498/113‑1 (MQ=255)                                                                          cACGACTGCCCAGTGACAGCGATTGACCCGTGCCTTGACGATGCCACAGCACCGGAGCTTCAGCAATCAGCTTGCCCGACTTCAACTTGGCATGCGGGGTCAGGGCGcca‑cca                                                                                                     >  2:145498/1‑113 (MQ=255)                                                                             gACTGCCCAGTGACAGCGATTGACCCGTGCCTTGACGATGCCACAGCACCGGAGCTTCAGCAATCAGCTTGCCCGACTTCAACTTGGCATGCGGGGTCAGGGCGCCA‑CCAGCACCCGTCGTACGGTCGCTCAAGAAGCGCACGGTACTg                                                              >  7:88051/1‑149 (MQ=255)                                                                                   ccAGTGACAGCGATTGACCCGTGCCTTGACGATGCCACAGCACCGGAGCTTCAGCAATCAGCTTGCCCGACTTCAACTTGGCATGCGGGGTCAGGGCGCCA‑CCAGCACCCGTCGTACGGTCGCACAAGAAGCGCACGGTACTGCCTGTc                                                        <  1:340966/149‑1 (MQ=255)                                                                                       tGACAGCGATTGACCCGTGCCTTGACGATGCCACAGCACCGGAGCTTCAGCAATCAGCTCGCCCGACTTCAACTTGGCATGCGGGGTCAGGGCGCCA‑CCAGCACCCGTCGTACGGTCGCTc                                                                                >  5:243870/1‑121 (MQ=255)                                                                                       tGACAGCGATTGACCCGTGCCTTGACGATGCCACAGCACCGGAGCTTCAGCAATCAGCTCGCCCGACTTCAACTTGGCATGCGGGGTCAGGGCGCCA‑CCAGCACCCGTCGTACGGTCGCTc                                                                                <  6:243870/121‑1 (MQ=255)                                                                                           aGCGATTGACCCGTGCCTTGACGATGCCACAGCACCGGAGCTTCAGCAATCAGCTTGCCCGACTTCAACTTGGCATGCGGGGTCAGGGCGCCA‑CCAGCACCCGTCGTACGGTCGCTCAAGAAGCGCACGGTACTGCCTGTCGCCGCggg                                                >  5:62137/1‑149 (MQ=255)                                                                                           aGCGATTGACCCGTGCCTTGACGATGCCACAGCACCGGAGCTTCAGCAATCAGCTTGCCCGACTTCAACTTGGCATGCGGGGTCAGGGCGCCA‑CCAGCACCCGTCGTACGGTCGCTCAAGAAGCGCACGGTACTGCCTGTCGCCGCggg                                                >  2:313812/1‑149 (MQ=255)                                                                                            gCGATTGACCCGTGCCTTGACGATGCCACAGCACCGGAGCTTCAGCAATCAGCTTGCCCGACTTCAACTTGGCATGCGGGGTCAGGGCGCCA‑CCAGCACCCGTCGTACGGTCGCTCAAGAAGCGCACGGTACTGCCTGTCGCCGCGGgc                                               <  1:292817/149‑1 (MQ=255)                                                                                            gCGATTGACCCGTGCCTTGACGATGCCACAGCACCGGAGCTTCAGCAATCAGCTCGCCCGACTTCAACTTGGCATGCGGGGTCAGGGCGCCA‑CCAGCACCCGTCGTACGGTCGCTCAAGAAGCGCACGGTACTGCCTGTCGCCGCGGgc                                               <  7:71569/149‑1 (MQ=255)                                                                                             cGATTGACCCGTGCCTTGACGATGCCACAGCACCGGAGCTTCAGCAATCAGCTTGCCCGACTTCAACTTGGCATGCGGGGTCAGGGCGCCA‑CCAGCACCCGTCGTACGGTCACTCAAGAAGCGCACGGTACTGCCTGTCGCCGCGGgcg                                              >  4:122814/1‑149 (MQ=255)                                                                                               aTTGACCCGTGCCTTGACGATGCCACAGCACCGGAGCTTCAGCAATCAGCTTGCCCGACTTCAACTTGGCATGCGGGGTCAGGGCGCCA‑CCAGCACCCGTCGTACGGTCGCTCaagaa                                                                           >  5:195176/1‑118 (MQ=255)                                                                                               aTTGACCCGTGCCTTGACGATGCCACAGCACCGGAGCTTCAGCAATCAGCTTGCCCGACTTCAACTTGGCATGCGGGGTCAGGGCGCCA‑CCAGCACCCGTCGTACGGTCGCTCaagaa                                                                           <  6:195176/118‑1 (MQ=255)                                                                                                ttGACCCGTGCCTTGACGATGCCACAGCACCGGAGCTTCAGCAATCAGCTTGCCCGACTTCAACTTGGCATGCGGGGTCAGGGCGCCA‑CCAGCACCCGTCGTACGGTCGCTCAAGAAGCGCACGGTACTGCCTGTCGCCGCGGGCGCTg                                           <  6:62137/149‑1 (MQ=255)                                                                                                     ccGTGCCTTGCCGATGCCACAGCACCGGAGCTTCAGCAATCAGCTTGCCCGACTTCAACTTGGCATGCGGGGTCAGGGAGCCA‑CCAGCACCCGTCGTACGGTCGCTCAAGAAGCGCACGGTACTGCCTGTCGCCGCGGGCGCTGAGTCg                                      >  3:286463/1‑149 (MQ=255)                                                                                                     ccGTGCCTTGACGATGCCACAGCACCGGAGCTTCAGCAATCAGCTTGCCCGACTTCAACTTGGCATGCGGGGTCAGGGCGCCA‑CCAGCACCCGTCGTACGGTCGCTCAAGAAGCGCACGGTACTGCCTGTCGCCGCGGGCGCTGAGTCg                                      >  5:108821/1‑149 (MQ=255)                                                                                                     ccGTGCCTTGACGATGCCACAGCACCGGAGCTTCAGCAATCAGCTCGCCCGACTTCAACTTGGCATGCGGGGTCAGGGCGCCA‑CCAGCACCCGTCGTACGGTCGCTCAAGAAGCGCACGGTACTGCCTGTCGCCGCGGGCGCTGAGTCg                                      >  6:28623/1‑149 (MQ=255)                                                                                                      cGTGCCTTGACGATGCCACAGCACCGGAGCTTCAGCAATCAGCTTGCCCGACTTCAACTTGGCATGCGGGGTCAGGGCGCCA‑CCAGCACCCGTCGTACGGTCGCTCAAGAAGCGCACGGTACTGCCTGTCGCCGCGGGCGCTGAGTCGt                                     <  6:197996/149‑1 (MQ=255)                                                                                                       gTGCCTTGACGATGCCACAGCACCGGAGCTTCAGCAATCAGCTCGCCCGACTTCAACTTg                                                                                                                              >  2:298300/1‑60 (MQ=255)                                                                                                          ccTTGACGATGCCACAGCACCGGAGCTTCAGCAATCAGCTTGCCCGACTTCAACTTGGCATGCGGGGTCAGGGCGCCA‑CCAGCACCCGTCGTACGGTCGCTCAAGAAGCGCACGGTACTGCCTGTCGCCGCGGGCGCTGAGTCGTTaaa                                 <  8:88051/149‑1 (MQ=255)                                                                                                          ccTTGACGATGCCACAGCACCGGAGCTTCAGCAATCAGCTTGCCCGACTTCAACTTGGCATGCGGGGTCAGGGCGCCA‑CCAGCACCCGTCGTACGGTCACTCAAGAAGCGCACGGTACTGCCTGTCGCCGCGGGCGCTGAGTCGTTaaa                                 <  3:122814/149‑1 (MQ=255)                                                                                                          ccTTGACGATGCCACAGCACCGGAGCTTCAGCAATCAGCTCGCCCGACTTCAACTTGGCATGCGGGGTCAGGGCGCCA‑CCAGCACCCGTCGTACGGTCGCTCAAGAAGCGCACGGTACTGCCTGTCGCCGCGGGCGCTGAGTCGTTaaa                                 >  6:344034/1‑149 (MQ=255)                                                                                                              gACGATGCCACAGCACCGGAGCTTCAGCAATCAGCTTGCCCGACTTCAACTTGGCATGCGGGGTCAGGGCGCCA‑CCAGCACCCGTCGTACGGTCGCTCAAGAAGCGCACGGTACTGCCTGTCGCCGCGGGCGCTGAGTCGTTAAACCCa                             >  3:178637/1‑149 (MQ=255)                                                                                                                      cacaGCACCGGAGCTTCAGCAATCAGCTTGCCCGACTTCAACTTGGCATGCGGGGTCAGGGCGCCA‑CCAGCACCCGTCGTACGGTCGCTCAAGAAGCGCACGGTACTGCCTGTCGCCGCGGGCGCTGAGTCGTTAAACCCATCGccac                      >  7:269706/1‑148 (MQ=255)                                                                                                                        caGCACCGGAGCTTCAGCAATCAGCTTGCGCGACTTCTACTTGGCATTCGGGGTCAGGGCGCCA‑CCAGCACCCGTCGTACGGTCGCTCAAGAAGCGCACGGTACTGCCTGTCGCcgcg                                                  <  4:245323/118‑1 (MQ=255)                                                                                                                         aGCACCGGAGCTTCAGCAATCAGCTTGCCCGACTTCAACTTGGCATGCGGGGTCAGGGCGCCA‑CCAGCACCCGTCGTACGGTCGCTCAAGAAGCGCACGGTACTGCCTGTCGCCGCGGGCGCt                                            <  5:274531/123‑1 (MQ=255)                                                                                                                         aGCACCGGAGCTTCAGCAATCAGCTTGCCCGACTTCAACTTGGCATGCGGGGTCAGGGCGCCA‑CCAGCACCCGTCGTACGGTCGCTCAAGAAGCGCACGGTACTGCCTGTCGCCGCGGGCGCt                                            >  6:274531/1‑123 (MQ=255)                                                                                                                         aGCACCGGAGCTTCAGCAATCAGCTTGCCCGACTTCAACTTGGCATGCGGGGTCAGGGCGCCA‑CCAGCACCCGTCGTACGGTCGCTCAAGAAGCGCACGGTACTGCCTGTCGCCGCGGGCGCTGAGTCGTTAAACCCATCGCCACCAgg                  >  5:91962/1‑149 (MQ=255)                                                                                                                            aCCGGAGCTTCAGCAATCAGCTTGCCCGACTTCAACTTGGCATGCGGGGTCAGGGCGCCATACAGCACCCGTCGTACGGTCGCTCAAGAAGCGCACGGTACTGCCTGTCGCCGCGGGCGCTGAGTCGTTAAACCCATCGCCACCAggg                 >  1:102610/1‑148 (MQ=255)                                                                                                                             ccGGAGCTTCAGCAATCAGCTTGCCCGACTTCAACTTGGCATGCGGGGTCAGGGCGCCA‑CCAGCACCCGTCGTACGGTCGCTCAAGAAGCGCACGGTACTGCCTGTCGCCGCGGGCGCTGAGTCGTTAAACCCATCGCCACCAGGGaa               >  2:45540/1‑148 (MQ=255)                                                                                                                               ggAGCTTCAGCAATCAGCTCGCCCGACTTCAACTTGGCATGCGGGGTCAGGGCGCCA‑CCAGCACCCGTCGTACGGTCGCTCAAGTAGCGCACGGTACTGCCTGTCGCCGCGGGCGCTGAGTCGTTAAACCCATCGCCACCAGGGAACAg            <  7:312938/149‑1 (MQ=255)                                                                                                                                caGCTTCAGCAATCAGCTTGCCCGACTTCAACTTGGCATGCGGGGTCAGGGCGCCA‑CCAGCACCCGTCGTACGGTCGCTCAAGAAGCGCACGGTACTGCCTGTCGCCGCGGGCGCTg                                           <  2:350329/116‑1 (MQ=255)                                                                                                                                     tCAGCAATCAGCTTGCCCGACTTCAACTTGGCATGCGGGGTCAGGGCGCCA‑CCAGCACCCGTCGTACGGTCGCTCAAGAAGCGCACGGTACTGCCTGTCGCCGCGGGCGCTGAGTCGt                                     >  5:125209/1‑118 (MQ=255)                                                                                                                                     tCAGCAATCAGCTTGCCCGACTTCAACTTGGCATGCGGGGTCAGGGCGCCA‑CCAGCACCCGTCGTACGGTCGCTCAAGAAGCGCACGGTACTGCCTGTCGCCGCGGGCGCTGAGTCGt                                     <  6:125209/118‑1 (MQ=255)                                                                                                                                       aGCAATCAGCTTGCCCGACTTCAACTTGGCATGCGGGGTCAGGGGGCCA‑CCAGCACCCGTCGTACGGTCGCTCAAGAAGCGCACGGTACTGCCTGTCGCCGCTGGCGCTGAGTCGTTAAACCCATCGCCACCAGGGAACAGGTCgcgc     <  4:286463/148‑1 (MQ=255)                                                                                                                                       aGCAATCAGCTTGCCCGACTTCAACTTGGCATGCGGGGTCAGGGCGCCA‑CCAGCACCCGTCGTACGGTCGCTCAAGAAGCGCACGGTACTGCCTGTCGCCGCGGGCGCTGAGTCGTTAAACCCATCGCCACCAGGGAACAGGTCgcgc     <  1:313812/148‑1 (MQ=255)                                                                                                                                         cAATCAGCTCGCCCGACTTCAACTTGGCATGCGGGGTCAGGGCGCCA‑CCAGCACCCGTCGTACGGTCGCTCAAGAAGCGCACGGTACTGCCTGTCGCCGCGGGCGCTGAGTCGTTAAACCCATCGCCACCAGGGAACAGGTCGCGCcac  <  5:344034/149‑1 (MQ=255)                                                                                                                                            tCAGCTTGCCCGACTTCAACTTGGCATGCGGGGTCAGGGCGCCA‑CCAGCACCCGTCGTACGGTCGCTCAAGAAGCGCACGGTACTGCCTGTCGCCGCGGGCGCTGAGTCGTTAAACCCATCGCCACCAGGGAACAGGTcgc       >  8:260226/1‑141 (MQ=255)                                                                                                                                            tCAGCTTGCCCGACTTCAACTTGGCATGCGGGGTCAGGGCGCCA‑CCAGCACCCGTCGTACGGTCGCTCAAGAAGCGCACGGTACTGCCTGTCGCCGCGGGCGCTGAGTCGTTAAACCCATCGCCACCAGGGAACAGGTcgc       <  7:260226/141‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                             GGAGCCTTGCTGAATACTTCGCCGAAAACGCCTTCTTCAAACAGCGGCAGGCCATTGGCGCCGATCATAAAGCCACGACTGCCCAGTGACAGCGATTGACCCGTGCCTTGACGATGCCACAGCACCGGAGCTTCAGCAATCAGCTCGCCCGACTTCAACTTGGCATGCGGGGTCAGGGCGCCA‑CCAGCACCCGTCGTACGGTCGCTCAAGAAGCGCACGGTACTGCCTGTCGCCGCGGGCGCTGAGTCGTTAAACCCATCGCCACCAGGGAACAGGTCGCGCCAC  >  NC_002947/4381627‑4381911

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