Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,004,143	T→C	8.9%	W37R (TGG→CGG) ‡	PP_3533 →	ornithine cyclodeaminase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,004,143	0	T	C	8.9%	184.7 / 7.1	67	W37R (TGG→CGG) ‡	PP_3533	ornithine cyclodeaminase
Reads supporting (aligned to +/- strand):  ref base T (26/35);  new base C (3/3);  total (29/38)
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.00e-01

GCAAGCGCACAATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCTGCATACGG  >  NC_002947/4004006‑4004287                                                                                                                                           |                                                                                                                                                 gCAAGCGCACAATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGttt                                                                                                                                       <  1:16734/149‑1 (MQ=255) gCAAGCGCACAATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGttt                                                                                                                                       >  7:150228/1‑149 (MQ=255)  cAAGCGCACAATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTACTTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGttt                                                                                                                                       <  1:85059/149‑1 (MQ=255)  cAAGCGCACAATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGc                                                                                                                                               >  2:275363/1‑140 (MQ=255)  cAAGCGCACAATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGc                                                                                                                                               <  1:275363/140‑1 (MQ=255)   aaGCGCACAATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGa                                                                                                                                     <  4:293773/149‑1 (MQ=255)     gcgcACAATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCCTGCCAGCGTTTGaaa                                                                                                                                   <  2:12881/149‑3 (MQ=255)     gcgcACAATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCCTGCCAGCGTTTGaaa                                                                                                                                   <  1:301433/149‑3 (MQ=255)      cgcACAATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACaa                                                                                                                                  <  8:150228/149‑1 (MQ=255)      cgcACAATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACaa                                                                                                                                  <  2:207734/149‑1 (MQ=255)        cacaATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGt                                                                                                                                >  6:322254/1‑149 (MQ=255)                    gTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTg                                                                                                                    >  1:106050/1‑149 (MQ=255)                       cAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTgcca                                                                                                                 <  2:95596/149‑1 (MQ=255)                        aaGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGAc                                                                                                                                    <  1:307848/129‑1 (MQ=255)                        aaGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGAc                                                                                                                                    >  2:307848/1‑129 (MQ=255)                        aaGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTgccag                                                                                                                <  8:221609/149‑1 (MQ=255)                        aaGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTgccag                                                                                                                >  6:156467/1‑149 (MQ=255)                                gACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCg                                                                                                         >  8:71638/1‑148 (MQ=255)                                   gTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAg                                                                                                     <  6:321128/149‑1 (MQ=255)                                          aTTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGtt                                                                                                                     <  4:12482/126‑1 (MQ=255)                                          aTTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGtt                                                                                                                     >  3:12482/1‑126 (MQ=255)                                          aTTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCg                                                                                               <  2:295357/148‑1 (MQ=255)                                            tGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTgccagcca                                                                                                             >  3:65223/1‑132 (MQ=255)                                            tGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTgccagcca                                                                                                             <  4:65223/132‑1 (MQ=255)                                              aTGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTc                                                                                                          >  3:215014/1‑133 (MQ=255)                                              aTGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTc                                                                                                          <  4:215014/133‑1 (MQ=255)                                                   ccAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAAc                                                                                      <  4:187883/148‑1 (MQ=255)                                                   ccAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCCTGCCAGCGTTTg                                                                                                                                      >  6:302900/1‑100 (MQ=255)                                                   ccAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCCTGCCAGCGTTTg                                                                                                                                      <  5:302900/100‑1 (MQ=255)                                                      aCCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGAt                                                                                  >  6:207431/1‑149 (MQ=255)                                                       ccATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATg                                                                                 <  2:106050/149‑1 (MQ=255)                                                       ccATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATg                                                                                 <  5:207431/149‑1 (MQ=255)                                                        cATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGc                                                                                <  5:156467/149‑1 (MQ=255)                                                        cATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGc                                                                                >  4:59228/1‑149 (MQ=255)                                                          tGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCg                                                                              <  5:308269/149‑1 (MQ=255)                                                           gTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCgg                                                                             >  4:49201/1‑149 (MQ=255)                                                              ggACTTGGTGCATGACATTGGCGTAGCGCCGTTTA‑TTGGCTAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTc                                                                           <  3:39472/148‑1 (MQ=255)                                                                aCTTGGTGCATGACATTGGCGTAGCGCCGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGcc                                                                        <  7:71638/149‑1 (MQ=255)                                                                        cATGACATTGGCGTAGCGCCGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACaaaa                                                                 <  8:186175/148‑1 (MQ=255)                                                                            aCATTGGCGTAGCGCCGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCg                                                             <  7:185753/148‑1 (MQ=255)                                                                                          ccGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGATAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATa                                              >  6:217262/1‑149 (MQ=255)                                                                                             tttA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACaaaa                                                                 <  5:233956/127‑1 (MQ=255)                                                                                             tttA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACaaaa                                                                 >  6:233956/1‑127 (MQ=255)                                                                                             tttA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGt                                           <  3:49201/149‑1 (MQ=255)                                                                                                  ttGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTAGGCCTTCAAATACGTCAAc                                       >  6:224453/1‑149 (MQ=255)                                                                                                     gCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGc                                    >  3:171769/1‑149 (MQ=255)                                                                                                      cGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGc                                    >  4:187048/1‑148 (MQ=255)                                                                                                          cTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCaaa                                                >  6:279060/1‑132 (MQ=255)                                                                                                          cTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCAccc                               >  8:169540/1‑149 (MQ=255)                                                                                                          cTTGCCGCTGCCCTGCGGGACGATATCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCaaa                                                <  5:279060/132‑1 (MQ=255)                                                                                                           ttGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCg                              >  5:131720/1‑149 (MQ=255)                                                                                                               cGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGcc                           >  2:238094/1‑148 (MQ=255)                                                                                                                 cTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCaa                         >  8:298560/1‑148 (MQ=255)                                                                                                                    cccTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATa                       <  3:212827/147‑1 (MQ=255)                                                                                                                         cGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACg                <  3:59228/149‑1 (MQ=255)                                                                                                                             aCGATTTCAAACGCTGGCAGGCGTTTGATAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAAc            <  5:217262/149‑1 (MQ=255)                                                                                                                             aCGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAAc            <  4:171769/149‑1 (MQ=255)                                                                                                                             aCGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGCGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAAc            <  7:169540/149‑1 (MQ=255)                                                                                                                                 tttCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCTg         >  7:50133/1‑148 (MQ=255)                                                                                                                                 tttCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCTg         <  8:50133/148‑1 (MQ=255)                                                                                                                                 tttCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCTGc        <  1:238094/149‑1 (MQ=255)                                                                                                                                  gtCAAACGCCTGCCAGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCTGCa       >  5:44534/2‑149 (MQ=255)                                                                                                                                   tCAAACGCTGGCAGGCGTTTGACAATTCCTCGCGCGATGCCAGCCAATAGTAAGTGTGAGTGATCGAACTTATGCCGGTCGCCGAAAAAATCCGCTATGCCTTCAAATACGTCAACGGAAACCAGGCCTATACTGAACTCAAcattca       >  3:264729/1‑143 (MQ=255)                                                                                                                                   tCAAACGCCTGCCAGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCTGCa       >  2:172906/1‑148 (MQ=255)                                                                                                                                    cAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGc                                                                                >  6:40474/1‑74 (MQ=255)                                                                                                                                    cAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGc                                                                                <  5:40474/74‑1 (MQ=255)                                                                                                                                       aCGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCTGCATACgg  <  3:187048/149‑1 (MQ=255)                                                                                                                                           |                                                                                                                                                 GCAAGCGCACAATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCTGCATACGG  >  NC_002947/4004006‑4004287

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