Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	979,599	A→T	6.6%	H350L (CAC→CTC)	iscS‑I →	cysteine desulfurase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	979,599	0	A	T	6.6%	184.0 / 2.4	66	H350L (CAC→CTC)	iscS‑I	cysteine desulfurase
Reads supporting (aligned to +/- strand):  ref base A (34/24);  new base T (2/2);  total (39/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

CTGAGCTTCAACTACGTCGAAGGCGAATCGCTGCTGATGTCGCTGAAAGACATCGCCGTATCGTCCGGTTCGGCCTGTACCTCCGCGTCGCTCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTG‑GCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGCGAGGCCGTCAATAAACTGCGTGAGCTGTCGCCGCTGTGGGATATGTACAAAGACGGCGTTGACATCTCCAAGATC  >  NC_002947/979457‑979747                                                                                                                                                |                                                                                                                                                     ctgaGCTTCAACTACGTCGAAGGCGAATCGCTGCTGATGTCGCTGAAAGACATCGCCGTATCGTCCGGTTCGGCCTGTACCTCCGCGTCGCTCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTG‑GCGCACAGCTc                                                                                                                                                >  1:79067/1‑149 (MQ=255)         cAACTACGTCGGAGGCGAATCGCTGCTGATGTCGCTGAAAGACATCGCCGTATCGTCCGGTTCGGCCTGTACCTCCGCGTCGCTCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGGAACGACGAGCTG‑GCGCACAGCTCGATCCGc                                                                                                                                         <  1:30589/148‑1 (MQ=255)            cTACGTCGAAGGCGAATCGCTGCTGATGTCGCTGAAAGACATCGCCGTATCGTCCGGCTCGGCCTGTACCTCCGCGTCGCTCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTGTGCGC‑CAGCTCgtcgt                                                                                                                                           >  5:40502/1‑139 (MQ=255)            cTACGTCGAAGGCGAATCGCTGCTGATGTCGCTGAAAGACATCGCCGTATCGTCCGGCTCGGCCTGTACCTCCGCGTCGCTCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTGTGCGC‑CAGCTCgtcgt                                                                                                                                           <  6:40502/143‑5 (MQ=255)                  cGAAGGCGAATCGCTGCTGATGTCGCTGAAAGACATCGCCGTATCGTCCGGTTCGGCCTGTACCTACGCGTCGCTCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTG‑GCGCACAGCTCGATCCGCTTCTCCTTCg                                                                                                                               >  5:122452/1‑149 (MQ=255)                      ggCGAATCGCTGCTGATGTCGCTGAAAGACATCGCCGTATCGTCCGGTTCGGCCTGTACCTCCGCGTCGCTCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTG‑GCGCACAGCTCGATCCGCTTCTCCTTCGGCCg                                                                                                                           >  4:145242/1‑149 (MQ=255)                         gAATCGCTGCTGATGTCGCTGAAAGACATCGCCGTATCGTCCGGTTCGGCCTGTACCTCCGCGTCGCTCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTG‑GCGCACAGCTCGATCCGCTTCTCCTTCGGCCGtt                                                                                                                         >  6:96691/1‑148 (MQ=255)                                tgctgATGTCGCTGAAAGACATCGCCGTATCGTCCGGTTCGGCCTGTACCTCCGCGTCGCTCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTG‑GCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCaccacc                                                                                                                 >  1:46494/1‑149 (MQ=255)                                  ctgATGTCGCTGAAAGACATCGCCGTATCGTCCGGTTCGGCCTGTACCTCCGCGTCGCTCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTG‑GCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCaccacc                                                                                                                 <  1:171646/147‑1 (MQ=255)                                     aTGTCTCTGAAAGACATCGCCGTATCGTCCGGTTCGGCCTGTACCTCCGCGTCGCTCGAGCCGTCCTACGTACTTCGTGCTCTGGGCCGCAACGACGAGCTG‑GCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCgaa                                                                                                              <  4:167233/147‑1 (MQ=255)                                        tCGCTGAAAGACATCGCCGTATCGTCCGGTTCGGCCTGTACCTCCGCGTCGCTCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTG‑GCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCgaagaag                                                                                                          >  6:195132/1‑148 (MQ=255)                                          gCTGAAAGACATCGCCGTATCGTCCGGTTCGGCCTGTACCTCCGCGTCGCTCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTG‑GCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGaagaagaag                                                                                                       >  8:76346/1‑149 (MQ=255)                                           cTGAAAGACATCGCCGTATCGTCCGGTTCGGCCTGTACCTCCGCGTCGCTCGAGCAGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTG‑GCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGt                                                                                                      >  5:308075/1‑149 (MQ=255)                                            tGAAAGACATCGCCGTATCGTCCGGTTCGGCCTGTACCTCCGCGTCGCTCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTG‑GCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGt                                                                                                      >  4:196638/1‑148 (MQ=255)                                            tGAAAGACATCGCCGTATCGTCCGGTTCGGCCTGTACCTCCGCGTCGCTCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTG‑GCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTc                                                                                                     >  3:197733/1‑149 (MQ=255)                                             gAAAGACATCGCCGTATCGTCCGGTTCGGCCTGTACCTCCGCGTCGCTCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTG‑GCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTc                                                                                                     >  1:68953/1‑148 (MQ=255)                                              aaaGACATCGCCGTATCGTCCGGTTCGGCCTGTACCTCCGCGTCGCTCGAGCCGTCCTACGTACTGCGTGCTCTTGGCCGCAACGACGAGCTG‑GCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGa                                                                                                   <  3:196638/149‑1 (MQ=255)                                                          gTATCGTCCGGTTCGGCCTGTACCTCCGCGTCGCTCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTG‑GCGCTCAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCcgcg                                                                                         >  8:40584/1‑147 (MQ=255)                                                          gTATCGTCCGGTTCGGCCTGTACCTCCGCGTCGCTCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTG‑GCGCTCAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCa                                                                                       >  8:40586/1‑149 (MQ=255)                                                             tCGTCCGGTTCGGCCTGTACCTCCGCGTCGCTCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTG‑GCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGa                                                                                     <  3:145242/148‑1 (MQ=255)                                                                    gctCGGCCTGTACCTCCGCGTCGCGCGCGCCGTCCTACGTACTGCGTGCTCTGGCCCGCAACGACGAGCTG‑GCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGc                                                                             <  4:197733/147‑1 (MQ=255)                                                                       cGGCCTGTACCTCCGCGTCGCTCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTG‑GCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGCGAg                                                                          <  2:46494/149‑1 (MQ=255)                                                                       cGGCCTGTACCTCCGCGTCGCTCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTG‑GCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGCGAg                                                                          <  6:222630/149‑1 (MQ=255)                                                                        ggCCTGTACCTCCGCGTCGCTCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTG‑GCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCTCGCAGACTTTCTGCGCgg                                                                         <  2:264836/149‑1 (MQ=255)                                                                          ccTGTACCTCCGCGTCGCTCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTG‑GCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGCGAGGcc                                                                       >  1:305724/1‑149 (MQ=255)                                                                          ccTGTACCTCCGCGTCGCTCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTG‑GCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGCGAGGcc                                                                       >  6:60916/1‑149 (MQ=255)                                                                          ccTGTACCTCCGCGTCGCTCGAGCCGTACTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTG‑GCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGCGAGGcc                                                                       <  4:12494/149‑1 (MQ=255)                                                                                ccTCCGCGTCGCTCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTG‑GCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCt                                                                               >  2:263135/1‑135 (MQ=255)                                                                                ccTCCGCGTCGCTCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTG‑GCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCt                                                                               <  1:263135/135‑1 (MQ=255)                                                                                    cgcgTCGCTCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTG‑GCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGCGAGGCCGTCAATaaa                                                              <  7:76346/148‑1 (MQ=255)                                                                                    cgcgTCGCTCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTG‑GCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGCGAGGCCGTCAATaaa                                                              <  5:195132/148‑1 (MQ=255)                                                                                     gcgTCGCTCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTG‑GCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGCGAGGCCGTCAATAAACt                                                            <  7:156803/149‑1 (MQ=255)                                                                                        tCGCTCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTG‑GCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGCGAGGCCGTCAATAAACTGCg                                                         <  5:96691/149‑1 (MQ=255)                                                                                           cTCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTG‑GCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGCGAGGCCGTCAATAAACTGCGTGa                                                      <  6:308075/149‑1 (MQ=255)                                                                                           cTCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTG‑GCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGCGAGGCCGTCAATAAACTGCGTGa                                                      <  1:316620/149‑1 (MQ=255)                                                                                            tCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAAAGACGAGCTG‑GCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGCGAGGCCGTCAATAAACTGCGTGAg                                                     >  2:52650/1‑149 (MQ=255)                                                                                              gAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTG‑GCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGCGAGGCCGTCAATAAACTGCGTGAGCt                                                   <  7:32148/149‑1 (MQ=255)                                                                                              gAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTG‑GCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGCGAGGCCGTCAATAAACTGCGTGAGCt                                                   >  8:258535/1‑149 (MQ=255)                                                                                                      cTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTG‑GCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGCGAGGCCGTCAATAAACTGCGTGAGCTGTcgccgc                                           >  2:56213/1‑149 (MQ=255)                                                                                                      cTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTG‑GCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGCGAGGCCGTCAATAAACTGCGTGAGCTGTcgccgc                                           >  2:275492/1‑149 (MQ=255)                                                                                                      cTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTG‑GCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGCGAGGCCGTCAATAAACTGCGTGAGCTGTcgcagc                                           >  1:253883/1‑149 (MQ=255)                                                                                                        aCGTACTGCGTGCTCTGGGCCGCAACGACGAGCTG‑GCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGCGAGGCCGTCAATAAACTGCGTGAGCTGTCGCCGCtg                                         >  7:136971/1‑149 (MQ=255)                                                                                                          gTACTGCGTGCTCTGGGCCGCAACGACGAGCTG‑GCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGCGAGGcc                                                                       >  4:165673/1‑117 (MQ=255)                                                                                                          gTACTGCGTGCTCTGGGCCGCAACGACGAGCTG‑GCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGCGAGGcc                                                                       <  3:165673/117‑1 (MQ=255)                                                                                                            aCTGCGTGCTCTGGGCCGCAACGACGAGCTG‑GCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGCGAGGCCGTCAATAAACTGCGTGAGCTGTCGCCGCTGTgg                                      >  3:110741/1‑148 (MQ=255)                                                                                                                cGTGCTCTGGGCCGCAACGACGAGCTG‑GCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGCGAGGCCGTCAATAAACTGCGTGAGCTGTCGCCGCTGTGGGata                                  <  7:258535/148‑1 (MQ=255)                                                                                                                cGTGCTCTGGGCCGCAACGACGAGCTG‑GCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCACAGAAAGTCTGCGAGGCCGTCAATAAACTGCGTGAGCTGTCGCCGCTGTGGGatat                                 >  1:178373/1‑149 (MQ=255)                                                                                                                 gTGCTCTGGGCCGCAACGACGAGCTG‑GCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGCGAGGCCGTCAATAAACTGCGTGAGCTGTCGCCGCTGTGGGATATg                                >  4:244922/1‑149 (MQ=255)                                                                                                                 gTGCTCTGGGCCGCAACGACGAGCTG‑GCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGCGAGGCCGTCAATAAACTGCGTGAGCTGTCGCCGCTGTGGGATATg                                >  1:34967/1‑149 (MQ=255)                                                                                                                  tGCTCTGGGCCGCAACGACGAGCTG‑GCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGCGAGGCCGTCAATAAACTGCGTGAGCTGTCGCCGCTGTGGGATATg                                >  7:129159/1‑148 (MQ=255)                                                                                                                     tctGGGCCGCAACGACGAGCTG‑GCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCG‑AGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGCGAGGCCGTCAATAAACTGCGTGAGCTGTCGCCGCTGTGGGATATGTACa                            >  1:190735/1‑148 (MQ=255)                                                                                                                     tctGGGCCGCAAAGACGAGCTT‑GCGCACAGCTCGATCCGCTTATCCTACGGCCGATTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTTCGAGGCCCTCAATAAAATGCCTGATCTGTCGCCGCTGTGGGATATGTAc                             >  4:100019/1‑148 (MQ=255)                                                                                                                            cGCAACGACGAGCTG‑GCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGCGAGGCCGTCAATAAACTGCGTGAGCTGTCGCCGCTGTGGGATATGTACAAAGACgg                     >  8:165503/1‑149 (MQ=255)                                                                                                                             gCAACGACGAGCTG‑GCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCACAGAAAGTCTGCGAGGCCGTCAATAAACTGCGTGAGCTGTCGCCGCTGTGGGATATGTACAAAGACGGc                    <  2:178373/149‑1 (MQ=255)                                                                                                                              cAACGACGAGCTG‑GCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGCGAGGCCGTCAATAAACTGCGTGAGCTGTCGCCGCTGTGGGATATGTACAAAGACGGCg                   >  4:301406/1‑149 (MQ=255)                                                                                                                                ggatcgaGCTGTGCGC‑CAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGCGAGGCCGTCAATAAACTGCGTGAGCTGTCGCCGCTGTGGGATATGTACAAAGACGGCGtt                 >  2:245591/5‑149 (MQ=255)                                                                                                                                acgacgAGCTG‑GCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGCGAGGCCGTCAATAAACTGCGTGAGCTGTCGCCGCTGTGGGATATGTACAAAGACGGCGtt                 <  2:68953/149‑1 (MQ=255)                                                                                                                                 gatcgaGCTGTGCGC‑CAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGCGAGGCCGTCAATAAACTGCGTGAGCTGTCGCCGCTGTGGGATATGTACAAAGACGGCGtt                 >  7:183981/4‑148 (MQ=255)                                                                                                                                  gacgaGCTG‑GCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGCGAGGCCGTCAATAAACTGCGTGAGCTGTCGCCGCTGTGGGATATGTACAAAGACGGCGTTGa               <  6:122452/149‑1 (MQ=255)                                                                                                                                   acgaGCTG‑GCGCTCAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGCGAGGCCGTCAATAAACTGCGTGAGCTGTCGCCGCTGTGGGATATGTACAAAGACGGCGTTGAc              <  7:40584/149‑1 (MQ=255)                                                                                                                                   acgaGCTG‑GCGCTCAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGCGAGGCCGTCAATAAACTGCGTGAGCTGTCGCCGCTGTGGGATATGTACAAAGACGGCGTTGAc              <  7:40586/149‑1 (MQ=255)                                                                                                                                   acgaGCTG‑GCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGCGAGGCCGTCAATAAACTGCGTGAGCTGTCGCCGCTGTGGGATATGTACAAAGACGGCGTTGAc              <  2:305724/149‑1 (MQ=255)                                                                                                                                       gCTG‑GCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGCGAGGCCGTCAATAAACTGCGTGAGCTGTCGCCTCTGTGGGATATGTACAAAGACGGCGTTGACAtct          >  1:176251/1‑149 (MQ=255)                                                                                                                                            gcgcACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGCGAGGCCGTCAATAAACTGCGTGAGCTGTCGCCGCTGTGGGATATGTACAAAGACGGCGTTGACATCTCCaa      >  6:303717/1‑149 (MQ=255)                                                                                                                                             cgcACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGCGAGGCCGTCAATAAACTGCGTGAGCTGTCGCCGCTGTGGGATATGTACAAAGACGGCGTTGACATCTcc        >  5:261122/1‑146 (MQ=255)                                                                                                                                              gcACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGCGAGGCCGTCAATAAACTGCGTGAGCTGTCGCCGCTGTGGGATATGTACAAAGACGGAGTTGACATCTCCAAg     >  8:80530/1‑148 (MQ=255)                                                                                                                                                acaGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGCGAGGCCGTCAATAAACTGCGTGAGCTGTCGCCGCTGTGGGATATGTACAAAGACGGCGTTGACATCTCCAAGATc  <  6:261122/149‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                     CTGAGCTTCAACTACGTCGAAGGCGAATCGCTGCTGATGTCGCTGAAAGACATCGCCGTATCGTCCGGTTCGGCCTGTACCTCCGCGTCGCTCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTG‑GCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGCGAGGCCGTCAATAAACTGCGTGAGCTGTCGCCGCTGTGGGATATGTACAAAGACGGCGTTGACATCTCCAAGATC  >  NC_002947/979457‑979747

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