Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,124,811	T→C	6.0%	Q312Q (CAA→CAG)	wcaI ←	putative glycosyl transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,124,811	0	T	C	6.0%	203.7 / 2.5	66	Q312Q (CAA→CAG)	wcaI	putative glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base T (33/29);  new base C (2/2);  total (35/31)
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.83e-01

GTTCAACGCAAACCGCAATGCCCGGAAAGGTTTCGCAAAGCTGCCCCAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGTTACCGCCTACTGCCAGAATATTGGTCAGTTTCGACGGCAATACGGCATCTGCCGCGCCGCGTTTTTGCACCACCAGATGGCAATCGCCCATCTTCAGCAGTGCGGGTAAAGCGTCATACGATTGCAGCGGGAAAAATTGCATGTTGCGCAGTCCACGCTGCTGCGCCATTTTTTCCAGCCGCGCTTTGCCGCCGCCCTGCCCGACAATG  >  NC_000913/2124666‑2124955                                                                                                                                                  |                                                                                                                                                 gTTCAACGCAAACCGCAATGCCCGGAAAGGTTTCGCAAAGCTGCCCCAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGTTACCGCCTACTGCCAGAATATTGGTCAGTTTCGCCGGCAATACGGCATCTGCCGCGCCGCGTTTTTGc                                                                                                                                               <  1:332984/149‑1 (MQ=255) gTTCAACGCAAACCGCAATGCCCGGAAAGGTTTCGCAAAGCTGCCCCAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGTTACCGCCTACTGCCAGAATATTGGTCAGTTTCGACGGCAATACGGCATCTGCCGCGCCGCGTTTTTGc                                                                                                                                               <  8:110362/149‑1 (MQ=255) gTTCAACGCAAACCGCAATGCCCGGAAAGGTTTCGCAAAGCTGCCCCAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGTTACCGCCTACTGCCAGAATATTGGTCAGTTTCGACGGCAATACGGCATCTGCCGCGCCGCGTTTTTGc                                                                                                                                               <  6:44689/149‑1 (MQ=255)  ttCAACGCAAACCGCAATGCCCGGAAAGGTTTCGCAAAGCTGCCCCAGTTCTGTGTAGGCTTCAGCAGTACTCACCGCGTTACCGCCTACTGCCAGAATATTGGTCAGTTTCGACGGCAATACGGCATCTGCCGCGCCGCGTTTTTGca                                                                                                                                              >  4:295370/1‑149 (MQ=255)     aaCGCAAACCGCAATGCCCGGAAAGGTTTCGCAAAGCTGCCCCAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGTTACCGCCTACTGCCAGAATATTGGTCAGTTTCGACGGCAATACGGCATCTGCCGCGCCGCGTTTCTGcacca                                                                                                                                           >  7:58978/1‑149 (MQ=255)     aaCGCAAACCGCAATGCCCGGAAAGGTTTCGCAAAGCTGCCCCAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGTTACCGCCTACTGCCAGAATATTGGTCAGTTTCGACGGCAATACGGCATCTGCCGCGCCGCGTTTCTGcacca                                                                                                                                           >  7:53003/1‑149 (MQ=255)      aCGCAAACCGCAATGCCCGGAAAGGTTTCGCAAAGCTGCCCCAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGTTACCGCCTACTGCCAGAATATTGGTCAGTTTCGACGGCAATACGGCATCTGCCGCGCCGCGTTTTTGcaccac                                                                                                                                          >  3:192463/1‑149 (MQ=255)      aCGCAAACCGCAATGCCCGGAAAGGTTTCGCAAAGCTGCCCCAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGTTACCGCCTACTGCCAGAATATTGGTCAGTTTCGACGGCAATACGGCATCTGCCGCGCCGCGTTTTTGcaccac                                                                                                                                          >  5:376521/1‑149 (MQ=255)            aCCGCAATGCCCGGAAAGGTTTCGCAAAGCTGCCCCAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGTTACCGCCTACTGCCAGAATATTGGTCAGTTTCGACGGCAATACGGCATCTGCCGCGCCGCGTTTTTGCACCACCAGATg                                                                                                                                    <  2:186873/149‑1 (MQ=255)            aCCGCAATGCCCGGAAAGGTTTCGCAAAGCTGCCCCAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGTTACCGCCTACTGCCAGAATATTGGTCAGTTTCGACGGCAATACGGCATCTGCCGCGCCGCGTTTTTGCACCACCAGATg                                                                                                                                    <  1:155374/149‑1 (MQ=255)            aCCGCAATGCCCGGAAAGGTTTCGCAAAGCTGCCCCAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGTTACCGCCTACTGCCAGAATATTGGTCAGTTTCGACGGCAATACGGCATCTGCCGCGCCGCGTTTTTGCACCACCAGATg                                                                                                                                    <  6:376521/149‑1 (MQ=255)             ccGCAATGCCCGGAAAGGTTTCGCAAAGCTGCCCCAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGTTACCGCCTACTGCCAGAATATTGGTCAGTTTCGACGGCAATACGGCATCTGCCTCGCCGCGTTTCTGCACCACCAGATgg                                                                                                                                   <  8:58978/149‑1 (MQ=255)             ccGCAATGCCCGGAAAGGTTTCGCAAAGCTGCCCCAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGTTACCGCCTACTGCCAGAATATTGGTCAGTTTCGACGGCAATACGGCATCTGCCGCGCCGCGTTTCTGCACCACCAGATgg                                                                                                                                   <  8:53003/149‑1 (MQ=255)              cGCAATGCCCGGAAAGGTTTCGCAAAGCTGCCCCAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGTTACCGCCTACTGCCAGAATATTGGTCAGTTTCGACGGCAATACGGCATCTGCCGCGCCGCGTTTTTGc                                                                                                                                               >  1:14557/1‑136 (MQ=255)              cGCAATGCCCGGAAAGGTTTCGCAAAGCTGCCCCAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGTTACCGCCTACTGCCAGAATATTGGTCAGTTTCGACGGCAATACGGCATCTGCCGCGCCGCGTTTTTGc                                                                                                                                               <  2:14557/136‑1 (MQ=255)                   tGCCCGGAAAGGTTTCGCAAAGCTGCCCCAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGTTACCGCCTACTGCCAGAATATTGGTCAGTTTCGCCGGCAATACGGCATCTGCCGCGCCGCGTTTTTGCACCACCAGATGGCAATCg                                                                                                                             >  2:244503/1‑149 (MQ=255)                     cccGGAAAGGTTTCGCAAAGCTGCCCCAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGTTACCGCCTACTGCCAGAATATTGGTCAGTTTCGACGGCAATACGGCATCTGCCGCGCCGCGTTTTTGCACCACCAGATGGCAATCGc                                                                                                                            <  4:246227/148‑1 (MQ=255)                      ccGGAAAGGTTTCGCAAAGCTGCCCCAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGTTACCGCCTACTGCCAGAATATTGGTCAGTTTCGACGGCAATACGGCATCTGCCGCGCCGCGTTTTTGCACCACCAg                                                                                                                                       <  1:145176/136‑1 (MQ=255)                      ccGGAAAGGTTTCGCAAAGCTGCCCCAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGTTACCGCCTACTGCCAGAATATTGGTCAGTTTCGACGGCAATACGGCATCTGCCGCGCCGCGTTTTTGCACCACCAg                                                                                                                                       >  2:145176/1‑136 (MQ=255)                         gAAAGGTTTCGCAAAGCTGCCCCAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGTTACCGCCTACTGCCAGAATATTGGTCAGTTTCGACGGCAATACGGCATCTGCCGCGCCGCGTTTTTGCACCACCAGATGGCAATCGCCCATc                                                                                                                       >  4:41488/1‑149 (MQ=255)                         gAAAGGTTTCGCAAAGCTGCCCCAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGTTACCGCCTACTGCCAGAATATTGGTCAGTTTCGACGGCAATACGGCATCTGCCGCGCCGCGTTTTTGCACCACCAGATGGCAATCGCCCATc                                                                                                                       >  6:16508/1‑149 (MQ=255)                            aGGTTTCGCAAAGCTGCCCCAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGGTACCGCCTACTGCCAGAATATTGGTCAGTTTCGACGGCAATACGGCATCTGCCGCGCCGCGTTTTTGCACCACCAGATGGCAATCGCCCATCTTc                                                                                                                    <  5:16508/149‑1 (MQ=255)                                 tCGCAAAGCTGCCCCAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGTTACCGCCTACTGCCAGAATATTGGTCAGTTTCGACGGCAATACGGCATCTGCCGCGCCGCGTTTTTGCACCACCAGATGGCAATCGCCCATCTTcagca                                                                                                                <  1:244503/148‑1 (MQ=255)                                   gCAAAGCTGCCCCAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGTTACCGCCTACTGCCAGAATATTGGTCAGTTTCGACGGCAATACGGCATCTGCCGCGCCGCGTTTTTg                                                                                                                                                <  1:358326/114‑1 (MQ=255)                                   gCAAAGCTGCCCCAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGTTACCGCCTACTGCCAGAATATTGGTCAGTTTCGACGGCAATACGGCATCTGCCGCGCCGCGTTTTTg                                                                                                                                                >  2:358326/1‑114 (MQ=255)                                   gCAAAGCTGCCCCAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGTTACCGCCTACTGCCAGAATATTGGTCAGTTTCGACGGCAATACGGCATCTGCCGCGCCGCGTTTTTGc                                                                                                                                               >  2:386609/1‑115 (MQ=255)                                   gCAAAGCTGCCCCAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGTTACCGCCTACTGCCAGAATATTGGTCAGTTTCGACGGCAATACGGCATCTGCCGCGCCGCGTTTTTGc                                                                                                                                               <  1:386609/115‑1 (MQ=255)                                             cccAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGTTACCGCCTACTGCCAGAATATTGGTCAGTTTCGACGGCAATACGGCATCTGCCGCGCCGCGTTTTTGCACCACCAGATGGCAATCGCCCATCTTCAGCAGTGCGGGTAACGc                                                                                                   >  5:209708/1‑149 (MQ=255)                                             cccAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGTTACCGCCTACTGCCAGAATATTGGTCAGTTTCGACGGCAATACGGCATCTGCCGCGCCGCGTTTTTGCACCACCAGATGGCAATCGCCCATCTTCAGCAGTGCGGGTAAAGc                                                                                                   <  8:385219/149‑1 (MQ=255)                                              ccAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGTTACCGCCTACTGCCAGAATATTGGTCAGTTTCGACGGCAATACGGCATCTGCCGCGCCGCGTTTTTGCACCACCAGATGGCAATCGCCCATCTTCAGCAGTGCGGGTAAAGCg                                                                                                  >  1:124683/1‑149 (MQ=255)                                               cAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGTTACCGCCTACTGCCAGAATATTGGTCAGTTTCGACGGCAATACGGCATCTGCCGCGCCGCGTTTTTGCACCACCAGATGGCAATCGCCCATCTTCAGCAGTGCgg                                                                                                          <  6:82295/140‑1 (MQ=255)                                               cAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGTTACCGCCTACTGCCAGAATATTGGTCAGTTTCGACGGCAATACGGCATCTGCCGCGCCGCGTTTTTGCACCACCAGATGGCAATCGCCCATCTTCAGCAGTGCgg                                                                                                          >  5:82295/1‑140 (MQ=255)                                               cAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGTTACCGCCTACTGCCAGAATATTGGTCAGTTTCGACGGCAATACGGCATCTGCCGCGCCGCGTTTTTGCACCACCAGATGGCAATCGCCCATCTTCAGCAGTGCGGGTAAAGc                                                                                                   >  3:136508/1‑147 (MQ=255)                                               cAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGTTACCGCCTACTGCCAGAATATTGGTCAGTTTCGACGGCAATACGGCATCTGCCGCGCCGCGTTTTTGCACCACCAGATGGCAATCGCCCATCTTCAGCAGTGCGGGTAAAGc                                                                                                   <  4:136508/147‑1 (MQ=255)                                                         tAGGCTTCAGCAGTAATCACCGCGTTACCGCCTACTGCCAGAATATTGGTCAGTTTCGACGGCAATACGGCATCTGCCGCGCCGCGTTTTTGCACCACCAGATGGCAATCGCCCATCTTCAGCAGTGc                                                                                                            >  6:190689/1‑128 (MQ=255)                                                         tAGGCTTCAGCAGTAATCACCGCGTTACCGCCTACTGCCAGAATATTGGTCAGTTTCGACGGCAATACGGCATCTGCCGCGCCGCGTTTTTGCACCACCAGATGGCAATCGCCCATCTTCAGCAGTGc                                                                                                            <  5:190689/128‑1 (MQ=255)                                                              cctagcagTAATCACCGCGTTACCGCCTACTGCCAGAATATTGGTCAGTTTCGACGGCAATACGGCATCTGCCGCGCCGCGTTTTTGCACCACCAGATGGCAATCGCCCATCTTCAGCAGTGCGGGTAAAGCGTCATACGATTGCAGCg                                                                                  <  6:118263/146‑1 (MQ=255)                                                                cagcagTAATCACCGCGTTACCGCCTACTGCCAGAATATTGGTCAGTTTCGACGGCAATACGGCATCTGCCGCGCCGCGTTTTTGCACCACCAGATGGCAATCGCCCATCTTCAGCAGTGCGGGTAAAGCGTCATACGATTGCAGCg                                                                                  <  1:382222/147‑1 (MQ=255)                                                                        aTCACCGCGTTACCGCCTACTGCCAGAATATTGGTCAGTTTCGACGGCAATACGGCATCTGCCGCGCCGCGTTTTTGCACCACCAGATGGCAATCGCCCATCTTCAGCAGTGCGGGTAAAGCGTCATACGATTGCAGCGGGAAAAAtt                                                                         <  2:124683/148‑1 (MQ=255)                                                                            ccGCGTTACCGCCTACTGCCAGAATATTGGTCAGTTTCGACGGCAATACGGCATCTGCCGCGCCGCGTTTTTGCACCACCAGATgg                                                                                                                                   >  7:293227/1‑86 (MQ=255)                                                                            ccGCGTTACCGCCTACTGCCAGAATATTGGTCAGTTTCGACGGCAATACGGCATCTGCCGCGCCGCGTTTTTGCACCACCAGATgg                                                                                                                                   <  8:293227/86‑1 (MQ=255)                                                                                      gCCTACTGCCAGAATATTGGTCAGTTTCGACGGCAATACGGCATCTGCCGCGCCGCGTTTTTGCACCACCAGATGGCAATCGCCCATCTTCAGCAGTGCGGGTAAAGCGTCATACGATTGCAGCGGGAAAAATTGc                                                                       >  4:280310/1‑136 (MQ=255)                                                                                      gCCTACTGCCAGAATATTGGTCAGTTTCGACGGCAATACGGCATCTGCCGCGCCGCGTTTTTGCACCACCAGATGGCAATCGCCCATCTTCAGCAGTGCGGGTAAAGCGTCATACGATTGCAGCGGGAAAAATTGc                                                                       <  3:280310/136‑1 (MQ=255)                                                                                      gCCTACTGCCAGAATATTGGTCAGTTTCGACGGCAATACGGCATCTGCCGCGCCGCGTTTTTGCACCACCAGATGGCAATCGCCCATCTTCAGCAGTGCGGGTAAAGCGTCATACGATTGCAGCGGGAAAAATTGCATGTTGCGCAGTc                                                          >  2:262931/1‑149 (MQ=255)                                                                                      gCCTACTGCCAGAATATTGGTCAGTTTCGACGGCAATACGGCATCTGCCGCGCCGCGTTTTTGCACCACCAGATGGCAATCGCCCATCTTCAGCAGTGCGGGTAAAGCGTCATACGATTGCAGCGGGAAAAATTGCATGTTGCGCAGTc                                                          >  6:55165/1‑149 (MQ=255)                                                                                        cTACTGCCAGAATATTGGTCAGTTTCGACGGCAATACGGCATCTGCCGCGCCGCGTTTTTGCACCACCAGATGGCAATCGCCCATCTTCAGCAGTGCGGGTAAAGCGTCATACGATTGCAGCGGGAAAAATTGCATGTTGCGCAGTCCa                                                        >  7:224176/1‑149 (MQ=255)                                                                                           cTGCCAGAATATTGGTCAGTTTCGACGGCAATACGGCATCTGCCGCGCCGCGTTTTTGCACCACCAGATGGCAATCGCCCATCTTCAGCAGTGCGGGTAAAGCGTCATACGATTGCAGCGGGAAAAATTGCATGTTGCGCAGTCCACgc                                                     <  1:262931/149‑1 (MQ=255)                                                                                            tGCCAGAATATTGGTCAGTTTCGACGGCAAAACGGCGACTGCAGCGCCGCGTTTTTGCACCACCATATCGCAATCGCCCATCTACACCAGAGCCGCTAAAGCCTCATACGATAGCAGAGGGAAAAATTGCATGTTGAGCAGTCCACgct                                                    >  4:276425/1‑149 (MQ=255)                                                                                              ccAGAATATTGGTCAGTTTCGACGGCAATACGGCATCTGCCGCGCCGCGTTTTTGCACCACCAGATGGCAATCGCCCATCTTCAGCAGTGCGGGTAAAGCGTCATACGATTGCAGCGGGAAAAATTGCATGTTGCGCAGTCCACgctgc                                                  >  4:249108/1‑149 (MQ=255)                                                                                                     atTGGTCAGTTTCGACGGCAATACGGCATCTGCCGCGCCGCGTTTTTGCACCACCAGATGGCAATCGCCCATCTTCAGCAGTGCGGGTAAAGCGTCATACGATTGCAGCGGGAAAAATCGCATGTTGCGCAGTCCACGCTGCTGCGCCa                                           >  5:164759/1‑149 (MQ=255)                                                                                                        ggTCAGTTTCGACGGCAATACGGCATCTGCCGCGCCGCGTTTTTGCACCACCAGATGGCAATCGCCCATCTTCAGCAGTGCGGGTAAAGCGTCATACGATTGCAGCGGGAAAAATTGCATGTTGCGCAGTCCACGCTGCTGCGCCAttt                                        >  1:354943/1‑149 (MQ=255)                                                                                                        ggTCAGTTTCGACGGCAATACGGCATCTGCCGCGCCGCGTTTTTGCACCACCAGATGGCAATCGCCCATCTTCAGCAGTGCGGGTAAAGCGTCATACGATTGCAGCGGGAAAAATTGCATGTTGCGCAGTCCACGCTGCTGCGCCAttt                                        >  3:192432/1‑149 (MQ=255)                                                                                                                  gACGGCAATACGGCATCTGCCGCGCCGCGTTTTTGCACCACCAGATGGCAATCGCCCATCTTCAGCAGTGCGGGTAAAGCGTCATACGATTGCa                                                                                     <  2:22125/94‑1 (MQ=255)                                                                                                                  gACGGCAATACGGCATCTGCCGCGCCGCGTTTTTGCACCACCAGATGGCAATCGCCCATCTTCAGCAGTGCGGGTAAAGCGTCATACGATTGCa                                                                                     >  1:22125/1‑94 (MQ=255)                                                                                                                        aaTACGGCATCTGCCGCGCCGCGTTTTTGCACCACCAGATGGCAATCGCCCATCTTCAGCAGTGCGGGTAAAGCGTCATACGATTGCAGCGGGAAAAATTGCATGTTGCGCAGTCCACGCTGCTGCGCCATTTTTTCCAGCCGCGCttt                        >  6:15717/1‑149 (MQ=255)                                                                                                                        aaTACGGCATCTGCCGCGCCGCGTTTTTGCACCACCAGATGGCAATCGCCCATCTTCAGCAGTGCGGGTAAAGCGTCATACGATTGCAGCGGGAAAAATTGCATGTTGCGCAGTCCACGCTGCTGCGCCATTTTTTCCAGCCGCGCttt                        >  1:184071/1‑149 (MQ=255)                                                                                                                        aaTACGGCATCTGCCGCGCCGCGTTTTTGCACCACCAGATGGCAATCGCCCATCTTCAGCAGTGCGGGTAAAGCGTCATACGATTGCAGCGGGAAAAATTGCATGTTGCGCAGTCCACGCTGCTGCGCCATTTTTTCCAGCCGCGCttt                        >  7:244501/1‑149 (MQ=255)                                                                                                                            cGGCATCTGCCGCGCCGCGTTTTTGCACCACCAGATGGCAATCGCCCATCTTCAGCAGTGCGGGTAAAGCGTCATACGATTGCAGCGGGAAAAATTGCATGTTGCGCAGTCCACGCTGCTGCGCCATTTTTTCCAGCCGCGCTTTgcc                     <  5:15717/148‑1 (MQ=255)                                                                                                                             ggCATCTGCCGCGCCTCGTTTTTGTACCACCAGATGGCAATCGCCCATCTTCAGCAGTGCGGGTAAAGCGTCATACGATTGAAGCGGGAAAAATTGCATGTTTAGCAGTCAACGATGATGAGCCAGTTTTGCCAGCCGCGCTTTgccgc                   >  8:85540/1‑149 (MQ=255)                                                                                                                                    gccgcgccgcGTTTTTGCACCACCAGATGGCAATCGCCCATCTTCAGCAGTGCGGGTAAAGCGTCATACGATTGCAGCGGGAAAAATTGCATGTTGCGCAGTc                                                          <  6:89551/103‑1 (MQ=255)                                                                                                                                    gccgcgccgcGTTTTTGCACCACCAGATGGCAATCGCCCATCTTCAGCAGTGCGGGTAAAGCGTCATACGATTGCAGCGGGAAAAATTGCATGTTGCGCAGTc                                                          >  5:89551/1‑103 (MQ=255)                                                                                                                                     ccgcgccgcgTTTTTGCACCACCAGATGGCAATCGCCCATCTTCAGCAGTGCGGGTAAAGCGTCATACGATTGCAGCGGGAAAAATTGCATGTTGCGCAGTCCACGCTGCTGCGCCATTTTTTCCAGCCGCGCTTTGCCGCCGCCCt             <  3:249108/147‑1 (MQ=255)                                                                                                                                       gcgccgcgTTTTTGCACCACCAGATGGCAATCGCCCATCTTCAGCAGTGCGGGTAAAGCGTCATACGATTGCAGCGGGAAAAATTGCATGTTGCGCAGTCCACGCTGCTGCGCCATTTTTTCCAGCCGCGCTTTGCCGCCGCCCTGccc         <  8:224176/149‑1 (MQ=255)                                                                                                                                           cgcgTTTTTGCACCACCAGATGGCAATCGCCCATCTTCAGCAGTGCGGGTAAAGCGTCATACGATTGCAGCGGGAAAAATTGCATGTTGCGCAGTCCACGCTGCTGCGCCATTTTTTCCAGCCGCGCTTTGCCGCCGCCCTGCCCGAc      >  7:337804/1‑148 (MQ=255)                                                                                                                                           cgcgTTTTTGCACCACCAGATGGCAATCGCCCATCTTCAGCAGTGCGGGTAAAGCGTCATACGATTGCAGCGGGAAAAATTGCATGTTGCGCAGTCCACGCTGCTGCGCCATTTTTTCCAGCCGCGCTTTGCCGCCGCCCTGCCCGACa     <  6:209708/149‑1 (MQ=255)                                                                                                                                              gTTTTTGCACCACCAGATGGCAATCGCCCATCTTCAGCAGTGCGGGTAAAGCGTCATACGATTGCAGCGGGAAAAATTGCATGTTGCGCAGTCCACGCTGCTGCGCCATTTTTTCCAGCCGCGCTTTGCCGCCGCCCTGCCCGACAATg  >  3:342684/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                 GTTCAACGCAAACCGCAATGCCCGGAAAGGTTTCGCAAAGCTGCCCCAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGTTACCGCCTACTGCCAGAATATTGGTCAGTTTCGACGGCAATACGGCATCTGCCGCGCCGCGTTTTTGCACCACCAGATGGCAATCGCCCATCTTCAGCAGTGCGGGTAAAGCGTCATACGATTGCAGCGGGAAAAATTGCATGTTGCGCAGTCCACGCTGCTGCGCCATTTTTTCCAGCCGCGCTTTGCCGCCGCCCTGCCCGACAATG  >  NC_000913/2124666‑2124955

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑