Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,545,220	G→A	5.8%	P113P (CCG→CCA)	gntX →	DNA catabolic protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,545,220	0	G	A	5.8%	218.6 / 2.1	69	P113P (CCG→CCA)	gntX	DNA catabolic protein
Reads supporting (aligned to +/- strand):  ref base G (34/31);  new base A (2/2);  total (36/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 = 1.00e+00

TCACGGTTGCCGACTATGCGCCGCCGTTAAGTCCGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCAGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCCGGATCGCATCGTCAGCGTTCCGTTATGGCAGCGGCGTCACTGGCGTCGGGGATTTAATCAGAGCGATTTGCTGTGTCAGCCGTTATCACGCTGGTTGCACTGCCAATGGGATAGCGAAGCCGTCACACGTACACGGGCCACTGCGAC  >  NC_000913/3545075‑3545366                                                                                                                                                  |                                                                                                                                                   tCACGGTTGCCGACTATGCGCCGCCGTTAAGTCCGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCAGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCCGGAt                                                                                                                                                 <  3:157459/149‑1 (MQ=255)  cACGGTTGCCGACTATGCGCCGCCGTTAAGTCCGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCAGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCCGGATc                                                                                                                                                <  7:212781/149‑1 (MQ=255)   aCGGTTGCCGACTATGCGCCGCCGTTAAGTCCGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCAGCGCCCTGTCCCGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCCGGATCg                                                                                                                                               >  1:164011/1‑149 (MQ=255)   aCGGTTGCCGACTATGCGCCGCCGTTAAGTCCGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCAGCGCCCTGTCATGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCCGGATc                                                                                                                                                <  3:107816/148‑1 (MQ=255)     ggTTGCCGACTATGCGCCGCCGTTAAGTCCGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCAGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCCGGATCGCa                                                                                                                                             >  1:86763/1‑149 (MQ=255)          ccGACTATGCGCCGCCGTTAAGTCCGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCAGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCCGGATCGCATCGt                                                                                                                                         >  8:29763/1‑148 (MQ=255)           cGACTATGCGCCGCCGTTAAGTCCGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCAGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCCGGATCGCATCGTCa                                                                                                                                       <  8:38615/149‑1 (MQ=255)            gACTATGCGCCGCCGTTAAGTCCGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCAGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCCGGATCGCATCGTCAg                                                                                                                                      >  2:229737/1‑149 (MQ=255)                   cgccgcCGTTAAGTCCGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCAGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCCGGATCGCATCGTCAGCGTTcc                                                                                                                                <  1:45976/148‑1 (MQ=255)                   cgccgcCGTTAAGTCCGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCAGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCCGGATCGCATCGTCAGCGTTCCg                                                                                                                               <  7:29763/149‑1 (MQ=255)                     ccgccgTTAAGTCCGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCAGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCCGGATCGCATCGTCAGCGTTCCGt                                                                                                                              >  8:58725/1‑148 (MQ=255)                       gccgTTAAGTCCGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCAGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCTGGATCGCATCGTCAGCGTTCCGTTAt                                                                                                                           <  5:156463/149‑1 (MQ=255)                        ccgTTAAGTCCGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCAGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCCGGATCGCATCGTCAGCGTTCCGTTATg                                                                                                                          <  1:229737/149‑1 (MQ=255)                        ccgTTAAGTCCGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCAGCGCCCTGTCACCTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCCGGATCGCATCGTCAGCGTTCCGTTATg                                                                                                                          >  4:182481/1‑149 (MQ=255)                         cgTTAAGTCCGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCAGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCCGGATCGCATCGTCAGCGTTCCGTTATgg                                                                                                                         >  4:49238/1‑149 (MQ=255)                           ttAAGTCCGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCAGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCCGGATCGCATCGTCAGCGTTCCGTTATgg                                                                                                                         >  8:141741/1‑147 (MQ=255)                             aaGTCCGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCAGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCCGGATCGCATCGTCAGCGTTCCGTTATGGCAgc                                                                                                                     >  8:137116/1‑149 (MQ=255)                              aGTCCGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCAGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCCGGATCGCATCGTCAGCGTTCCGTTATGGCAgc                                                                                                                     <  2:86763/148‑1 (MQ=255)                              aGTCCGCTTATCCACCACCTTAAATTTTCCCGGCGCAGTGAAATCGCCAGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCCGGATCGCATCGTCAGCGTTCCGTTATGGCAgcg                                                                                                                    <  8:234985/149‑1 (MQ=255)                               gTCCGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCAGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCCGGATCGCATCGTCa                                                                                                                                       >  6:27371/1‑129 (MQ=255)                               gTCCGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCAGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCCGGATCGCATCGTCAGCGTTCCGTTATGGCAgcgg                                                                                                                   >  1:219718/1‑149 (MQ=255)                                  cGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCAGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCCGGATCGCATCGTCAGCGTTCCGTTATGGCAGCGGCGt                                                                                                                <  4:116883/149‑1 (MQ=255)                                   gCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCAGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCCAGATCGCATCGTCAGCGTTCCGTTATGGCAGCGGCGt                                                                                                                >  6:193314/1‑148 (MQ=255)                                       aTCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCAGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCCGGATCGCATCGTCAGCGTTCCGTTATGGCAGCGGCGTCAc                                                                                                             <  5:178267/147‑1 (MQ=255)                                          caccaGCTTAAATTTTCCCGGCGCAGTGAAATCGCCAGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCCGGATCGCATCGTCAGCGTTCCGt                                                                                                                              <  1:35104/127‑1 (MQ=255)                                          caccaGCTTAAATTTTCCCGGCGCAGTGAAATCGCCAGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCCGGATCGCATCGTCAGCGTTCCGt                                                                                                                              >  2:35104/1‑127 (MQ=255)                                           accaGCTTAAATTTTCCCGGCGCAGTGAAATCGCCAGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCCGGATCGCATCGTCAGCGTTCCGTTATGGCAGCGGCGTCACTGGCGt                                                                                                       >  1:20768/1‑149 (MQ=255)                                           accaGCTTAAATTTTCCCGGCGCAGTGAAATCGCCAGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCCGGATCGCATCGTCAGCGTTCCGTTATGGCAGCGGCGTCACTGGCGt                                                                                                       >  2:50485/1‑149 (MQ=255)                                             ccgCTTAAATTTTCCCGGCGCAGTGAAATCGCCCGCGCCCTGTCCCGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTGGCCATTGCCGGATCGCATCGTCCGCGTTCCGTTATGGCAGCGGCGTCACTGGCGTCg                                                                                                     >  7:174393/3‑149 (MQ=255)                                                cTTAAATTTTCCCGGCGCAGTGAAATCGCCAGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCCGGATCGCATCGTCGGCGTTCCGTTATGGCAGCGGCGTCACTGGCGTCggg                                                                                                   >  2:222114/1‑148 (MQ=255)                                                cTTAAATTTTCCCGGCGCAGTGAAATCGCCAGCGCCCTGTCACGGCTGTTACTATTGGAAGTCTTACCCGCTCGTCGCACCACCGGGTTGCAATTGCCAGATCGCATCGTCAGCGTTCCGTTATGGCAGCGGCGTCACTGGCGTCggg                                                                                                   >  1:110946/1‑148 (MQ=255)                                                   aaaTTTTCCCGGCGCAGTGAAATCGCCAGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCCGGATCGCATCGTCAGCGTTCCGTTATGGCAGCGGCGTCACTGGCGTCGGGGa                                                                                                 <  3:49238/147‑1 (MQ=255)                                                    aaTTTTCCCGGCGCAGTGAAATCGCCAGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCCGGATCGCATCGTCAGCGTTCCGTTATGGCAGCGGCGTCACTGGCGTCGGGGAttt                                                                                              <  1:214020/149‑1 (MQ=255)                                                     aTTTTCCCGGCGCAGTGAAATCGCCAGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCCGGATCGCATCGTCAGCGTTCCGTTATGGCAGCGGCGTCACTGGCGTCGGGGATTTa                                                                                             <  1:222114/149‑1 (MQ=255)                                                        ttCCCGGCGCAGTGAAATCGCCAGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCGCGTCGCACCACCGGGGTGCAATTGCCGGATCGCATCGTCAGCGTTCCGTTATGGCAGCGGCGTCACTGGCGTCGGGGATTTAAt                                                                                           <  2:120849/148‑1 (MQ=255)                                                          cccGGCGCAGTGAAATCGCCAGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCCGGATCGCATCGTCAGCGTTCCGTTATGGCAGCGGCGTCACTGGCGTCGGGGATTTAATCa                                                                                         <  2:8553/148‑1 (MQ=255)                                                           ccGGCGCAGTGAAATCGCCAGCGCCCTGTCAAGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCCGGAGCGCATCGTCAGCGTTCCGTTATGGCAGCGGCGTCACTGGCGTCGGGGATTTAATCa                                                                                         >  8:43297/1‑147 (MQ=255)                                                                   gTGAAATCGCCAGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCCGGATCGCATCGTCAGCGTTCCGTTATGGCAGCGGCGTCCCTGGCGTCGGGGATTTAATCAGAGCGATTTg                                                                               >  5:128502/1‑149 (MQ=255)                                                                     gAAATCGCCAGCGCCCTGTCACGTCTGTTACTATTGGATGTCTTAGACTCTCGTCGCACCACCGGGTTGCAATTGCCGGATCGCATCGTCTGGGTTCCGTTATGGCAGCGGCGTCAATTGCGTCGTGGATTTAATCAGAGCTATTTGCt                                                                             <  8:109809/149‑1 (MQ=255)                                                                       aaTCGCCAGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCCGGATCGCATCGTCAGCGTTCCGTTATGGCAGCGGCGTCACTGGCGTCGGGGATTTAATCAGAGCGATTTGCtgt                                                                           >  6:141767/1‑149 (MQ=255)                                                                       aaTCGCCAGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCCGGATCGCATCGTCAGCGTTCCGTTATGGCAGCGGCGTCACTGGCGTCGGGGATTTAATCAGAGCGATTTGCtgt                                                                           <  2:20768/149‑1 (MQ=255)                                                                         tCGCCAGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCCGGATCGCATCGTCAGCGTTCCGTTATGGCAGCGGCGTCACTGGCGTCGGGGATTTAATCAGAGCGATTTGCTgtgt                                                                         <  7:137116/149‑1 (MQ=255)                                                                            ccAGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCCGGATCGCATCGTCAGCGTTCCGTTATGGCAGCGGCGTCACTGGCGTCGGGGATTTAATCAGAGCGATTTGCTGTGTCAg                                                                      <  7:141741/149‑1 (MQ=255)                                                                                cgcCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCCGGATCGCATCGTCAGCGTTCCGTTATGGCAGCGGCGTCACTGGCGTCGGGGATTTAATCAGAGCGATTTGCTGTGTCAGCCGt                                                                  <  1:50485/149‑1 (MQ=255)                                                                                cgcCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGGCGCACCACCGGGTTGCAATTGCCAGATCGCATCGTCAGCGTTCCGTTATGGCAGCGGCGTCACTGGCGTCGGGGATTTAATCAGAGCGATTTGCTGTGTCAGCCGt                                                                  <  2:110946/149‑1 (MQ=255)                                                                                         aCGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCCGGATCGCATCGTCAGCGTTCCGTTATGGCAGCGGCGTCACTGGCGTCGGGGATTTAATCAGAGCGATTTGCTGTGTCAGCCGTTATCACGCt                                                         >  3:225755/1‑149 (MQ=255)                                                                                           gTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCCGGATCGCATCGTCAGCGTTCCGTTATGGCAGCGGCGTCACTGGCGTCGGGGATTTAATCAGAGCGATTTGCTGTGTCAGCCGt                                                                  <  3:174193/138‑1 (MQ=255)                                                                                           gTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCCGGATCGCATCGTCAGCGTTCCGTTATGGCAGCGGCGTCACTGGCGTCGGGGATTTAATCAGAGCGATTTGCTGTGTCAGCCGt                                                                  >  4:174193/1‑138 (MQ=255)                                                                                           gTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCCGGATCGCATCGTCAGCGTTCCGTTATGGCAGCGGCGTCACTGGCGTCGGGGATTTAATCAGAGCGATTTGCTGTGTCAGCCGTTATCACGCTgg                                                       <  4:225755/149‑1 (MQ=255)                                                                                               gTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCCGGATCGCATCGTCAGCGTTCCGTTATGGCAGCGGCGTCCCTGGCGTCGGGGATTTAATCCGAGCGATTTGCTGTGTCAGCCGTGATCACGCTGGTTg                                                    >  6:188376/1‑148 (MQ=255)                                                                                                ttACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCCGGATCGCATCGTCAGCGTTCCGTTATGGCAGCGGCGTCACTGGCGTCGGGGATTTAATCAGAGCGATTTGCTGTGTCAGCCGTTATCACGCTGGTTGc                                                   <  5:141767/148‑1 (MQ=255)                                                                                                  aCTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCCGGATCGCATCGTCAGCGTTCCGTTATGGCAGCGGCGTCACTGGCGTCGGGGATTTAATCAGAGCGATTTGCTGTGTCAGCCGTTATCACGCTGGTTGCACt                                                >  7:203694/1‑149 (MQ=255)                                                                                                    tATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCCGGATCGCATCGTCAGCGTTCCGTTATGGCAGCGGCGTCACTGGCGTCGGGGATTTAATCAGAGCGATTTGCTGTGTCAGCCGTTATCACGCTGGTTGCACTGc                                              >  7:30492/1‑149 (MQ=255)                                                                                                    tATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCCGGATCGCATCGTCAGCGTTCCGTTATGGCAGCGGCGTCACTGGCGTCGGGGATTTAAGCAGAGCTATTTGCTTTGTCAGCCGTTTTCACGATGGTTGCACTGc                                              >  8:171437/1‑149 (MQ=255)                                                                                                        ggAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCCAGATCGCATCGTCAGCGTTCCGTTATGGCAGCGGCGTCACTGGCGTCGGGGATTTAATCAGAGCGATTTGCTGTGTCAGCCGTTATCACGCTGGTTGCACTGCCaa                                           <  5:193314/148‑1 (MQ=255)                                                                                                            gTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCCGGATCGCATCGTCAGCGTTCCGTTATGGCAGCGGCGTAACTGGCGTCGGGGATTTAATCAGAGCGATTTGCTGTGTCAGCCGTTATCACGCTGGTTGCACTGCCAATGGGa                                      >  1:216290/1‑149 (MQ=255)                                                                                                                 acacGCTCGTCGCACCACCGGGTTGCAATTGCCGGATCGCATCGTCAGCGTTCCGTTATGGCAGCGGCGTCACTGGCGTCGGGGATTTAATCAGAGCGATTTGCTGTGTCAGCCGTTATCACGCTGGTTGCACTGCCAATGGGATa                                    <  7:41172/146‑1 (MQ=255)                                                                                                                 acacGCTCGTCGCACCACCGGGTTGCAATTGCCGGATCGCATCGTCAGCGTTCCGTTATGGCAGCGGCGTCACTGGCGTCGGGGATTTAATCAGAGCGATTTGCTGTGTCAGCCGTTATCACGCTGGTTGCACTGCCAATGGGATa                                    <  8:203694/146‑1 (MQ=255)                                                                                                                   acGCTCGTCGCACCACCGGGTTGCAATTGCCGGATCGCATCGTCAGCGTTCCGTTATGGCAGCTGCGTCACTGGCGTCGGGGATTTAATCAGAGCGATTTGATGTTTCAGCCGTTATCACGCTGGTTGCACTGCCAATGGGATAGCg                                 >  4:11873/1‑147 (MQ=255)                                                                                                                   acGATCGTCGCCCCACCGGGTTGCAATTGCCGGAGCGCAGCGTCAGCGTTCCGTTATGGCAGCGGCGTCACTGGCGTCGGGGATTTAATCAGAGCGATTTGCTGTGTCAGCCGTTATCACGCTGGTTGCACTGCCAATGGGATAGCGaa                               <  8:174393/149‑1 (MQ=255)                                                                                                                     gCTCGTCGCACCACCGGGTTGCAATTGCCGGATCGCATCGTCAGCGTTCCg                                                                                                                               >  8:167118/1‑51 (MQ=255)                                                                                                                     gCTCGTCGCACCACCGGGTTGCAATTGCCGGATCGCATCGTCAGCGTTCCg                                                                                                                               <  7:167118/51‑1 (MQ=255)                                                                                                                      cTCGTCGCACCACCGGGTTGCAATTGCCGGATCGCATCGTCAGCGTTCCGTTATGGCAGCGGCGTCACTGGCGTCGGGGATTTAATCAGAGCGATTTGCTGTGTCAGCCGTTATCACGCTGGTTGCACTGCCAATGGGATAGCGAAGcc                            <  2:219718/149‑1 (MQ=255)                                                                                                                              accaccGGGTTGCAATTGCCGGATCGCATCGTCAGCGTTCCGTTATGGCAGCGGCGTCACTGGCGTCGGGGATTTAATCAGAGCGATTTGCTGTGTCAGCCGTTATCACGCTGGTTGCACTGCCAATGGGATAGCGAAGCCGTCACACg                    >  2:68953/1‑149 (MQ=255)                                                                                                                              accaccGGGTTGCAATTGCCGGATCGCATCGTCAGCGTTCCGTTATGGCAGCGGCGTCACTGGCGTCGGGGATTTAATCAGAGCGATTTGCTGTGTCAGCCGTTATCACGCTGGTTGCAATGCCAATGGGATAGCGAAGCCGTCACACg                    >  5:158633/1‑149 (MQ=255)                                                                                                                                           aaTTGCCGGATCGCATCGTCAGCGTTCCGTTATGGCAGCGGCGTCACTGGCGTCGGGGATTTAATCAGAGCGATTTGCTGTGTCAGCCGTTATCACGCTGGTTGCACTGCCAATGGGATAGCGAAGCCGTCACACGTACACGGGCCACt       >  7:251664/1‑149 (MQ=255)                                                                                                                                           aaTTGCCGGATCGCATCGTCAGCGTTCCGTTATGGCAGCGGCGTCACTGGCGTCGGGGATTTAATCAGAGCGATTTGCTGTGTCAGCCGTTATCACGCTGGTTGCACTGCCAATGGGATAGCGAAGCCGTCACACGTACACGGGCCACt       >  4:228970/1‑149 (MQ=255)                                                                                                                                               gCCGGATCGCATCGTCAGCGTTCCGTTATGGCAGCGGCGTCACTGGCGTCGGGGATTTAATCAGAGCGATTTGCTGTGTCAGCCGTTATCACGCTGGTTGCACTGCCAATGGGATAGCGAAGCCGTCACACGTACACGGGCCACTGCGa   >  1:34812/1‑149 (MQ=255)                                                                                                                                                ccGGATCGCATCGTCAGCGTTCCGTTATGGCAGCGGCGTCACTGGCGTCGGGGATTTAATCAGAGCGATTTGCTGTGTCAGCCGTTATCACGCTGGTTGCACTGCCAATGGGATAGCGAAGCCGTCACACGTACACGGGCCACTGCGAc  >  7:229982/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                   TCACGGTTGCCGACTATGCGCCGCCGTTAAGTCCGCTTATCCACCAGCTTAAATTTTCCCGGCGCAGTGAAATCGCCAGCGCCCTGTCACGTCTGTTACTATTGGAAGTCTTACACGCTCGTCGCACCACCGGGTTGCAATTGCCGGATCGCATCGTCAGCGTTCCGTTATGGCAGCGGCGTCACTGGCGTCGGGGATTTAATCAGAGCGATTTGCTGTGTCAGCCGTTATCACGCTGGTTGCACTGCCAATGGGATAGCGAAGCCGTCACACGTACACGGGCCACTGCGAC  >  NC_000913/3545075‑3545366

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