Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,431,270	T→G	6.1%	L347W (TTG→TGG) ‡	PP_3925 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,431,270	0	T	G	6.1%	258.8 / 6.3	82	L347W (TTG→TGG) ‡	PP_3925	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (40/37);  new base G (3/2);  total (43/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.26e-01

ACAACTTCGATTACTACTATTGCTATGGCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCGGCGAG  >  NC_002947/4431123‑4431409                                                                                                                                                    |                                                                                                                                            aCAACTTCGATTACTACTATTGCTATGGCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTg                                                                                                                                            <  1:125386/149‑1 (MQ=255)      ttCGATTACTACTATTGCTATGGCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAgg                                                                                                                                        >  3:165576/1‑148 (MQ=255)      ttCGATTACTACTATTGCTATGGCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGt                                                                                                                                       >  8:231660/1‑149 (MQ=255)          aTTACTACTATTGCTATGGCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCAtt                                                                                                                                   <  5:245289/149‑1 (MQ=255)           ttACTACTATTGCTATGGCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCAttt                                                                                                                                  >  7:181477/1‑149 (MQ=255)                 ctaTTGCTATGGCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCaa                                                                                                                              <  5:316758/147‑1 (MQ=255)                  taTTGCTATGGCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTc                                                                                                                           >  3:221846/1‑149 (MQ=255)                        tATGGCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCct                                                                                                                     >  2:177122/1‑149 (MQ=255)                         aTGGCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCctc                                                                                                                    >  7:259383/1‑149 (MQ=255)                          tGGCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCctct                                                                                                                   >  1:340501/1‑149 (MQ=255)                            gCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGa                                                                                                                 <  3:177114/149‑1 (MQ=255)                            gCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGa                                                                                                                 <  1:74898/149‑1 (MQ=255)                            gCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGa                                                                                                                 <  8:259383/149‑1 (MQ=255)                               tGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGca                                                                                                              >  1:322959/1‑149 (MQ=255)                               tGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGc                                                                                                               >  5:23344/1‑148 (MQ=255)                                    cGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGcc                                                                                                         >  4:317992/1‑149 (MQ=255)                                     gCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGcc                                                                                                         >  1:345478/1‑148 (MQ=255)                                       cTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCac                                                                                                      <  2:340501/149‑1 (MQ=255)                                                  cGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCacgacgacg                                                                                               <  3:222468/145‑1 (MQ=255)                                                  cGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCacgacgacg                                                                                               >  4:222468/1‑145 (MQ=255)                                                    ttCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAAc                                                                                         <  2:345478/149‑1 (MQ=255)                                                       agcagATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACaa                                                                                       <  7:322146/148‑1 (MQ=255)                                                         cagATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGa                                                                                    <  4:221846/149‑1 (MQ=255)                                                          agATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGaa                                                                                   <  6:23344/149‑1 (MQ=255)                                                             tCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGa                                                                                             >  4:183110/1‑136 (MQ=255)                                                             tCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGa                                                                                             <  3:183110/136‑1 (MQ=255)                                                             tCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGaaaa                                                                                 >  2:200111/1‑148 (MQ=255)                                                              ctactaCCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGTAAAACa                                                                               <  8:181477/149‑1 (MQ=255)                                                               tactacCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAg                                                                              >  5:376704/1‑149 (MQ=255)                                                                actacCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGAAGACGGATAACAAGGAAAACAg                                                                              >  8:314448/1‑148 (MQ=255)                                                                 ctacCGCTACTCCCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAg                                                                              <  7:150048/147‑1 (MQ=255)                                                                 ctacCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCa                                                                            <  1:177122/149‑1 (MQ=255)                                                                 ctacCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCa                                                                            <  4:165576/149‑1 (MQ=255)                                                                    ccGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGt                                                                         >  3:330087/1‑149 (MQ=255)                                                                     cGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCAAGTTCAATCACATGAACCGCTTGATGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACCAAGACGGATAACAAGGAAAACAGCATGTc                                                                        >  1:288540/1‑149 (MQ=255)                                                                       ctactaCCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTcc                                                                       >  3:208481/1‑148 (MQ=255)                                                                             ccATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTGCAAGGTCAtttgttc                                                                                                                              <  7:243275/87‑5 (MQ=38)                                                                             ccATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTGCAAGGTCAtttgttc                                                                                                                              >  8:243275/1‑83 (MQ=38)                                                                               aTGGCCAGACCCAGGAGAAGCGTTTAGCCGAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAATTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACTAGGAAAACAGGATGTCCAAGACCGAc                                                              <  2:288540/149‑1 (MQ=255)                                                                                   ccAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGt                                                          >  4:4988/1‑149 (MQ=255)                                                                                    cAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGtt                                                         >  1:116411/1‑149 (MQ=255)                                                                                    cAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGtt                                                         >  5:81166/1‑149 (MQ=255)                                                                                     aGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGAAGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTc                                                        <  7:314448/149‑1 (MQ=255)                                                                                      gACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCg                                                       <  4:208481/149‑1 (MQ=255)                                                                                        cccAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGAc                                                     <  2:111250/149‑1 (MQ=255)                                                                                        cccAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGAc                                                     >  6:366879/1‑149 (MQ=255)                                                                                         ccAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGAcc                                                    >  6:220586/1‑149 (MQ=255)                                                                                           aGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGa                                                                                                                 <  6:35945/86‑1 (MQ=255)                                                                                           aGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGa                                                                                                                 >  5:35945/1‑86 (MQ=255)                                                                                           aGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACa                                                                               >  4:345121/1‑120 (MQ=255)                                                                                           aGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACa                                                                               <  3:345121/120‑1 (MQ=255)                                                                                              aCAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACaa                                                                                       >  5:255394/1‑109 (MQ=255)                                                                                              aCAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACaa                                                                                       <  6:255394/109‑1 (MQ=255)                                                                                                aaGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACgg                                             >  8:316667/1‑149 (MQ=255)                                                                                                aaGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACgg                                             >  2:104395/1‑149 (MQ=255)                                                                                                 aGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACgg                                             >  3:54347/1‑148 (MQ=255)                                                                                                   cGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCaa                                          <  6:376704/149‑1 (MQ=255)                                                                                                          cccAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGcc                                   <  3:338663/149‑1 (MQ=255)                                                                                                            cAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCtt                                 <  3:126448/149‑1 (MQ=255)                                                                                                             aGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCttt                                <  1:104395/149‑1 (MQ=255)                                                                                                                  attcaCATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGttt                            <  2:322959/148‑1 (MQ=255)                                                                                                                   ttcaCATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTcc                          >  4:109750/1‑149 (MQ=255)                                                                                                                      acaTGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTc                                                                        <  7:142584/100‑1 (MQ=255)                                                                                                                      acaTGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGgc                       >  1:278956/1‑149 (MQ=255)                                                                                                                        aTGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTc                                                        <  5:102147/114‑1 (MQ=255)                                                                                                                        aTGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTc                                                        >  6:102147/1‑114 (MQ=255)                                                                                                                          gAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCGTTGTTTCCGgcgc                     >  6:367647/1‑147 (MQ=255)                                                                                                                             ccGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGcc                >  8:258315/1‑149 (MQ=255)                                                                                                                              cGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGGTTGCCTTTGTTTCCGGCGCCAGCCg               >  8:282693/1‑149 (MQ=255)                                                                                                                               gCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGt              <  1:200111/149‑1 (MQ=255)                                                                                                                                cTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTg             <  7:316667/149‑1 (MQ=255)                                                                                                                                   gCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGa                                                <  2:31934/111‑1 (MQ=255)                                                                                                                                   gCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGa                                                >  1:31934/1‑111 (MQ=255)                                                                                                                                    cTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCAt         >  6:338290/1‑149 (MQ=255)                                                                                                                                       ttggaaaTGACCTGCAAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTg             >  2:286956/5‑142 (MQ=255)                                                                                                                                       ttggaaaTGACCTGCAAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTg             <  1:286956/138‑1 (MQ=255)                                                                                                                                       gAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCgg      >  1:176723/1‑149 (MQ=255)                                                                                                                                       gAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCgg      <  2:116411/149‑1 (MQ=255)                                                                                                                                       gAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCgg      >  2:75261/1‑149 (MQ=255)                                                                                                                                       gAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCg       <  4:54347/148‑1 (MQ=255)                                                                                                                                        aaCAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCGGc     <  3:317992/149‑1 (MQ=255)                                                                                                                                           aaaTGACCTGCAAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTCCCGGCGCCAGCCGTGGCATCGGCGAg  >  7:212397/1‑149 (MQ=255)                                                                                                                                                   ttGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACTGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCGGc     <  6:95700/138‑1 (MQ=255)                                                                                                                                                   ttGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCGGc     >  5:95700/1‑138 (MQ=255)                                                                                                                                                    |                                                                                                                                            ACAACTTCGATTACTACTATTGCTATGGCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCGGCGAG  >  NC_002947/4431123‑4431409

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