Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,431,271	G→C	6.2%	L347F (TTG→TTC) ‡	PP_3925 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,431,271	0	G	C	6.2%	257.5 / 5.8	80	L347F (TTG→TTC) ‡	PP_3925	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (42/33);  new base C (2/3);  total (44/36)
Fisher's exact test for biased strand distribution p-value = 6.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.14e-01

ACAACTTCGATTACTACTATTGCTATGGCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCGGCGAGGCCATCGCC  >  NC_002947/4431123‑4431418                                                                                                                                                     |                                                                                                                                                    aCAACTTCGATTACTACTATTGCTATGGCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTg                                                                                                                                                     >  2:13060/1‑149 (MQ=255)   aaCTTCGATTACTACTATTGCTATGGCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGc                                                                                                                                                    >  3:311178/1‑148 (MQ=255)   aaCTTCGATTACTACTATTGCTATGGCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGc                                                                                                                                                    <  4:311178/148‑1 (MQ=255)    aCTTCGATTACTACTATTGCTATGGCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAg                                                                                                                                                  <  8:194560/149‑1 (MQ=255)      ttCGATTACTACTATTGCTATGGCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGt                                                                                                                                                >  8:220993/1‑149 (MQ=255)         gATTACTACTATTGCTATGGCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCa                                                                                                                                              <  3:170372/148‑1 (MQ=255)          aTTACTACTATTGCTATGGCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTGCAAGGTCAtt                                                                                                                                            <  1:243805/149‑1 (MQ=255)           ttACTACTATTGCTATGGCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTGCAAGGTCAttt                                                                                                                                           <  8:17066/149‑1 (MQ=255)               tactaTTGCTATGGCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTGCAAGGTCAtttgttc                                                                                                                                       <  8:61754/149‑5 (MQ=255)                  taTTGCTATGGCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTc                                                                                                                                    >  7:148974/1‑149 (MQ=255)                    ttGCTATGGCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTc                                                                                                                                    >  1:333626/1‑147 (MQ=255)                    ttGCTATGGCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCa                                                                                                                                  <  8:271931/149‑1 (MQ=255)                          tGGCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCctct                                                                                                                            <  1:261994/149‑1 (MQ=255)                            gCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTg                                                                                                                           <  1:13060/148‑1 (MQ=255)                              cTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAg                                                                                                                         >  6:334167/1‑148 (MQ=255)                               tGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTcc                                                                                                                                   <  2:48022/137‑1 (MQ=255)                               tGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTcc                                                                                                                                   >  1:48022/1‑137 (MQ=255)                                    cGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGcc                                                                                                                  >  2:250163/1‑149 (MQ=255)                                     gCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGcc                                                                                                                  >  4:258563/1‑148 (MQ=255)                                        tGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCac                                                                                                               >  8:249882/1‑148 (MQ=255)                                          gcaggcaTCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCacg                                                                                                              <  4:203202/147‑1 (MQ=255)                                             ggcaTCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCacgacgac                                                                                                         <  8:106776/149‑1 (MQ=255)                                              gcaTCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCacgacgacg                                                                                                        <  1:256759/149‑1 (MQ=255)                                               caTCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCacgacgacg                                                                                                        >  2:205572/1‑148 (MQ=255)                                                   gTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATaa                                                                                                   <  5:238678/149‑1 (MQ=255)                                                    ttCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATaa                                                                                                   >  2:159507/1‑148 (MQ=255)                                                      cagcagATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACa                                                                                                 <  7:220993/148‑1 (MQ=255)                                                        gcagATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAgg                                                                                                                                                 <  3:119931/98‑1 (MQ=255)                                                        gcagATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAgg                                                                                                                                                 >  4:119931/1‑98 (MQ=255)                                                         cagATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAgg                                                                                              >  7:95287/1‑148 (MQ=255)                                                          agATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCTTTTCGCACAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCAAGCATCTGAGCACAGCCCCCGACGACTGATAACAAGGa                                                                                             >  3:237004/1‑148 (MQ=255)                                                               tactacCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAg                                                                                       >  8:296169/1‑149 (MQ=255)                                                                 ctacCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGc                                                                                      <  8:95287/148‑1 (MQ=255)                                                                    ccGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAAc                                                                                         <  3:90965/142‑1 (MQ=255)                                                                    ccGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAAc                                                                                         >  4:90965/1‑142 (MQ=255)                                                                     cGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTc                                                                                 <  8:148974/149‑1 (MQ=255)                                                                          ctacCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGAGCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGa                                                                            <  7:105934/149‑1 (MQ=255)                                                                          ctacCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGa                                                                            <  1:250163/149‑1 (MQ=255)                                                                           tacCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAACCAGCATGTCCAAGAc                                                                           >  4:29252/1‑149 (MQ=255)                                                                             ccATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGAt                                                                                                     <  8:212156/121‑1 (MQ=255)                                                                             ccATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGAcc                                                                          >  6:21001/1‑148 (MQ=255)                                                                             ccATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCAGTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGAt                                                                                                     >  7:212156/1‑121 (MQ=255)                                                                               aTGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCAc                                                                       <  5:71965/149‑1 (MQ=255)                                                                                   ccAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGc                                                                                                                                >  2:207631/1‑88 (MQ=255)                                                                                   ccAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGc                                                                                                                                <  1:207631/88‑1 (MQ=255)                                                                                   ccAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGt                                                                   >  8:305583/1‑149 (MQ=255)                                                                                    cAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGtt                                                                  <  5:21001/149‑1 (MQ=255)                                                                                        cccAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGAc                                                              >  2:268988/1‑149 (MQ=255)                                                                                         ccAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGAcc                                                             >  8:98853/1‑149 (MQ=255)                                                                                         ccAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGAc                                                              >  5:228918/1‑148 (MQ=255)                                                                                         ccAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGAc                                                              >  8:94097/1‑148 (MQ=255)                                                                                           aGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTc                                                           >  3:247072/1‑149 (MQ=255)                                                                                             gACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGa                                                         <  1:205572/149‑1 (MQ=255)                                                                                                 aGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGc                                                     <  1:268988/149‑1 (MQ=255)                                                                                                  gCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGCCGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGc                                                     >  6:182000/1‑148 (MQ=255)                                                                                                  gCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGAc                                                              <  4:161219/139‑1 (MQ=255)                                                                                                  gCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGAc                                                              >  3:161219/1‑139 (MQ=255)                                                                                                    gTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCaa                                                   <  6:228918/148‑1 (MQ=255)                                                                                                          cccAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGcc                                            >  6:320147/1‑149 (MQ=255)                                                                                                           ccAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCt                                           >  3:217661/1‑149 (MQ=255)                                                                                                             aGTTCATTCACATGAACCGCTTGCTTGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCttt                                         <  5:334167/149‑1 (MQ=255)                                                                                                                  attcaCATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTc                                    <  1:159507/149‑1 (MQ=255)                                                                                                                   ttgacaTGAACGGCTTGCTGGAACAATTGACCTTGCAGGTCATTTCCAAGTCCTGCCTCTGAGCACAGCCCAGGACGACGGATAACAAGGAAATGAGCATGTGCAAGACGCTCCTGTTCGTCCTCGACGGGAAGATTGCCTTTGTTTcc                                   <  6:187312/146‑1 (MQ=255)                                                                                                                        aTGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGgcgc                              >  3:134793/1‑149 (MQ=255)                                                                                                                           aaCCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCa                            <  7:249882/148‑1 (MQ=255)                                                                                                                            aCCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCAACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGc                          >  8:233403/1‑149 (MQ=255)                                                                                                                              cGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCg                        >  4:33991/1‑149 (MQ=255)                                                                                                                                cTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTg                      <  5:320147/149‑1 (MQ=255)                                                                                                                                   gCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGCGCACAGCCAACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCACGACGGCAAGAATGCCTAAGATTCCGGCGCCAGCCGTGGc                    >  7:296205/1‑148 (MQ=255)                                                                                                                                     tGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATc                 >  3:195323/1‑149 (MQ=255)                                                                                                                                       ttggaaaTGACCTGCAAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCgg               >  6:255332/5‑149 (MQ=255)                                                                                                                                       ttggaaaTGACCTGCAAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCgg               >  1:143470/5‑149 (MQ=255)                                                                                                                                       gAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGa                                                         >  2:16002/1‑107 (MQ=255)                                                                                                                                       gAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGa                                                         <  1:16002/107‑1 (MQ=255)                                                                                                                                           aaaTGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCGGCGAg           <  4:223562/149‑1 (MQ=255)                                                                                                                                           aaaTGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCGGCGAg           >  4:257488/1‑149 (MQ=255)                                                                                                                                            aaTGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCGGCGAgg          >  6:32449/1‑149 (MQ=255)                                                                                                                                              tGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCGGCGAGGc         >  4:153813/1‑148 (MQ=255)                                                                                                                                               gACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCGGCGAGGcc        >  8:265609/1‑148 (MQ=255)                                                                                                                                               gACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCGGCGAGGCCa       >  2:55271/1‑149 (MQ=255)                                                                                                                                                aCCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCGGCGAGGCCAt      >  4:17478/1‑149 (MQ=255)                                                                                                                                                     gCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCCTCGCCGATGCCacggct  <  4:247072/148‑6 (MQ=255)                                                                                                                                                     |                                                                                                                                                    ACAACTTCGATTACTACTATTGCTATGGCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCGGCGAGGCCATCGCC  >  NC_002947/4431123‑4431418

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