Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,638,690	C→T	12.5%	G171G (GGC→GGT)	creD →	inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,638,690	0	C	T	12.5%	194.2 / 18.0	72	G171G (GGC→GGT)	creD	inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (25/38);  new base T (2/7);  total (27/45)
Fisher's exact test for biased strand distribution p-value = 4.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.11e-01

CGATGTTTCGCGTCTTAGCGAACTCAACGCGCCAAATATCACCTTAGGCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGA  >  NC_000913/4638546‑4638838                                                                                                                                                 |                                                                                                                                                     cGATGTTTCGCGTCTTAGCGAACTCAACGCGCCAAATATCACCTTAGGCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCg                                                                                                                                                  <  3:101750/149‑1 (MQ=255)    tGTTTCGCGTCTTAGCGAACTCAACGCGCCAAATATCACCTTAGGCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGt                                                                                                                                               <  3:231450/149‑1 (MQ=255)        tCGCGTCTTAGCGAACTCAACGCGCCAAATATCACCTTAGGCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTaga                                                                                                                                           <  1:384032/149‑1 (MQ=255)         cgcgTCTTAGCGAACTCAACGCGCCAAATATCACCTTAGGCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGTGCCGGGttcaat                                                                                                                                          <  2:678209/149‑5 (MQ=255)         cgcgTCTTAGCGAACTCAACGCGCCAAATATCACCTTAGGCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGTGCCGGGttcaat                                                                                                                                          <  4:435517/149‑5 (MQ=255)         cgcgTCTTAGCGAACTCAACGCGCCAAATATCACCTTAGGCAAGCCATCTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGTGCCGGGttcaat                                                                                                                                          <  3:624290/149‑5 (MQ=255)          gcgTCTTAGCGAACTCAACGCGCCAAATATCACCTTAGGCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGTGCCGGGttcaatg                                                                                                                                         <  4:679141/149‑6 (MQ=255)          gcgTCTTAGCGAACTCAACGCGCCAAATATCACCTTAGGCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGTGCCGGGttcaatg                                                                                                                                         <  1:334185/149‑6 (MQ=255)          gcgTCTTAGCGAACTCAACGCGCCAAATATCACCTTAGGCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGc                                                                                                                                         >  3:134839/1‑149 (MQ=255)               ttAGCGAACTCAACGCGCCAAATATCACCTTAGGCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGTGCCGGGttcaa                                                                                                                                           <  1:361387/142‑4 (MQ=255)               ttAGCGAACTCAACGCGCCAAATATCACCTTAGGCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGTGCCGGGttcaa                                                                                                                                           >  2:361387/1‑139 (MQ=255)               ttAGCGAACTCAACGCGCCAAATATCACCTTAGGCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGc                                                                                                                                    >  3:398305/1‑149 (MQ=255)                tAGCGAACTCAACGCGCCAAATATCACCTTAGGCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCg                                                                                                                                   <  2:88734/149‑1 (MQ=255)                tAGCGAACTCAACGCGCCAAATATCACCTTAGGCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCg                                                                                                                                   <  2:702537/149‑1 (MQ=255)                tAGCGAACTCAACGCGCCAAATATCACCTTAGGCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCg                                                                                                                                   <  2:142746/149‑1 (MQ=255)                  gCGAACTCAACGCGCCAAATATCACCTTAGGCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTaga                                                                                                                                           <  1:679287/139‑1 (MQ=255)                  gCGAACTCAACGCGCCAAATATCACCTTAGGCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTaga                                                                                                                                           >  2:679287/1‑139 (MQ=255)                          aaCGCGCCAAATATCACCTTAGGCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTaga                                                                                                                                           >  3:648736/1‑131 (MQ=255)                          aaCGCGCCAAATATCACCTTAGGCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTaga                                                                                                                                           <  4:648736/131‑1 (MQ=255)                                 cAAATATCACCTTAGGCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATc                                                                                                                  <  1:246483/149‑1 (MQ=255)                                  aaaTATCACCTTAGGCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGCGTTGTATTGGTGTGGTGAAAGCGCCTGAAGTGAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATcc                                                                                                                 <  1:414993/149‑1 (MQ=255)                                    atatCACCTTAGGCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCt                                                                                                               <  1:543478/149‑1 (MQ=255)                                               ggCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAgg                                                                                                    <  4:398305/149‑1 (MQ=255)                                                gCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAggg                                                                                                   <  1:54475/149‑1 (MQ=255)                                                gCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAggg                                                                                                   >  4:540883/1‑149 (MQ=255)                                                gCAAGCCATTGATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAggg                                                                                                   >  4:493385/1‑149 (MQ=255)                                                  aaGCCATTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGa                                                                                                 >  2:577084/1‑149 (MQ=255)                                                  aaGCCATTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGa                                                                                                 <  1:202940/149‑1 (MQ=255)                                                     ccATTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTg                                                                                              >  4:362201/1‑149 (MQ=255)                                                     ccATTTATTCTGATTAGCGTCGGGGATGCGCGTGGGATTGGTGTGGTGAAAGCGCCTGAAGATAACGGTATGGCGCTGACCTTTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGTCGTGCATATCCCGTTACCTGAAGCGGTCTg                                                                                              >  4:350483/1‑149 (MQ=255)                                                       aTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGc                                                                                            <  3:42748/149‑1 (MQ=255)                                                        tttATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCg                                                                                           <  3:401799/149‑1 (MQ=255)                                                           aTTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGaa                                                                                        <  1:577084/149‑1 (MQ=255)                                                           aTTGTGATTAGCGTCGGGGATGCGCCAGGTATTGGTCTGGTGATAGCGCCTGAAGTTAATGGATCGGCGGTGACCATTGATCCCGGCACCGGGATAGAGCAAGGCGGGCAGGGCGTGCAGATCCATTTACCTGAAGGGGACTGGCGaaa                                                                                        <  4:470841/149‑1 (MQ=255)                                                            ttGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTaga                                                                                                                                           <  1:132051/97‑1 (MQ=255)                                                            ttGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTaga                                                                                                                                           >  2:132051/1‑97 (MQ=255)                                                            ttGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAg                                                                                       >  4:143533/1‑149 (MQ=255)                                                            ttGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAg                                                                                       <  1:121412/149‑1 (MQ=255)                                                                 aTTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGaa                                                                                  >  2:481433/1‑149 (MQ=255)                                                                       gTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAGCCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTTCCTGAAGGGGACTGGCGGAAGCAGAACCTGaa                                                                            <  4:170918/149‑1 (MQ=255)                                                                          ggggATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCt                                                                         >  4:643491/1‑149 (MQ=255)                                                                             gATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGaa                                                                      <  3:355304/149‑1 (MQ=255)                                                                                  gcgTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATgg                                                                 >  1:228915/1‑149 (MQ=255)                                                                                  gcgTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATgg                                                                 <  3:673678/149‑1 (MQ=255)                                                                                  gcgTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATgg                                                                 <  4:294199/149‑1 (MQ=255)                                                                                   cgTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGc                                                                >  2:587287/1‑149 (MQ=255)                                                                                   cgTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGc                                                                <  4:134839/149‑1 (MQ=255)                                                                                          ttGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAAt                                                         >  4:609813/1‑149 (MQ=255)                                                                                          ttGGTGTCGTGAAAGCGCCAGAAGTTAACGGAACGGCGCAGACCATTGAACCCGGCACCGGGTTAGAGCAACGCCGGCAGGGCGTGCATATCTCTATTCCTCAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATAAGGCCGTGAAt                                                         <  2:235912/149‑1 (MQ=255)                                                                                              tgtggtgAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTaa                                                     >  3:595061/1‑149 (MQ=255)                                                                                              tgtggtgAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTaa                                                     <  1:648950/149‑1 (MQ=255)                                                                                              tgtggtgAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTaa                                                     <  3:698190/149‑1 (MQ=255)                                                                                                tggtgAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGc                                                   <  1:504278/149‑1 (MQ=255)                                                                                                tggtgAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGc                                                   <  4:198148/149‑1 (MQ=255)                                                                                                  gtgAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAgg                                                                                                    <  1:667146/98‑1 (MQ=255)                                                                                                  gtgAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAgg                                                                                                    >  2:667146/1‑98 (MQ=255)                                                                                                  gtgAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCgg                                                 >  4:95652/1‑149 (MQ=255)                                                                                                          gcCTGAATTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCg                                         >  4:15424/1‑149 (MQ=255)                                                                                                                  ttAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCt                                 >  2:170841/1‑149 (MQ=255)                                                                                                                    aacggaacggCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCtgt                               <  1:170841/149‑1 (MQ=255)                                                                                                                       ggaacggCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTgtggt                            >  3:632612/1‑149 (MQ=255)                                                                                                                       ggaacggCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTgtggt                            >  4:586701/1‑149 (MQ=255)                                                                                                                          acggCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGcc                         <  1:10542/149‑1 (MQ=255)                                                                                                                                 tGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTggcggg                  <  1:451608/149‑1 (MQ=255)                                                                                                                                      aTTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTaa             >  2:145120/1‑149 (MQ=255)                                                                                                                                      aTTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTaa             <  4:593123/149‑1 (MQ=255)                                                                                                                                       ttGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAAt            >  2:264490/1‑149 (MQ=255)                                                                                                                                       tctaaCCCGGTGCCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAg                                                                                                     >  4:484001/4‑60 (MQ=38)                                                                                                                                       tctaaCCCGGTGCCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAg                                                                                                     <  3:484001/57‑1 (MQ=255)                                                                                                                                         gAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAg          <  1:587287/149‑1 (MQ=255)                                                                                                                                          aaCCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGc         <  2:493540/149‑1 (MQ=255)                                                                                                                                          aaCCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGc         <  3:35552/149‑1 (MQ=255)                                                                                                                                                 cACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGa  >  2:657580/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                     CGATGTTTCGCGTCTTAGCGAACTCAACGCGCCAAATATCACCTTAGGCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGA  >  NC_000913/4638546‑4638838

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑