Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,338,916	C→T	V536I (GTC→ATC)	adiA ←	arginine decarboxylase, degradative

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,338,916	0	C	T	97.2%	129.5 / NA	36	V536I (GTC→ATC)	adiA	arginine decarboxylase, degradative
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (17/18);  minor base G (0/1);  total (17/19)
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.98e-01

CACCTGCGCCAGCGGTGTGTTGGCGTCATAGTGGCGTTTGAAGGAGCAAAGGGTGTTAACCAGAGTTCCCCATTTCCCACGGGTTACGCCCATAGAGAACAGGAACATAATTTGGAAGTCAGTGGTGCGGGTAGGTACAATGCCGTGGCGACCAAGCCAGGCAGTGACCAGCGCCGCCGGAACACCGGTTTCTTCCAGTTCACCATCTTCACCCATTCCCGGAGCAAGGATGCTGACTTTAATCGGGTCGAGCATACTCCAGTTATCCGGAATATCTTTGAAGCCGTGCCAGCTTTCGCCCGGATGCATTACCCAGCAGTCCTGAACGGTGGTCAGCAGTTTGGTTGGTGCGTCAGCAAAGTCATAGGTTTTGCCGGTTTGTGGGTCGGTGACGACTTCTTTGTTCCACGGTTTGAAGAACCAGCTACCGTCAGCGGTGAACTCTTTATATAGCCGCGCCATCGCCTGACGGAAATC  >  NC_000913/4338680‑4339156                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                 cACCTGCGCCAGCGGTGTGTTGGCGTCATAGTGGCGTTTGAAGGAGCAAAGGGTGTTAACCAGAGTTCCCCATTTCCCACGGGTTACGCCCATAGAGAACAGGAACATAATTTGGAAGTCAGTGGTGCGGGTAGGTACAATGCCGTGGCGACCAAGCCAGGCAGTGACCAGCGCCGCCGGAACACCGGTTTCTTCCAGTTCACCATCTTCACCCATTCCCGGAGCAAGGATGCTGATTTTAATCGGGTCg                                                                                                                                                                                                                                     <  2:120506/250‑1 (MQ=255)   ccTGCGCCAGCGGTGTGTTGGCGTCATAGTGGCGTTTGAAGGAGCAAAGGGTGTTAACCAGAGTTCCCCATTTCCCACGGGTTACGCCCATAGAGAACAGGAACATAATTTGGAAGTCAGTGGTGCGGGTAGGTACAATGCCGTGGCGACCAAGCCAGGCAGTGACCAGCGCCGCCGGAACACCGGTTTCTTCCAGGTCACCATCTTCACCCATTCCCGGAGCAAGGATGCTGATTGTAATCGGTTCGAGc                                                                                                                                                                                                                                  <  2:168294/251‑1 (MQ=255)       cgcCAGCGGTGTGTTGGCGTCATAGTGGCGTTTGAAGGAGCAAAGGGTGTTAACCAGAGTTCCCCATTTCCCACGGGTTACGCCCATAGAGAACAGGAACATAATTTGGAAGTCAGTGGTGCGGGTAGGTACAATGCCGTGGCGACCAAGCCAGGCAGTGACCAGCGCCGCCGGAACACCGGTTTCTTCCAGTTCACCATCTTCACCCATTCCCGGAGCAAGGATGCTGATTTTAATCGGGTCGAGCATAc                                                                                                                                                                                                                              >  1:157856/1‑251 (MQ=255)                            aTAGTGGCGTTTGAAGGAGCAAAGGGTGTTAACCAGAGTTCCCCATTTCCCACGGGTTACGCCCATAGAGAACAGGAACATAATTTGGAAGTCAGTGGTGCGGGTAGGTACAATGCCGTGGCGACCAAGCCAGGCAGTGACCAGCGCCGCCGGAACACCGGTTTCTTCCAGTTCACCATCTTCACCCATTCCCGGAGCAAGGATGCTGATTTTAATCGGGTCGAGCATACTCCAGTTATCCGGAATATCtt                                                                                                                                                                                                         >  2:295185/1‑251 (MQ=255)                                         aaGGAGCAAAGGGTGTTAACCAGAGTTCCCCATTTCCCACGGGTTACGCCCATAGAGAACAGGAACATAATTTGGAAGTCAGTGGTGCGGGTAGGTACAATGCCGTGGCGACCAAGCCAGGCAGTGACCAGCGCCGCCGGAACACCGGTTTCTTCCAGTTCACCATCTTCACCCATTCCCGGAGCAAGGATGCTGATTTTAATCGGGTCGAGCATACt                                                                                                                                                                                                                             >  1:252935/1‑218 (MQ=255)                                         aaGGAGCAAAGGGTGTTAACCAGAGTTCCCCATTTCCCACGGGTTACGCCCATAGAGAACAGGAACATAATTTGGAAGTCAGTGGTGCGGGTAAGTACAATGCCGTGGCGACCAAGCCAGGCAGTGACCAGCGCCGCCGGAACACCGGTTTCTTCCAGTTCACCATCTTCACCCATTCCCGGAGCAAGGATGCTGATTTTAATCGGGTCGAGCATACt                                                                                                                                                                                                                             <  2:252935/218‑1 (MQ=255)                                           ggAGCAAAGGGTGTTAACCAGAGTTCCCCATTTCCCACGGGTTACGCACATAGAGAACAGGAACATAATTTGGAAGTCAGTGGTGCGGGTAGGTACAATGCCGTGGCGACCAAGCCAGGCAGTGACCAGCGCCGCCGGAACACCGGTTTCTTCCAGTTCACCATCTTCACCCATTCCCGGAGCAAGGATGCTGATTTTAATCGGGTCGAGCATACTCCAGTTATCCGGAATATCTTTGAAGCCGTGCCAGc                                                                                                                                                                                          >  2:145700/1‑251 (MQ=255)                                              gCAAAGGGTGGTAACCAGAGTTCCCCATGTCCCACGGGTTACGCCCATAGAGAACAGGAACAGAATTTGGAAGTCAGTGGTGCGGGTAGTTACAATGCCGTGGCGACCAAGCCAGGCAGTGACCAGCGCCGCCGGAACACCGGTTTCTTCCAGTTCACCATCTTCACCCATTCCCGGAGCAAGGATGCTGATTTTAATCGGGTCGAGCATACTCCAGGTAGCCGGAATGTCTTTGAAGCCGTGCCAGCttt                                                                                                                                                                                       <  2:228378/251‑1 (MQ=255)                                                    ggTGTTAACCAGAGTTCCCCATTTCCCACGGGTTACGCCCATAGAGAACAGGAACATAATTTGGAAGTCAGTGGTGCGGGTAGGTACAATGCCGTGGCGACCAAGCCAGGCAGTGACCAGCGCCGCCGGAACACCGGTTTCTTCCAGTTCACCATCTTCACCCATTCCCGGAGCAAGGATGCTGATTTTAATCGGGTCGAGCATACTCCAGTTATCCGGAATATCTTTGAAGCCGTGCCAGCTTTCGCCCg                                                                                                                                                                                 <  1:92105/251‑1 (MQ=255)                                                      tgtTAACCAGAGTTCCCCATTTCCCACGGGTTACGCCCATAGAGAACAGGAACATAATTTGGAAGTCAGTGGTGCGGGTAGGTACAATGCCGTGGCGACCAAGCCAGGCAGTGACCAGCGCCGCCGGAACACCGGTTTCTTCCAGTTCACCATCTTCACCCATTCCCGGAGCAAGGATGCTGATTTTAATCGGGTCGAGCATACTCCAGTTATCCGGAATATCTTTGAAGCCGTGCCAGCTTTCGCCCgg                                                                                                                                                                                >  2:216033/1‑250 (MQ=255)                                                       gtTAACCAGAGTTCCCCATTTCCCACGGGTTACGCCCATAGAGAACAGGAACATAATTTGGAAGTCAGTGGTGCGGGTAGGTACAATGCCGTGGCGACCAAGCCAGGCAGTGACCAGCGCCGCCGGAACACCGGTTTCTTCCAGTTCACCATCTTCACCCATTCCCGGAGCAAGGATGCTGATTTTAATCGGGTCGAGCATACTCCAGTTATCCGGAATATCTTTGAAGCCGTGCCAGCTTTCGCCCGGAt                                                                                                                                                                              >  1:323102/1‑251 (MQ=255)                                                           aCCAGAGTTCCCCATTTCCCACGGGTTACGCCCATAGAGAACAGGAACATAATTTGGAAGTCAGTGGTGCGGGTAGGTACAATGCCGTGGCGACCAAGCCAGGCAGTGACCAGCGCCGCCGGAACACCGGTTTCTTCCAGTTCACCATCTTCACCCATTCCCGGAGCAAGGATGCTGATTTTAATCGGGTCGAGCATACTCCAGTTATCCGGAATATCTTTGAAGCCGTGCCAGCTTTCGCCCGGATGCAt                                                                                                                                                                          <  1:288497/251‑1 (MQ=255)                                                                              cACGGGTTACGCCCATAGAGAACAGGAACATAATTTGGAAGTCAGTGGTGCGGGTAGGTACAATGCCGTGGCGACCAAGCCAGGCAGTGACCAGCGCCGCCGGAACACCGGTTTCTTCCAGTTCACCATCTTCACCCATTCCCGGAGCAAGGATGCTGATTTTAATCGGGTCGAGCATACTCCAGTTATCCGGAATATCTTTGAAGCCGTGCCAGCTTTCGCCCGGATGCATTACCCAGCAGTCCTGAACg                                                                                                                                                       >  1:115650/1‑251 (MQ=255)                                                                                    ttACGCCCATAGAGAACAGGAACATAATTTGGAAGTCAGTGGTGCGGGTAGGTACAATGCCGTGGCGACCAAGCCAGGCAGTGACCAGCGCCGCCGGAACACCGGTTTCTTCCAGTTCACCATCTTCACCCATTCCCGGAGCAAGGATGCTGATTTTAATCGGGTCGAGCATACTCCAGTTATCCGGAATATCTTTGAAGCCGTGCCAGCTTTCGCCCGGCTGCATTACCCAGCAGTCCTGAACGGTGGTc                                                                                                                                                 >  2:188539/1‑251 (MQ=255)                                                                                              agagAACAGGAACATAATTTGGAAGTCAGTGGTGCGGGTAGGTACAATGCCGTGGCGACCAAGCCAGGCAGTGACCAGCGCCGCCGGAACACCGGTTTCTTCCAGTTCACCATCTTCACCCATTCCCGGAGCAAGGATGCTGAGTTTAATCGGGTCGAGCATACTCCAGTTATCCGGAACATCTTTGAAGCCGTGCCAGCTTTCGCCCGGATGCATTACCCAGCAGTCCTGAACGGTGGTCAGCAGTttg                                                                                                                                        <  2:282872/250‑1 (MQ=255)                                                                                                      ggAACATAATTTGGAAGTCAGTGGTGCGGGTAGGTACAATGCCGTGGCGACCAAGCCAGGCAGTGACCAGCGCCGCCGGAACACCGGTTTCTTCCAGTTCACCATCTTCACCCATTCCCGGAGCAAGGATGCTGATTTTAATCGGGTCGAGCATACTCCAGTTATCCGGAATATCTTTGAAGCCGTGCCAGCTTTCGCCCGGATGCATTACCCAGCAGTCCTGAACGGTGGTCAGCAGTTTGGTTGGTGCg                                                                                                                               >  2:80931/1‑251 (MQ=255)                                                                                                         aCATAATTTGGAAGTCAGTGGTGCGGGTAGGTACAATGCCGTGGCGACCAAGCCAGGCAGTGACCAGCGCCGCCGGAACACCGGTTTCTTCCAGTTCACCATCTTCACCCATTCCCGGAGCAAGGATGCTGATTTTAATCGGGTCGAGCATACTCCAGTTATCCGGAATATCTTTGAAGCCGTGCCAGCTTTCGCCCGGATGCATTACCCAGCAGTCCTGAACGGTGGTCAGCAGTTTGGTTGGTGCGTCa                                                                                                                            <  1:351369/251‑1 (MQ=255)                                                                                                                ttGGAAGTCAGTGGTCCGGGTAGGTACAATGCCGTGGCGACCAAGCCAGGCAGTGACCCGCGCCGCCGGAACACCGGTTTCTTCCAGTTCACCATCTTCACCCATTCCCGGAGCAAGGATGCTGATTTTAATCGGGTCGAGCATACTCCAGTTATCCGGAATATCTTTGAAGCCGTGCCAGCTTTCGCCCGGATGCATTACCCAGCAGTCATGAACGGTGGTCAGCAGTTTGGTTGGTGCGTCAGCAAAGt                                                                                                                     >  1:46472/1‑251 (MQ=255)                                                                                                                 tGGAAGTCAGTGGTGCGGGTAGGTACAATGCCGTGGCGACCAAGCCAGGCAGTGACCAGCGCCGCCGGAACACCGGTTTCTTCCAGTTCACCATCTTCACCCATTCCCGGAGCAAGGATGCTGATTTTAATCGGGTCGAGCATACTCCAGTTATCCGGAATATCTTTGAAGCCGTGCCAGCTTTCGCCCGGATGCATTACCCAGCAGTCCTGAACGGTGGTCAGCAGTTTGGTTGGTGCGTCAGCAAAGt                                                                                                                     <  2:24729/250‑1 (MQ=255)                                                                                                                         aGTGGTGCGGGTAGGTACAATGCCGTGGCGACCAAGCCAGGCAGTGACCAGCGCCGCCGGAACACCGGTTTCTTCCAGTTCACCATCTTCACCCATTCCCGGAGCAAGGATGCTGATTTTAATCGGGTCGAGCATACTCCAGTTATCCGGAATATCTTTGAAGCCGTGCCAGCTTTCGCCCGGATGCATTACCCAGCAGTCCTGAACGGTGGTCAGCAGTTTGGTTGGTGCGTCAGCAAAGTCATAGGttt                                                                                                            >  2:22418/1‑251 (MQ=255)                                                                                                                                          cAATGCCGTGGCGACCAAGCCAGGCAGTGACCAGCGCCGCCGGAACACCGGTTTCTTCCAGTTCACCATCTTCACCCATTCCCGGAGCAAGGATGCTGATTTTAATCGGGTCGAGCATACTCCAGTTATCCGGAATATCTTTGAAGCCGTGCCAGCTTTCGCCCGGATGCATTACCCAGCAGTCCTGAACGGTGGTCAGCAGTTTGGTTGGTGCGTCAGCAAAGTCATAGGTTTTGCCGGTTTGTGGGTCg                                                                                           <  2:323102/251‑1 (MQ=255)                                                                                                                                          cAATGCCGTGGCGACCAAGCCAGGCAGTGACCAGCGCCGCCGGAACACCGGTTTCTTCCAGTTCACCATCTTCACCCATTCCCGGAGCAAGGATGCTGATTTTAATCGGGTCGAGCATACTCCAGTTATCCGGAATATCTTTGAAGCCGTGCCAGCTTTCGCCCGGATGCATTACCCAGCAGTCCTGAACGGTGGTCAGCAGTTTGGTTGGTGCGTCAGCAAAGTCATAGGTTTTGCCGGTTTGTGGGTCg                                                                                           <  2:15651/251‑1 (MQ=255)                                                                                                                                                  tGGCGACCAAGCCAGGCAGTGACCAGCGCCGCCGGAACACCGGTTTCTTCCAGTTCACCATCTTCACCCATTCCCGGAGCAAGGATGCTGATTTTAATCGGGTCGAGCATACTCCAGTTATCCGGAATATCTTTGAAGCCGTGCCAGCTTTCGCCCGGATGCATTACCCAGCAGTCCTGAACGGTGGTCAGCAGTTTGGTTGGTGCGTCAGCAAAGTCATAGGTTTTGCCGGTTTGTGGGTCGGTgacgac                                                                                   >  2:124137/1‑251 (MQ=255)                                                                                                                                                                            cgccgcCGGAACACCGGTTTCTTCCAGTTCACCATCTTCACCCATTCCCGGAGCAAGGATGCTGATTTTAATCGGGTCGAGCATACTCCAGTTATCCGGAATATCTTTGAAGCCGTGCCAGCTTTCGCCCGGATGCATTACCCAGCAGTCCTGAACGGTGGTCAGCAGTTTGGTTGGTGCGTCAGCAAAGTCATAGGTTTTGCCGGTTTGTGGGTCGGTGACGACTTCTTTGTTCCACGGTTTGAAGAAcc                                                         <  1:21345/251‑1 (MQ=255)                                                                                                                                                                                  cgGAACACCGGTTTCTTCCAGTTCACCATCTTCACCCATTCCCGGAGCAAGGATGCTGATTTTAATCGGGTCGAGCATACTCCAGTTATCCGGAATATCTTTGAAGCCGTGCCAGCTTTCGCCCGGATGCATTACCCAGCAGTCCTGAACGGTGGTCAGCAGTTTGGTTGGTGCGTCAGCAAAGTCATAGGTTTTGCCGGTTTGTGGGTCGGTGACGACTTCTTTGTTCCACGGTTTGAAGAACCAGCTAc                                                   <  2:93867/251‑1 (MQ=255)                                                                                                                                                                                        acCGGTTTCTTCCAGTTCACCATCTTCACCCATTCCCGGAGCAAGGATGCTGATTTTAATCGGGTCGAGCATACTCCAGTTATCCGGAATATCTTTGAAGCCGTGCCAGCTTTCGCCCGGATGCATTACCCAGCAGTCCTGAACGGTGGTCAGCAGTTTGGTTGGTGCGTCAGCAAAGTCATAGGTTTTGCCGGTTTGTGGGTCGGTGACGACTTCTTTGTTCCACGGTTTGAAGAACCAGCTACCGTCAg                                             <  1:145700/251‑1 (MQ=255)                                                                                                                                                                                               tcttcCAGTTCACCATCTTCACCCATTCCCGGAGCAAGGATGCTGATTTTAATCGGGTCGAGCATACTCCAGTTATCCGGAATATCTTTGAAGCCGTGCCAGCTTTCGCCCGGATGCATTACCCAGCAGTCCTTAACGGTGGTCAGCAGTTTGGTTGGTGCGTCAGCAAAGTCATAGGTTTTGCCGGTTTGTGGGTCGGTGACGACTTCTTTGTTCCACGGTTTGAAGAACCAGCTACCGTCAGCGGTGa                                       <  1:44451/250‑1 (MQ=255)                                                                                                                                                                                                  tcCAGTTCACCATCTTCACCCATTCCCGGAGCAAGGATGCTGATTTTAATCGGGTCGAGCATACTCCAGTTATCCGGAATATCTTTGAAGCCGTGCCAGCTTTCGCCCGGATGCATTACCCAGCAGTCCTGAACGGTGGTCAGCAGTTTGGTTGGTGCGTCAGCAAAGTCATAGGTTTTGCCGGTTTGTGGGTCGGTGACGACTTCTTTGTTCCACGGTTTGAAGAACCAGCTACCGTCAGCGGTGaa                                      <  1:63870/248‑1 (MQ=255)                                                                                                                                                                                                    cAGTTCACCATCTTCACCCATTCCCGGAGCAAGGATGCTGATTTTAATCGGGTCGAGCATACTCCAGTTATCCGGAATATCTTTGAAGCCGTGCCAGCTTTCGCCCGGATGCATTACCCAGCAGTCCTGAACGGTGGTCAGCAGTTTGGTTGGTGCGTCAGCAAAGTCATAGGTTTTGCCGGTTTGTGGGTCGGTGACGACTTCTTTGTTCCACGGTTTGAAGAACCAGCTACCGTCAGCGGTGAACTCtt                                 <  1:267564/251‑1 (MQ=255)                                                                                                                                                                                                     aGTTCACCATCTTCACCCATTCCCGGAGCAAGGATGCTGATTTTAATCGGGTCGAGCATACTCCAGTTATCCGGACTATCTTTGAAGCCGTGCCAGCTTTCGCCCGGATGCATTACCCAGCAGTCCTGAACGGTGGTCAGCAGTTTGGTTGGTGCGTCAGCAAAGTCATAGGTTTTGCCGGTTTGTGGGTCGGTGACGACTTCTTTGTTCCACGGTTTGAAGAACCAGCTACCGTCAGCGGTGAACTCttt                                >  2:226396/1‑251 (MQ=255)                                                                                                                                                                                                            cATCTTCACCCATTCCCGGAGCAAGGATGCTGATTTTAATCGGGTCGAGCATACTCCAGTTATCCGGAATATCTTTGAAGCCGTGCCAGCTTTCGCCCGGATGCATTACCCAGCAGTCCTGAACGGTGGGCAGCAGTTTGGTTGGTGCGTCAGCAAAGTCATAGGTTTTGCCGGTTTGTGGGTCGGTGACGACTTCTTTGTTCCACGGTTTGAAGAACCAGCTACCGGCAGCGGTGAACTCTTTATATAGc                         <  2:84854/251‑1 (MQ=255)                                                                                                                                                                                                                 tCACCCATTCCCGGAGCAAGGATGCTGATTTTAATCGGGTCGAGCATACTCCAGTTATCCGGAATATCTTTGAAGCCGTGCCAGCTTTCGCCCGGATGCATTACCCAGCAGTCCTGAACGGTGGTCAGCAGTTTGGTTGGTGCGTCAGCAAAGTCATAGGTTTTGCCGGTTTGTGGGTCGGTGACGACTTCTTTGTTCCACGGTTTGAAGAACCAGCTACCGTCAGCGGTGAACTCTTTATATAGCCgcgc                    >  1:73773/1‑251 (MQ=255)                                                                                                                                                                                                                          cccGGAGCAAGGATGCTGATTTTAATCGGGTCGAGCATACTCCAGTTATCCGGAATATCTTTGAAGCCGTGCCAGCTTTCGCCCGGATGCATTACCCAGCAGTCCTGAACGGTGGTCAGCAGTTTGGTTGGTGCGTCAGCAAAGTCATAGGTTTTGCCGGTTTGTGGGTCGGTGACGACTTCTTTGTTCCACGGTTTGAAGAACCAGCTACCGTCAGCGGTGAACTCTTTATATAGCCGCGCCATCGCCTg           >  2:356714/1‑251 (MQ=255)                                                                                                                                                                                                                               aGCAAGGATGCTGATTTTAATCGGGTCGAGCATACTCCAGTTATCCGGAATATCTTTGAAGCCGTGCCAGCTTTCGCCCGGATGCATTACCCAGCAGTCCTGAACGGTGGTCAGCAGTTTGGTTGGTGCGTCAGCAAAGTCATAGGTTTTGCCGGTTTGTGGGTCGGTGACGACTTCTTTGTTCCACGGTTTGAAGAACCAGCTACCGTCAGCGGTGAACTCTTTATATAGCCGCGCCATCGCCTGACgg       >  2:287124/1‑250 (MQ=255)                                                                                                                                                                                                                                 cAAGGATGCTGATTTTAATCGGGTCGAGCATACTCCAGTTATCCGGAATATCTTTGAAGCCGTGCCAGCTTTCGCCCGGATGCATTACCCAGCAGTCCTGAACGGTGGTCAGCAGTTTGGTTGGTGCGTCAGCAAAGTCATAGGTTTTGCCGGTTTGTGGGTCGGTGACGACTTCTTTGTTCCACGGTTTGAAGAACCAGCTACCGTCAGCGGTGAACTCTTTATATAGCCGCGCCATCGCCTGACGGaaa    <  1:207632/251‑1 (MQ=255)                                                                                                                                                                                                                                    ggATGCTGATTTTAATCGGGTCGAGCATACTCCAGTTATCCGGAATATCTTTGAAGCCGTGCCAGCTTTCGCCCGGATGCATTACCCAGCAGTCCTGAACGGTGGTCAGCAGTTTGGTTGGTGCGTCAGCAAAGTCATAGGTTTTGCCGGTTTGTGGGTCGGTGACGACTTCTTTGTTCCACGGTTTGAAGAACCAGCTACCGTCAGCGGTGAACTCTTTATATAGCCGCGCCATCGCCTGACGGAAATc  >  2:91768/1‑250 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                 CACCTGCGCCAGCGGTGTGTTGGCGTCATAGTGGCGTTTGAAGGAGCAAAGGGTGTTAACCAGAGTTCCCCATTTCCCACGGGTTACGCCCATAGAGAACAGGAACATAATTTGGAAGTCAGTGGTGCGGGTAGGTACAATGCCGTGGCGACCAAGCCAGGCAGTGACCAGCGCCGCCGGAACACCGGTTTCTTCCAGTTCACCATCTTCACCCATTCCCGGAGCAAGGATGCTGACTTTAATCGGGTCGAGCATACTCCAGTTATCCGGAATATCTTTGAAGCCGTGCCAGCTTTCGCCCGGATGCATTACCCAGCAGTCCTGAACGGTGGTCAGCAGTTTGGTTGGTGCGTCAGCAAAGTCATAGGTTTTGCCGGTTTGTGGGTCGGTGACGACTTCTTTGTTCCACGGTTTGAAGAACCAGCTACCGTCAGCGGTGAACTCTTTATATAGCCGCGCCATCGCCTGACGGAAATC  >  NC_000913/4338680‑4339156

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