Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	585,109	T→A	6.9%	Q49H (CAA→CAT)	ansA ←	type 1 L‑asparaginase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	585,109	0	T	A	6.9%	174.0 / 2.8	58	Q49H (CAA→CAT)	ansA	type 1 L‑asparaginase
Reads supporting (aligned to +/- strand):  ref base T (23/31);  new base A (2/2);  total (25/33)
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.60e-01

GGCCATGCTGTCGGTACCGTGCAGCACCAACACGCCATCGTGGCCCGCGACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGTGCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATGGTGCCACCGGTATAGAGAACGAGAAGATTCTTTACTGCGCGCATG  >  NC_002947/584962‑585256                                                                                                                                                    |                                                                                                                                                    ggCCATGCTGTCGGTACCGTGCAGCACCAACACGCCATCGTGGCCCGCGACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGTGCAGTGTCCAtt                                                                                                                                                    <  4:208087/149‑1 (MQ=255)    cATGCTGTCGGTACCGTGCAGCACCAACACGCCATCGTGGCCCGCGACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGTGCAGTGTCCATTGCa                                                                                                                                                 >  6:165091/1‑149 (MQ=255)          gTCGGTACCGTGCAGCACCAACACGCCATCGTGGCCCGCGACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGTGCAGTGTCCATTGCACCTggg                                                                                                                                           >  6:309887/1‑149 (MQ=255)           tCGGTACCGTGCAGCACCAACACGCCATCGTGGCCCGCGACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGTGCAGTGTCCATTGCACCTggg                                                                                                                                           <  3:260817/148‑1 (MQ=255)           tCGGTACCGTGCAGCACCAACACGCCATCGTGGCCCGCGACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGTGCAGTGTCCATTGCACCTggg                                                                                                                                           <  3:255856/148‑1 (MQ=255)                 ccGTGCAGCACCAACACGCCATCGTGGCCCGCGACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGTGCAGTGTCCATTGCACCCCAGGTGCaat                                                                                                                                    <  4:73726/149‑3 (MQ=255)                  cgggcagcaCCAACACGCCATCGTGGCCCGCGACATCCACCCCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGTGCAGTGTCCATTGCACCTGGGGTGCATTg                                                                                                                                   <  8:208850/146‑1 (MQ=255)                  cGTGCAGCACCAACACGCCATCGTGGCCCGCGACATCCACCGCCGCGACGATGGCATCGCGCATCGCTAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGTGCAGTGTCCATTGCACCTGGGGTGCATTg                                                                                                                                   <  3:217488/149‑1 (MQ=255)                    tGCAGCACCAACACGCCATCGTGGCCCGCGACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGTGCAGTGTCCATTGCACCTGGGGTGCATTGGc                                                                                                                                 <  4:294476/149‑1 (MQ=255)                      cagcaCCAACACGCCATCGTGGCCCGCGACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGTGCAGTGTCCAATGCACCCCAGGTGCAATGGACa                                                                                                                               <  3:3158/149‑1 (MQ=255)                           ccAACACGCCATCGTGGCCCGCGACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGTGCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTgg                                                                                                                          <  7:163558/149‑1 (MQ=255)                               cacGCCATCGTGGCCCGCGACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGTGCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCa                                                                                                                       >  2:314907/1‑148 (MQ=255)                                acGCCATCGTGGCCCGCGACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGTGCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCaaa                                                                                                                     <  6:165562/149‑1 (MQ=255)                                acGCCATCGTGGCCCGCGACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGTGCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCaaa                                                                                                                     <  2:246829/149‑1 (MQ=255)                                acGCCATCGTGGCCCGCGACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGTGCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCaaa                                                                                                                     <  2:140847/149‑1 (MQ=255)                                acGCCATCGTGGCCCGCGACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGTGCAGTGGCCATTGCACCTGGGGTGCATTGGCCATCTGGGCaaa                                                                                                                     <  5:97114/149‑1 (MQ=255)                                  gCCATCGTGGCCCGCGACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGTGCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAgt                                                                                                                   <  5:309887/149‑1 (MQ=255)                                        gTGGCCCGCGACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGTGCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCa                                                                                                              <  1:329973/148‑1 (MQ=255)                                        gTGGCCCGCGACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGTGCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCAc                                                                                                             <  6:93907/149‑1 (MQ=255)                                            cccGCGACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGTGCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCAt                                                                                                         >  8:253063/1‑149 (MQ=255)                                             ccGCGACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGTGCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATc                                                                                                        >  5:368667/1‑149 (MQ=255)                                             ccGCGACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGTGCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATc                                                                                                        >  2:163627/1‑149 (MQ=255)                                               gcgACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGTGCAGTGTCCAATGCACCCCAGGTGCAATgga                                                                                                                                 <  4:381936/122‑2 (MQ=255)                                               gcgACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGTGCAGTGTCCAATGCACCCCAGGTGCAATgga                                                                                                                                 >  3:381936/1‑121 (MQ=255)                                                  aCATCCCCCCCCGCGACGCTGGCATCGCGCATCGCCAGCCAGCTGTGCGGGTGCATATTGGCGCTGTCGAGCAGCGGGGTCAGCTCGTGCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGc                                                                                                   <  8:316969/149‑1 (MQ=255)                                                       cACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGTGCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCaa                                                                                              <  4:59077/149‑1 (MQ=255)                                                         ccgccgCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGTGCATTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCaaaa                                                                                            <  1:311701/149‑1 (MQ=255)                                                         ccgccgCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGTGCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCaaaa                                                                                            <  2:381682/149‑1 (MQ=255)                                                                gacgaTGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGTGCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCAccgg                                                                                     >  2:40167/1‑149 (MQ=255)                                                                 acgaTGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGTGCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCAccggc                                                                                    >  5:4429/1‑149 (MQ=255)                                                                   gaTGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGTGCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCAccggccg                                                                                  >  8:318889/1‑149 (MQ=255)                                                                       gcatcgcgcatcgcCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGTGCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAATCACCGGCCGGTGc                                                                              <  7:317382/149‑1 (MQ=255)                                                                        catcgcgcatcgcCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGTGCAGTGTCCATTGCACCTGGGGTGCATTg                                                                                                                                   >  2:271215/1‑95 (MQ=255)                                                                        catcgcgcatcgcCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGTGCAGTGTCCATTGCACCTGGGGTGCATTg                                                                                                                                   <  1:271215/95‑1 (MQ=255)                                                                            gcgcatcgcCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGTGCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAgg                                                                          <  3:15759/148‑1 (MQ=255)                                                                                   gccagccaGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGTGCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCg                                                                   <  3:230362/148‑1 (MQ=255)                                                                                    ccagccagTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGTGCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCggg                                                                 >  7:128769/1‑149 (MQ=255)                                                                                        ccagTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGTGCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGttt                                                             >  7:237474/1‑149 (MQ=255)                                                                                          agTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGTGCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAACCCACCGGCCGGTGCCAGGCCTTCCGGGGTTTc                                                            >  8:250877/1‑148 (MQ=255)                                                                                            ttGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGTGCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAg                                                         <  7:220081/149‑1 (MQ=255)                                                                                                    ttGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGTGCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGa                                                 >  7:46389/1‑149 (MQ=255)                                                                                                      gCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGTGCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATg                                               >  6:242212/1‑149 (MQ=255)                                                                                                      gCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGTGCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATg                                               >  5:279860/1‑149 (MQ=255)                                                                                                       cATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGTGCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATgg                                              >  8:205856/1‑149 (MQ=255)                                                                                                              gcgcTGTCGAGCAGCGGGTTCAGCTCGTGCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATGGTGCCAcc                                       <  1:40167/149‑1 (MQ=255)                                                                                                                gcTGTCGAGCAGCGGGTTCAGCTCGTGCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATGGTGCCACCgg                                     <  7:253063/149‑1 (MQ=255)                                                                                                                gcTGTCGAGCAGCGGGTTCAGCTCGTGCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATGGTGCCACCgg                                     >  7:42517/1‑149 (MQ=255)                                                                                                                   gTCGAGCAGCGGGTTCAGCTCGTGCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATGGTGCCACCGGta                                   >  8:25752/1‑148 (MQ=255)                                                                                                                    tCGAGCAGCGGGTTCAGCTCGTGCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATGGTGCCACCgg                                     <  6:115817/145‑1 (MQ=255)                                                                                                                    tCGAGCAGCGGGTTCAGCTCGTGCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATGGTGCCACCgg                                     >  5:115817/1‑145 (MQ=255)                                                                                                                    tCGAGCAGCGGGTTCAGCTCGTGCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATGGTGCCACCGGtata                                 <  6:279860/149‑1 (MQ=255)                                                                                                                    tCGAGCAGCGGGTTCAGCTCGTGCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATGGTGCCACCGGtata                                 <  5:165091/149‑1 (MQ=255)                                                                                                                     cGAGCAGCGGGTTCAGCTCGTGCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCTGATGGTGCCACCGGTATag                                <  8:353308/149‑1 (MQ=255)                                                                                                                          agcGGGTTCAGCTCGTGCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTGGAGCATGCCGATGGTGCCACCGGTATAGAGAAc                           >  6:227136/1‑149 (MQ=255)                                                                                                                            cGGGTTCAGCTCGTGCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTAACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATGGTGCCACCGGTATAGAGAAc                           >  3:52019/1‑147 (MQ=255)                                                                                                                                          gCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAGAACCTCCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATGGTGCCACCGGTATAGAGAACGAGAAGATTCTTTACt           <  7:250877/149‑1 (MQ=255)                                                                                                                                           cAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATGGTGCCACCGGTATAGAGAACGAGAAGATTCTTTACTg          >  4:88791/1‑149 (MQ=255)                                                                                                                                                tCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATGGTGCCACCGGTATAGAGAACGAGAAGATTCTTTACTgcgcg      <  7:318889/148‑1 (MQ=255)                                                                                                                                                 ccAATGCACCCCAGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACAACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATGGTGACACCGGTATAGAGCACGAGAAGATTCTTTACTGCGCGCa    >  2:296996/1‑149 (MQ=255)                                                                                                                                                    ttGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATGGTGCCACCGGTATAGAGAACGAGAAGATTCTTTACt           <  6:316237/139‑1 (MQ=255)                                                                                                                                                    ttGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATGGTGCCACCGGTATAGAGAACGAGAAGATTCTTTACt           >  5:316237/1‑139 (MQ=255)                                                                                                                                                    ttGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATGGTGCCACCGGTATAGAGAACGAGAAGATTCTTTACTGCGCGCATg  >  6:33228/1‑148 (MQ=255)                                                                                                                                                    ttGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATGGTGCCACCGGTATAGAGAACGAGAAGATTCTTTACTGCGCGCATg  >  2:293692/1‑148 (MQ=255)                                                                                                                                                    |                                                                                                                                                    GGCCATGCTGTCGGTACCGTGCAGCACCAACACGCCATCGTGGCCCGCGACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGTGCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATGGTGCCACCGGTATAGAGAACGAGAAGATTCTTTACTGCGCGCATG  >  NC_002947/584962‑585256

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