Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,866,099	Δ1 bp	5.4%	coding (604/645 nt)	PP_1669 →	NLP/P60 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,866,099	0	C	.	5.4%	242.2 / 6.9	74	coding (604/645 nt)	PP_1669	NLP/P60 family protein
Reads supporting (aligned to +/- strand):  ref base C (41/29);  new base . (2/2);  total (43/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.86e-01

CAGCACCAATGGCCGCGGCCGGGTAAGCCATGCCGGCATCTACCTGGGTGACAACCAGTTCATCCACTCCAGCAGCCGCCGCAGTGGTGGTGTGCGCATCGACAGCCTCGGTGATCGCTACTGGAGCAAGACCTTCATTGAAGCCAA‑GCGCGC‑CTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCGGTGATGCTTTGAAAAAAAAGCTCACCGCCGTGCGCATTTACAGAAAGCGTCTAATCTAAATTCTCAACTCTTTTTCGG  >  NC_002947/1865951‑1866246                                                                                                                                                      |                                                                                                                                                     cAGCACCAATGGCCGCGGCCGGGTAAGCCATGCCGGCATCTACCTGGGTGACAACCAGTTCATCCACTCCAGCAGCCGCCGCAGTGGTGGTGTGCGCATCGACAGCATCGGTGATCGCTACTGGAGCAAGACCTTCATTGAAGCCAA‑gc                                                                                                                                                      >  1:338976/1‑149 (MQ=255)  aGCACCAATGGCCGCGGCCGGGTAAGCCATGCCGGCATCTACCTGGGTGACAACCAGTTCATCCACTCCAGCAGCCGCCGCAGTGGTGGTGTGCGCATCGACAGCCTCGGTGATCGCTACTGGAGCAAGACCTTCATTGAAGCCAA‑gc                                                                                                                                                      <  6:259875/148‑1 (MQ=255)          tGGCCGCGGCCGGGTAAGCCATGCCGGCATCTACCTGGGTGACAACCAGTTCATCCACTCCAGCAGCCGCCGCAGTGGTGGTGTGCGCATCGACAGCCTCGGTGATCGCTACTGGAGCAAGACCTTCATTGAAGCCAA‑GCGCGC‑CTTgg                                                                                                                                            <  5:44419/149‑1 (MQ=255)              cgcgGCCGGGTAAGCCATGCCGGCATCTACCTGGGTGACAACCAGTTCATCCACTCCAGCAGCCGCCGCAGTGGTGGTGTGCGCATCGACAGCCTCGGTGATCGCTACTGGAGCAAGACCTTCATTGAAGCCAA‑GCGCGC‑CTTGGCCAt                                                                                                                                        <  4:366022/149‑1 (MQ=255)                            ccATGCCGGCATCTACCTGGGTGACAACCAGTTCATCCACTCCAGCAGCCGCCGCAGTGGTGGTGTGCGCATCGACAGCCTCGGTGATCGCTACTGGAGCAAGACCTTCATTGAAGCCAA‑GCGCGC‑CTTGGCCATGGCGCCGACCAAc                                                                                                                           >  3:73716/1‑148 (MQ=255)                                  cGGCATCTACCTGGGTGACAACCAGTTCATCCACTCCAGCAGCCGCCGCAGTGGTGGTGTGCGCATCGACAGCCTCGGTGATCGCTACTGGAGCAAGACCTTCATTGAAGCCAA‑GCGCGC‑CTTGGCCATGGCGCCGACCAACATcgcga                                                                                                                    <  4:73716/149‑2 (MQ=255)                                   ggCATCTACCTGGGTGACAACCAGTTCATCCACTCCAGCAGCCGCCGCAGTGGTGGTGTGCGCATCGACAGCCTCGGTGATCGCTACTGGAGCAAGACCTTCATTGAAGCCAA‑GCGCGC‑CTTGGCCATGGCGCCGACCAACATcgcgcg                                                                                                                   >  3:201552/1‑149 (MQ=255)                                      aTCTACCTGGGTGACAACCAGTTCATCCACTCCAGCAGCCGCCGCAGTGGTGGTGTGCGCATCGACAGCCTCGGTGATCGCTACTGGAGCAAGACCTTCATTGAAGCCAA‑GCGCGC‑CTTGGCCATGGCGCCGACCAACATCGCGCGCaa                                                                                                                >  5:97324/1‑149 (MQ=255)                                          aCCTGGGTGACAACCAGTTCATCCACTCCAGCAGCCGCCGCAGTGGTGGTGTGCGCATCGACAGCCTCGGTGATCGCTACTGGAGCAAGACCTTCATTGAAGCCAA‑GCGCGC‑CTTGGCCATgg                                                                                                                                      >  2:345951/1‑123 (MQ=255)                                           ccTGGGTGACAACCAGTTCATCCACTCCAGCAGCCGCCGCAGTGGTGGTGTGCGCATCGACAGCCTCGGTGATCGCTACTGGAGCAAGACCTTCATTGAAGCCAA‑GCGCGC‑CTTGGCCATGGCGCCGACCAACATCGCGCGCAAc                                                                                                               <  7:158970/145‑1 (MQ=255)                                           ccTGGGTGACAACCAGTTCATCCACTCCAGCAGCCGCCGCAGTGGTGGTGTGCGCATCGACAGCCTCGGTGATCGCTACTGGAGCAAGACCTTCATTGAAGCCAA‑GCGCGC‑CTTGGCCATGGCGCCGACCAACATCGCGCGCAACt                                                                                                              <  2:249849/146‑1 (MQ=255)                                                   aCAACCAGTTCATCCACTCCAGCAGCCGCCGCAGTGGTGGTGTGCGCATCGACAGCCTCGGTGATCGCTACTGGAGCAAGACCTTCATTGAAGCCAA‑GCGCGC‑CTTGGCCATGGCGCCGACCAACATcgcg                                                                                                                     >  4:13675/1‑131 (MQ=255)                                                   aCAACCAGTTCATCCACTCCAGCAGCCGCCGCAGTGGTGGTGTGCGCATCGACAGCCTCGGTGATCGCTACTGGAGCAAGACCTTCATTGAAGCCAA‑GCGCGC‑CTTGGCCATGGCGCCGACCAACATcgcg                                                                                                                     <  3:13675/131‑1 (MQ=255)                                                   aCAACCAGTTCATCCACTCCAGCAGCCGCCGCAGTGGTGGTGTGCGCATCGACAGCCTCGGTGATCGCTACTGGAGCAAGACCTTCATTGAAGCCAA‑GCGCGC‑CTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCaaaa                                                                                                    <  1:264500/148‑1 (MQ=255)                                                   aCAACCAGTTCATCCACTCCAGCAGCCGCCGCAGTGGTGGTGTGCGCATCGACAGCCTCGGTGATCGCTACTGGAGCAAGACCTTCATTGAAGCCAA‑GCGCGC‑CTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCaaaa                                                                                                    <  4:57181/148‑1 (MQ=255)                                                     aaCCAGTTCATCCACTCCAGCAGCCGCCGCAGTGGTGGTGTGCGCATCGACAGCCTCGGTGATCGCTACTGGAGCAAGACCTTCATTGAAGCCAA‑G‑GCGCGCTTGGctt                                                                                                                                         <  3:196280/109‑3 (MQ=255)                                                     aaCCAGTTCATCCACTCCAGCAGCCGCCGCAGTGGTGGTGTGCGCATCGACAGCCTCGGTGATCGCTACTGGAGCAAGACCTTCATTGAAGCCAA‑G‑GCGCGCTTGGctt                                                                                                                                         >  4:196280/1‑107 (MQ=255)                                                      aCCAGTTCATCCACTCCAGCAGCCGCCGCAGTGGTGGTGTGCGCATCGACAGCCTCGGTGATCGCTACTGGAGCAAGACCTTCATTGAAGCCAA‑GCGCGC‑CTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAatgat                                                                                                <  3:189601/149‑1 (MQ=255)                                                          gTTCATCCACTCCAGCAGCCGCCGCAGTGGTGGTGTGTGCATCGACAGCCTCGGTGATCGCTACTGGAGCAAGACCTTCATTGAAGCCAA‑GCGCGC‑CTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTAc                                                                                            >  6:422522/1‑149 (MQ=255)                                                          gTTCATCCACTCCAGCAGCCGCCGCAGTGGTGGTGTGCGCATCGACAGCCTCGGTGATCGCTACTGGAGCAAGACCTTCATTGAAGCCAA‑G‑GCGCGCTTGGctt                                                                                                                                         >  8:240839/1‑102 (MQ=255)                                                          gTTCATCCACTCCAGCAGCCGCCGCAGTGGTGGTGTGCGCATCGACAGCCTCGGTGATCGCTACTGGAGCAAGACCTTCATTGAAGCCAA‑G‑GCGCGCTTGGctt                                                                                                                                         <  7:240839/104‑3 (MQ=255)                                                               tCCACTCCAGCAGCCGCCGCAGTGGTGGTGTGCGCATCGACAGCCTCGGTGATCGCTACTGGAGCAAGACCTTCATTGAAGCCAA‑GCGCGC‑CTTg                                                                                                                                             >  5:270854/1‑95 (MQ=255)                                                               tCCACTCCAGCAGCCGCCGCAGTGGTGGTGTGCGCATCGACAGCCTCGGTGATCGCTACTGGAGCAAGACCTTCATTGAAGCCAA‑GCGCGC‑CTTg                                                                                                                                             <  6:270854/95‑1 (MQ=255)                                                               tCCACTCCAGCAGCCGCCGCAGTGGTGGTGTGCGCATCGACAGCCTCGGTGATCGCTACTGGAGCAAGACCTTCATTGAAGCCAA‑GCGCGC‑CTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGc                                                                                       >  3:297561/1‑149 (MQ=255)                                                                 cACTCCAGCAGCCGCCGCAGTGGTGGTGTGCGCATCGACAGCCTCGGTGATCGCTACTGGAGCAAGACCTTCATTGAAGCCAA‑GCGCGC‑CTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCcg                                                                                     >  2:127164/1‑149 (MQ=255)                                                                  aCTCCAGCAGCCGCCGCAGTGGTGGTGTGCGCATCGACAGCCTCGGTGATCGCTACTGGAGCAAGACCTTCATTGAAGCCAA‑GCGCGC‑CTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCcgc                                                                                    >  2:123412/1‑149 (MQ=255)                                                                     ccAGCAGCCGCCGCAGTGGTGGTGTGCGCATCGACAGCCTCGGTGATCGCTACTGGAGCAAGACCTTCATTGAAGCCAA‑GCGCGC‑CTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCgcggc                                                                                 >  5:143822/1‑149 (MQ=255)                                                                      cagcagCCGCCGCAGTGGTGGTGTGCGCATCGACAGCCTCGGTGATCGCTACTGGAGCAAGACCTTCATTGAAGCCAA‑GCGCGC‑CTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCgcggcg                                                                                >  5:90840/1‑149 (MQ=255)                                                                      cagcagCCGCCGCAGTGGTGGTGTGCGCATCGACAGCCTCGGTGATCGCTACTGGAGCAAGACCTTCATTGAAGCCAA‑GCGCGC‑CTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCgcgg                                                                                  <  4:201552/147‑1 (MQ=255)                                                                         cagcCGCCGCAGTGGTGGTGTGCGCATCGACAGCCTCGGTGATCGCTACTGGAGCAAGACCTTCATTGAAGCCAA‑GCGCGC‑CTTGGCCATGGCGTCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCGGt                                                                              >  2:45020/1‑148 (MQ=255)                                                                         cagcCGCCGCAGTGGTGGTGTGCGCATCGACAGCCTCGGTGATCGCTACTGGAGCAAGACCTTCATTGAAGCCAA‑GCGCGC‑CTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCGGt                                                                              >  2:305141/1‑148 (MQ=255)                                                                                cgcAGTGGTGGTGTGCGCATCGACAGCCTCGGTGATCGCTACTGGAGCAAGACCTTCATTGAAGCCAA‑GCGCGC‑CTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCGGTGATGCttt                                                                      <  7:227951/149‑1 (MQ=255)                                                                                cgcAGTGGTGGTGTGCGCATCGACAGCCTCGGTGATCGCTACTGGAGCAAGACCTTCATTGAAGCCAA‑GCGCGC‑CTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCGGTGATGCttt                                                                      >  3:365584/1‑149 (MQ=255)                                                                                 gcAGTGGTGGTGTGCGCATCGACAGCCTCGGTGATCGCTACTGGAGCAAGACCTTCATTGAAGCCAA‑GCGCGC‑CTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCGGTGATGCTTTg                                                                     <  1:305141/149‑1 (MQ=255)                                                                                       gtggtgTGCGCATCGACAGCCTCGGTGATCGCTACTGGAGCAAGACCTTCATTGAAGCCAA‑GCGCGC‑CTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCGGTGATGCTTTGaaaaa                                                                >  6:167404/1‑148 (MQ=255)                                                                                        tggtgTGCGCATCGACAGCCTCGGTGATCGCTACTGGAGCAAGACCTTCATTGAAGCCAA‑GCGCGC‑CTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCGGTGATGCTTTGaaaaaa                                                               >  4:329275/1‑148 (MQ=255)                                                                                         ggtgTGCGCATCGACAGCCTCGGTGATCGCTACTGGAGCAAGACCTTCATTGAAGCCAA‑GCGCGC‑CTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCGGTGATGCTTTGaaaaaaaa                                                             <  7:79864/149‑1 (MQ=255)                                                                                             tgCGCATCGACAGCCTCGGTGATCGCTACTGGAGCAAGACCTTCATTGAAGCCAA‑GCGCGC‑CTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCGGTGATGCTTTGAAAAAAAAGCTc                                                         <  6:90840/149‑1 (MQ=255)                                                                                              gcgcATCGACAGCCTCGGTGATCGCTACTGGAGCAAGACCTTCATTGAAGCCAA‑GCGCGC‑CTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCGGTGATGCTTTGAAAAAAAAGCTCa                                                        >  8:21605/1‑149 (MQ=255)                                                                                                  aTCGACAGCCTCGGTGATCGCTACTGGAGCAAGACCTTCATTGAAGCCAA‑GCGCGC‑CTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCAGCGGTGATGCTTTGAAAAAAAAGCTCAccgc                                                    >  6:78447/1‑149 (MQ=255)                                                                                                     gACAGCCTCGGTGATCGCTACTGGAGCAAGACCTTCATTGAAGCCAA‑GCGCGC‑CTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCGGTGATGCTTTGAAAAAAAAGCTCAccgccg                                                  <  5:422522/148‑1 (MQ=255)                                                                                                     gACAGCCTCGGTGATCGCTACTGGAGCAAGACCTTCATTGAAGCCAA‑GCGCGC‑CTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCGGTGATGCTTTGAAAAAAAAGCTCACCGCCGt                                                 >  8:47433/1‑149 (MQ=255)                                                                                                       cAGCCTCGGTGATCGCTACTGGAGCAAGACCTTCATTGAAGCCAA‑GCGCGC‑CTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCGGTGATGCTTTGAAAAAAAAGCTCACCGCCGTgc                                               <  5:167404/149‑1 (MQ=255)                                                                                                       cAGCCTCGGTGATCGCTACTGGAGCAAGACCTTCATTGAAGCCAA‑GCGCGC‑CTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCGGTGATGCTTTGAAAAAAAAGCTCACCGCCGTgc                                               <  3:329275/149‑1 (MQ=255)                                                                                                          ccTCGGTGATCGCTACTGGAGCAAGACCTTCATTGAAGCCAA‑GCGCGC‑CTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCaaa                                                                                                     >  3:189725/1‑92 (MQ=255)                                                                                                          ccTCGGTGATCGCTACTGGAGCAAGACCTTCATTGAAGCCAA‑GCGCGC‑CTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCaaa                                                                                                     <  4:189725/92‑1 (MQ=255)                                                                                                          ccTCGGTGATCGCTACTGGAGCAAGACCTTCATTGAAGCCAA‑GCGCGC‑CTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCGGTGATGCTTTGAAAAAAAAGCTCACCGCCGTGCGCa                                            >  7:41761/1‑149 (MQ=255)                                                                                                            tCGGTGATCGCTACTGGAGCAAGACCTTCATTGAAGCCAA‑GCGCGC‑CTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCGGTGATGCTTTGAAAAAAAAGCTCACCGCCGTGCGCAtt                                          >  8:272899/1‑149 (MQ=255)                                                                                                                    cGCTACTGGAGCAAGACCTTCATTGAAGCCAA‑GCGCGC‑CTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCGGTGATGCTTTGAAAAAAAAGCTCACCGCCGTGCGCATTTACAGa                                    >  5:217666/1‑147 (MQ=255)                                                                                                                     gCTACTGGAGCAAGACCTTCATTGAAGCCAA‑GCGCGC‑CTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCGGTGATGCTTTGAAAAAAAAGCTCACCGCCGTGCGCAtt                                          >  1:379561/1‑140 (MQ=255)                                                                                                                     gCTACTGGAGCAAGACCTTCATTGAAGCCAA‑GCGCGC‑CTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCGGTGATGCTTTGAAAAAAAAGCTCACCGCCGTGCGCAtt                                          <  2:379561/140‑1 (MQ=255)                                                                                                                       tACTGGAGCAAGACCTTCATTGAAGCCAA‑GCGCGC‑CTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCGGTGATGCTTTGAAAAAAAAGCTCACCGCCGTGCGCATTTACAGAAAGCg                               >  7:135232/1‑149 (MQ=255)                                                                                                                          tGGAGCAAGACCTTCATTGAAGCCAA‑GCGCGC‑CTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCGGTGATGCTTTGAAAAAAAAGCTCACCGCCGTGCGCATTTACAGAAAGCGtc                             >  4:212075/1‑148 (MQ=255)                                                                                                                            gAGCAAGACCTTCATTGAAGCCAA‑GCGCGC‑CTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCGGTGATGCTTTGAAAAAAAAGCTCACCGCCGTGCGCATTTACAGAAAGCGtcta                           <  7:21605/148‑1 (MQ=255)                                                                                                                               cAAGACCTTCATTGAAGCCAA‑GCGCGC‑CTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCGGTGATGCTTTGAAAAAAAAGCTCACCGCCGTGCGCATTTACAGAAAGCGtctaatc                        >  4:362749/1‑148 (MQ=255)                                                                                                                               cAAGACCTTCATTGAAGCCAA‑GCGCGC‑CTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCGGTGATGCTTTGAAAAAAAAGCTCACCGCCGTGCGCATTTACAGAAAGCGtctaatc                        >  5:87533/1‑148 (MQ=255)                                                                                                                                aaGACCTTCATTGAAGCCAA‑GCGCGC‑CTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCGGTGATGCTTTGAAAAAAAAGCTCACCGCCGTGCGCATTTACAGAAAGCGtctaatct                       >  3:402413/1‑148 (MQ=255)                                                                                                                                     cTTCATTGAAGCCAA‑GCGCGC‑CTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCGGTGATGCTTTGAAAAAAAAGCTCACCGCCGTGCGCATTTACAGaa                                   <  4:255564/131‑1 (MQ=255)                                                                                                                                     cTTCATTGAAGCCAA‑GCGCGC‑CTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCGGTGATGCTTTGAAAAAAAAGCTCACCGCCGTGCGCATTTACAGaa                                   >  3:255564/1‑131 (MQ=255)                                                                                                                                     cTTCATTGAAGCCAA‑GCGCGC‑CTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCGGTGATGCTTTGAAAAAAAAGCTCACCGCCGTGCGCATTTACAGAAAGc                                <  3:275848/134‑1 (MQ=255)                                                                                                                                     cTTCATTGAAGCCAA‑GCGCGC‑CTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCGGTGATGCTTTGAAAAAAAAGCTCACCGCCGTGCGCATTTACAGAAAGc                                >  4:275848/1‑134 (MQ=255)                                                                                                                                     cTTCATTGAAGCCAA‑GCGCGC‑CTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCGGTGATGCTTTGAAAAAAAAGCTCACCGCCGTGCGCATTTACAGAAAGCGTCTAATCTAAATtc                 >  3:67007/1‑149 (MQ=255)                                                                                                                                      ttCATTGAAGCCAA‑GCGCGC‑CTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCGGTGATGCTTTGAAAAAAAAGCTCACCGCCGTGCGCATTTACAGAAAGCGTCTAATCTAAATtc                 >  8:31876/1‑148 (MQ=255)                                                                                                                                      ttCATTGAAGCCAA‑GCGCGC‑CTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCGGTGATGCTTTGAAAAAAAAGCTCACCGCCGTGCGCATTTACAGAAAGCGTCTAATCTAAATtc                 >  2:66265/1‑148 (MQ=255)                                                                                                                                        cATTGAAGCCAA‑GCGCGC‑CTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCGGTGATGCTTTGAAAAAAAAGCTCACCGCCGTGCGCATTTACAGAAAGCGTCTAATCTAAATtct                <  6:87533/147‑1 (MQ=255)                                                                                                                                         aTTGAAGCCAA‑GCGCGC‑CTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCGGTGATGCTTTGAAAAAAAAGCTCACCGCCGTGCGCATTTACAGAAAGCGTCTAATCTAAATTCTCaa             <  8:135232/149‑1 (MQ=255)                                                                                                                                         aTTGAAGCCAA‑GCGCGC‑CTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCGGTGATGCTTTGAAAAAAAAGCTCACCGCCGTGCGCATTTACAGAAAGCGTCTAATCTAAATTCTCaa             <  4:67007/149‑1 (MQ=255)                                                                                                                                           catggCCAAGGCGCG‑‑CTTGGCTATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCGGTGATGCTTTGAAAAAAAAGCTCACCGCCGTGCGCATTTACAGAAAGCGTCTAATCTAAATTCTCAAct           >  7:313869/5‑149 (MQ=255)                                                                                                                                               gCCAA‑GCGCGC‑CTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCTGCGGTGATGCTTTGAAAAAAAAGCTCACCGCCGTGCGCATTTACAGAAAGCGTCTAATCTAAATTCTCAACTCttt       <  5:78447/149‑1 (MQ=255)                                                                                                                                               gCCAA‑GCGCGC‑CTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCGGTGATGCTTTGAAAAAAAAGCTCACCGCCGTGCGCATTTACAGAAAGCGTCTAATCTAAATTCTCAACTCttt       >  2:219691/1‑149 (MQ=255)                                                                                                                                               gCCAA‑GCGCGC‑CTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGAACCCTGCCGCGGCGGTGATGCTTTGAAAAAAAAGCTCACCGCCGTGCGCATTTACAGAAAGCGTCTAATCTAAATTCTCAACTCttt       >  5:89633/1‑149 (MQ=255)                                                                                                                                                 cAA‑GCGCGC‑CTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCGGTGATGCTTTGAAAAAAAAGCTCACCGCCGTGCGCATTTACAGAAAGCGTCTAATCTAAATTCTCAACTCttttt     >  6:353540/1‑149 (MQ=255)                                                                                                                                                  aa‑GCGCGC‑CTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGGAATCAAAATGATGTACCCTGCCGCGGCGGTGATGCTTTGAAAAAAAAGCTCACCGCCGTGCGCATTTACAGAAAGCGTCTAATCTAAATTCTCAACTCTTTTTc    >  1:291505/1‑149 (MQ=255)                                                                                                                                                  aa‑GCGCGC‑CTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGTCGCGGCGGTGATGCTTTGAAAAAAAAGCTCACCGCCGTGCGCATTTACAGAAAGCGTCTAATCTAAATTCTCAACTCTTTTTc    <  2:5924/149‑1 (MQ=255)                                                                                                                                                  aa‑GCGCGC‑CTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCGGTGATGCTTTGAAAAAAAAGCTCACCGCCGTGCGCATTTACAGAAAGCGTCTAATCTAAATTCTCAACTCTTTTTc    >  4:260987/1‑149 (MQ=255)                                                                                                                                                  aa‑GCGCGC‑CTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCGGTGATGCTTTGAAAAAAAAGCTCACCGCCGTGCGCATTTACAGAAAGCGTCTAATCTAAATTCTCAACTCTTTTTc    <  6:143822/149‑1 (MQ=255)                                                                                                                                                     gcgcgc‑CTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCGGTGATGCTTTGAAAAAAAAGCTCACCGCCGTGCGCATTTACAGAAAGCGTCTAATCTAAATTCTCAACTCTTTTTCgg  <  1:123412/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                     CAGCACCAATGGCCGCGGCCGGGTAAGCCATGCCGGCATCTACCTGGGTGACAACCAGTTCATCCACTCCAGCAGCCGCCGCAGTGGTGGTGTGCGCATCGACAGCCTCGGTGATCGCTACTGGAGCAAGACCTTCATTGAAGCCAA‑GCGCGC‑CTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCGGTGATGCTTTGAAAAAAAAGCTCACCGCCGTGCGCATTTACAGAAAGCGTCTAATCTAAATTCTCAACTCTTTTTCGG  >  NC_002947/1865951‑1866246

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