Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,269,314	A→T	6.6%	L342H (CTC→CAC) ‡	radA ←	DNA repair enzyme

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,269,314	0	A	T	6.6%	220.0 / 3.7	75	L342H (CTC→CAC) ‡	radA	DNA repair enzyme
Reads supporting (aligned to +/- strand):  ref base A (34/36);  new base T (3/2);  total (37/38)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-01

AAACCATGGGCCAGTGGCCGGTTGCGCAAGCTCGACATCACCGCCGCCAGCAGGGCCAGGTCCGAAGCCGTTTCCAGCACCTTCACCCCGCCGACCACGTTGAGAAACACATCCTGGTCATGGGTCGGGATACCGCCATGGCGGTGGAGCACCGCCAACAACATGGCCAGACGGTTCTGGTCCAGGCCCAAGGTTACCCGCCGCGGGTTGGCCAGATGGCTATCGTCGACCAGCGCCTGCACCTCGACCAGCATTGGCCGGGTACCCTCCCAGGTTGCCATCACCACGC  >  NC_002947/5269167‑5269455                                                                                                                                                    |                                                                                                                                              aaaCCATGGGCCAGTGGCCGGTTGCGCAAGCTCGACATCACCGCCGCCAGCAGGGCCAGGTCCGAAGCCGTTTCCAGCACCTTCACCCCGCCGACCACGTTGAGAAACACATCCTGGTCATGGGTCGGGATACCGCCATGGCGGTGGAg                                                                                                                                              >  7:26752/1‑149 (MQ=255)    ccATGGGCCAGTGGCCGGTTGCGCAAGCTCGACATCACCGCCGCCAGCAGGGCCAGGTCCGAAGCCGTTTCCAGCACCTTCACCCCGCCGACCACGTTGAGAAACACATCCTGGTCATGGGTCGGGATACCGCCATGGCGGTGGAGCAc                                                                                                                                           >  5:121121/1‑149 (MQ=255)          gCCAGTGGCCGGTTGCGCAAGCTCGACATCCCCGCCGCCAGCAGGGCCAGGTCCGAAGCCGTTTCCAGCACCTTCACCCCGCCGACCACGTTGAGAAACACATCCTGGTCATGGGTCGGGATACCGCCATGGCGGTGCTCCACCGCcat                                                                                                                                     <  8:57528/149‑2 (MQ=255)           ccAGTGGCCGGTTGCGCAAGCTCGACATCACCGCCGCCAGCAGGGCCAGGTCCGAAGCCGTTTCCAGCACCTTCACCCCGCCGACCACGTTGAGAAACACATCCTGGTCATGGGTCGGGATACCGCCATGGCGGTGCTCCACCGCcatg                                                                                                                                    <  4:61724/149‑3 (MQ=255)            cAGTGGCCGGTTGCGCAAGCTCGACATCACCGCCGCCAGCAGGGCCAGGTCCGAAGCCGTTTCCAGCACCTTCACCCCGCCGACCACGTTGAGAAACACATCCTGGTCATGGGTCGGGATACCGCCATGGCGGTGGAGCACCGCcaaca                                                                                                                                   <  6:121121/149‑1 (MQ=255)              gTGGCCGGTTGCGCAAGCTCGACATCACCGCCGCCAGCAGGGCCAGGTCCGAAGCCGTTTCCAGCACCTTCACCCCGCCGACCACGTTGAGAAACACATCCTGGTCATGGGTCGGGATACCGCCATGGCGGTGGAGCACCGCCaacaac                                                                                                                                 <  1:317895/149‑1 (MQ=255)              gTGGCCGGTTGCGCAAGCTCGACATCACCGCCGCCAGCAGGGCCAGGTCCGAAGCCGTTTCCAGCACCTTCACCCCGCCGACCACGGTGAGAAACACATCCTGGTCATGGGTCGGGATACCGCCATGGCGGTGGAGCACCGCcaacaa                                                                                                                                  >  2:317895/1‑148 (MQ=255)               tGGCCGGTTGCGCAAGCTCGACATCACCGCCGCCAGCAGGGCCAGGTCCGAAGCCGTTTCCAGCACCTTCACCCCGCCGACCACGTTGAGAAACACATCCTGGTCATGGGTCGGGATACCGCCATGGCGGTGGAGCACCGCCAacaaca                                                                                                                                <  6:126163/149‑1 (MQ=255)                        gcgcAAGCTCGACATCACCGCCGCCAGCAGGGCCAGGTCCGAAGCCGTTTCCAGCACCTTCACCCCGCCGACCACGTTGAGAAACACATCCTGGTCATGGGTCGGGATACCGCCATGGCGGTGGAGCACCGCCAACAACATGGCCAGa                                                                                                                        <  1:234928/148‑1 (MQ=255)                          gcAAGCTCGACATCACCGCCGCCAGCAGGGCCAGGTCCGAAGCCGTTTCCAGCACCTTCACCCCGCCGACCACGTTGAGAAACACATCCTGGTCATGGGTCGGGATACCGCCATGGCGGTGGAGCACCGCCAACAACATGGCCAGACgg                                                                                                                     >  4:121339/1‑149 (MQ=255)                          gcAAGCTCGACATCACCGCCGCCAGCAGGGCCAGGTCCGAAGCCGTTTCCAGCACCTTCACCCCGCCGACCACGTTGAGAAACACATCCTGGTCATGGGTCGGGATACCGCCATGGCGGTGGAGCACCGCCAACAACATGGCCAGACgg                                                                                                                     >  6:145015/1‑149 (MQ=255)                          gcAAGCTCGACATCACCGCCGCCAGCAGGGCCAGGTCCGAAGCCGTTTCCAGCACCTTCACCCCGCCGACCACGTTGAGAAACACATCCTGGTCATGGGTCGGGATACCGCCATGGCGGTGGAGCACCGCCAACAACATGGCCAGACgg                                                                                                                     >  6:8057/1‑149 (MQ=255)                            aaGCTCGACATCACCGCCGCCAGCAGGGCCAGGTCCGAAGCCGTTTCCAGCACCTTCACCCCGCCGACCACGTTGAGAAACACATCCTGGTCATGGGTCGGGATACCGCCATGGCGGTGGAGCACCGCCAACAACATGGCCAGACGGtt                                                                                                                   <  3:116867/149‑1 (MQ=255)                             aGCTCGACATCACCGCCGCCAGCAGGGCCAGGTCCGAAGCCGTTTCCAGCACCTTCACCCCGCCGACCACGTTGAGAAACACATCCTGGTCATGGGTCGGGATACCGCCATGGCGGTGGAGCACCGCCAACAACATGGCCAGACGGTTc                                                                                                                  <  3:121339/149‑1 (MQ=255)                              gCTCGACATCACCGCCGCCAGCAGGGCCAGGTCCGAAGCCGTTTCCAGCACCTTCACCCCGCCGACCACGTTGAGAAACACATCCTGGTCATGGGTCGGGATACCGCCATGGCGGTGGAGCACCGCCAACAACATGGCCAGACGGTTCt                                                                                                                 >  5:82016/1‑149 (MQ=255)                               cTCGACATCACCGCCGCCAGCAGGGCCAGGTCCGAAGCCGTTTCCAGCACCTTCACCCCGCCGACCACGTTGAGAAACACATCCTGGTCATGGGTCGGGATACCGCCATGGCGGTGGAGCACCGCCACCAACATGGCCAGACGGTTCTg                                                                                                                >  6:224716/1‑149 (MQ=255)                               cTCGACATCACCGCCGCCAGCAGGGCCAGGTCCGAAGCCGTTTCCAGCACCTTCACCCCGCCGACCACGTTGAGAAACACATCCTGGTCATGGGTCGGGATACCGCCATGGCGGTGGAGCACCGCCAACAACATGGCCAGACGGTTCTg                                                                                                                <  6:82016/149‑1 (MQ=255)                                 cGACATCACCGCCGCCAGCAGGGCCAGGTCCGAAGCCGTTTCCAGCAACTTCACCCCGCCGACCACGTTGAGAAACACATCCTGGTCATGGGTCGGGATACCGACATGGCGGTGGAGCACCGCCAACAAGATGGCCAGACGGTTCTGGt                                                                                                              >  8:49631/1‑149 (MQ=255)                                  gACAACACCGCCGCCAGCAGGGCCAGGTCCGAAGCCGTTTCCAGCACCTTCACCCCGCCGACCACGTTGAGAAACACATCCTGGTCATGGGTCGGGATACCGCCATGGCGGTGGAGCACCGCCAACAACATGGCCAGACGGTTCTgg                                                                                                               <  6:125770/147‑1 (MQ=255)                                     aTCACCGCCGCCAGCAGGGCCAGGTCCGAAGCCGTTTCCAGCACCTTCACCCCGCCGACCACGTTGAGAAACACATCCTGGTCATGGGTCGGGATACCGCCATGGCGGTGGAGCACCGCCAACAACATGGCCAGACGGTTCTGGTCCAg                                                                                                          >  5:290210/1‑149 (MQ=255)                                     aTCACCGCCGCCAGCAGGGCCAGGTCCGAAGCCGTTTCCAGCACCTTCACCCCGCCGACCACGTTGAGAAACACATCCTGGTCATGGGTCGGGATACCGCCATGGCGGTGGAGCACCGCCAACAACATGGCCAGACGGTTCTGGTCCAg                                                                                                          >  3:47995/1‑149 (MQ=255)                                        aCCGCCGCCTGCAGGGCCAGGTCCGAAGCCGTTTCCAGCACCTTCACCCCGCCGACCACGTTGAGAAACACATCCTGGTCATGGGTCGGGATACCGCCATGGCGGTGGAGCACCGCCAACAACATGGCCAGACGGTTCTGGTCCAGGc                                                                                                        <  7:130870/148‑1 (MQ=255)                                                        ccAGGTCCGAAGCCGTTTCCAGCACCTTCACCCCGCCGACCACGTTGAGAAACACATCCTGGTCATGGGTCGGGATACCGCCATGGCGGTGGAGCACCGCCAACAACATGGCCAGACGGTTCTGGTCCAGGCCCAAgg                                                                                                  <  6:207820/138‑1 (MQ=255)                                                        ccAGGTCCGAAGCCGTTTCCAGCACCTTCACCCCGCCGACCACGTTGAGAAACACATCCTGGTCATGGGTCGGGATACCGCCATGGCGGTGGAGCACCGCCAACAACATGGCCAGACGGTTCTGGTCCAGGCCCAAgg                                                                                                  >  5:207820/1‑138 (MQ=255)                                                          aGGTCCGAAGCCGTTTCCAGCACCTTCACCCCGCCGACCACGTTGAGAAACACATCCTGGTCATGGGTCGGGATACCGCCATGGCGGTGGAGCACCGCCAACAACATGGCCAGACGGTTCTGGTCCAGGCCCAAGGTTACCCGCCGCgg                                                                                     <  5:8057/149‑1 (MQ=255)                                                             tCCGAAGCCGTTTCCAGCACCTTCACCCCGCCGACCCCGTTGAGAAACACATGCTGGTCATGGGTCGGGATACCGCCATGGCGGTGGAGCACCGCCAACAACATGGCCAGACGGTTCTGGTCCAGGCCCAAGGTTACCCGCCGCGGGtt                                                                                  <  8:91/149‑1 (MQ=255)                                                             tCCGAAGCCGTTTCCAGCACCTTCACCCCGCCGACCACGTTGAGAAACACATCCTGGTCATGGGTCGGGATACCGCCATGGCGGTGGAGCACCGCCAACAACATGGCCAGACGGTTCTGGTCCAGGCCCAAGGTTACCCGCCGCGGGtt                                                                                  <  8:104962/149‑1 (MQ=255)                                                               cGAAGCCGTTTCCAGCACCTTCACCCCGCCGACCACGTTGAGAAACACATCCTGGTCATGGGTCGGGATACCGCCATGGCGGTGGAGCACCGCCAACAACATGGCCAGACGGTTCTGGTCCAGGCCCAAGGTTACCCGCCGCGGGTTgg                                                                                <  6:96093/149‑1 (MQ=255)                                                                gAAGCCGTTTCCAGCACCTTCACCCCGCCGACCACGTTGAGAAACACATCCTGGTCATGGGTCGGGATACCGCCATGGCGGTGGAGCACCGCCAACAACATGGCCAGACGGtt                                                                                                                   >  5:103233/1‑113 (MQ=255)                                                                gAAGCCGTTTCCAGCACCTTCACCCCGCCGACCACGTTGAGAAACACATCCTGGTCATGGGTCGGGATACCGCCATGGCGGTGGAGCACCGCCAACAACATGGCCAGACGGtt                                                                                                                   <  6:103233/113‑1 (MQ=255)                                                                      gTTTCCAGCACCTTCACCCCGCCGACCACGTTGAGAAACACATCCTGGTCATGGGTCGGGATACCGCCATGGCGGTGGAGCACCGCCAACAACATGGCCAGACGGTTCTGGTCCAGGCCCAAGGTTACCCGCCGCGGGTTGGCCAGATg                                                                         <  7:60818/149‑1 (MQ=255)                                                                           cAGCACCTTCACCCCGCCGACCACGTTGAGAAACACATCCTGGTCATGGGTCGGGATACCGCCATGGCGGTGGAGCACCGCCAACAACATGGCCAGACGGTTCTGGTCCAGGCCCAAGGTTACCCGCCGCGGGTTGGc                                                                               <  5:7350/138‑1 (MQ=255)                                                                           cAGCACCTTCACCCCGCCGACCACGTTGAGAAACACATCCTGGTCATGGGTCGGGATACCGCCATGGCGGTGGAGCACCGCCAACAACATGGCCAGACGGTTCTGGTCCAGGCCCAAGGTTACCCGCCGCGGGTTGGc                                                                               >  6:7350/1‑138 (MQ=255)                                                                           cAGCACCTTCACCCCGCCGACCACGTTGAGAAACACATCCTGGTCATGGGTCGGGATACCGCCATGGCGGTGGAGCACCGCCAACAACATGGCCAGACGGTTCTGGTCCAGGCCCAAGGTTACCCGCCGCGGGTTGGCCAGATGGCTa                                                                     <  1:329741/148‑1 (MQ=255)                                                                                 cTTCACCCCGCCGACCACGTTGAGAAACACATCCTGGTCATGGGTCGGGATACCGCCATGGCGGTGGAGCACCGCCAACAACATGGCCAGACGGTTCTGGTCCAGGCCCAAGGTTACCCGCCGCGGGTTGGCCAGATGGCTATCGTCGa                                                              <  1:5203/149‑1 (MQ=255)                                                                                 cTTCACCCCGCCGACCACGTTGAGAAACACATCCTGGTCATGGGTCGGGATACCGCCATGGCGGTGGAGCACCGCCAACAACATGGCCAGACGGTTCTGGTCCAGGCCCAAGGTTACCCGCCGCGGGTTGGCCAGATGGCTATCGTCGa                                                              >  7:210604/1‑149 (MQ=255)                                                                                 cTTCACCCCGCCGACCACGTTGAGAAACACATCCTGGTCATGGCTCGGGATACCGCCATGGCGGTGGAGCACCGCCAACAACATGGCCAGACGGTTCTGGTCCAGGCCCAAGGTTACCCGCCGCGGGTTGGCCAGATGGCTATCGTCGa                                                              >  7:291566/1‑149 (MQ=255)                                                                                     aCCCCGCCGACCACGTTGAGAAACACATCCTGGTCATGGGTCGGGATACCGCCATGGCGGTGGAGCACCGCCAACAACATGGCCAGACGGTTCTGGTCCAGGCCCAAGGTTACCCGCCGCGGGTTGGCCAGATGGCTATCGTCGACCAg                                                          >  7:115477/1‑149 (MQ=255)                                                                                      ccccGCCGACCACGTTGAGAAACACATCCTGGTCATGGGTCGGGATACCGCCATGGCGGTGGAGCACCGCCAACAACATGGCCAGACGGTTCTGGTCCAGGCCCAAGGTTACCCGCCGCggg                                                                                    >  3:161455/1‑122 (MQ=255)                                                                                      ccccGCCGACCACGTTGAGAAACACATCCTGGTCATGGGTCGGGATACCGCCATGGCGGTGGAGCACCGCCAACAACATGGCCAGACGGTTCTGGTCCAGGCCCAAGGTTACCCGCCGCggg                                                                                    <  4:161455/122‑1 (MQ=255)                                                                                         cgccgACCACGTTGAGAAACACATCCTGGTCATGGGTCGGGATACCGCCATGGCGGTGGAGCACCGCCAACAACATGGCCAGACGGTTCTGGTCCAGGCCCAAGGTTACCCGCCGCGGGTTGGCCAGATGGCTATCGTCGACCAGCGcc                                                      >  7:173202/1‑149 (MQ=255)                                                                                            cgACCACGTTGAGAAACACATCCTGGTCATGGGTCGGGATACCGCCATGGCGGTGGAGCACCGCCAACAACATGGCCAGACGGTTCTGGTCCAGGCCCAAGGTTACCCGCCGCGGGTTGGCCAGATGGCTATCGTCGACCAGCGCCTg                                                    >  3:43428/1‑148 (MQ=255)                                                                                               ccACGTTGAGAAACACATCCTGGTCATGGCTCGGGATACCGCCATGGCGGTGGAGCACCGCCAACAACATGGCCAGACGGTTCTGGTCCAGGCCCAATGTTACCCGCAGCGCGTTGGCCACATGGCTATCGTCGACCAGCGCCTGCAc                                                 >  1:342001/1‑148 (MQ=255)                                                                                                 aCGTTGAGAAACACATCCTGGTCATGGGTCGGGATACCGCCATGGCGGTGGAGCACCGCCAACAACATGGCCAGACGGTTCTGGTCCAGGCCCAAGGTTACCCGCCGCGGGTTGGCCAGATGGCTATCGTCGACCAGCGCCTGCACCTc                                              <  8:115477/149‑1 (MQ=255)                                                                                                    ttGAGAAACACATCCTGGTCATGGGTCGGGATACCGCCATGGCGGTGGATCACCGCCAACAACATGGCCAGACGGTTCTGGTCCAGGCCCAAGGTTACCCGCCGCGGGTTGGCCAGATGGCTATCGTCGACCAGCGCCTGCACCTCGAc                                           <  8:160971/149‑1 (MQ=255)                                                                                                    ttGAGAAACACATCCTGGTCATGGGTCGGGATACCGCCATGGCGGTGGAGCACCGCCAACAACATGGCCAGACGGTTCTGGTCCAGGCCCAAGGTTACCCGCCGCGGGTTGGCCAGATGGCTATCGTCGACCAGCGCCTGCACCTCGAc                                           <  8:210604/149‑1 (MQ=255)                                                                                                           acacATCCTGGTCATGGGTCGGGATACCGCCATGGCGGTGGAGCACCGCCAACAACATGGCCAGACGGTTCTGGTCCAGGCCCAAGGTTACCcgccgc                                                                                       <  7:298091/98‑1 (MQ=255)                                                                                                           acacATCCTGGTCATGGGTCGGGATACCGCCATGGCGGTGGAGCACCGCCAACAACATGGCCAGACGGTTCTGGTCCAGGCCCAAGGTTACCcgccgc                                                                                       >  8:298091/1‑98 (MQ=255)                                                                                                           acacATCCTGGTCATGGGTCGGGATACCGCCATGGCGGTGGAGCACCGCCAACAACATGGCCAGACGGTTCTGGTCCAGGCCCAAGGTTACCCGCCGCGGGTTGGCCAGATGGCTATCGTCGACCAGCGCCTGCACCTCGACCAGCAtt                                    <  8:173202/149‑1 (MQ=255)                                                                                                               aTCCTGGTCATGGGTCGGGATACCGCCATGGCGGTGGAGCACCGCCAACAACATGGCCAGACGGTTCTGGTCCAGGCCCAAGGTTACCCGCCGCGGGTTGGCCAGATGGCTAtc                                                                   >  8:262618/1‑114 (MQ=255)                                                                                                               aTCCTGGTCATGGGTCGGGATACCGCCATGGCGGTGGAGCACCGCCAACAACATGGCCAGACGGTTCTGGTCCAGGCCCAAGGTTACCCGCCGCGGGTTGGCCAGATGGCTAtc                                                                   <  7:262618/114‑1 (MQ=255)                                                                                                                tCCTGGTCATGGGTTGGGATACCGCCATGGCGGTGGAGCACCGCCAACAACATGGCCAGACGGTTCTGGTCCAGGCCCAAGGTTACCCGCCGCGGGTTGGCCAGATGGCTATCGTCGACCAGCGCCTGCACCTCGACCAGCATTGGcc                                >  4:156102/1‑148 (MQ=255)                                                                                                                tCCTGGTCATGGGTTGGGATACCGCCATGGCGGTGGAGCACCGCCAACAACATGGCCAGACGGTTCTGGTCCAGGCCCAAGGTTACCCGCCGCGGGTTGGCCAGATGGCTATCGTCGACCAGCGCCTGCACCTCGACCAGCATTGGcc                                <  3:156102/148‑1 (MQ=255)                                                                                                                  cTGGTCATGGGTCGGGATACCGCCATGGCGGTGGAGCACCGCCAACAACATGGCCAGACGGTTCTGGTCCAGGCCCAAGGTTACCCGCCGCGGGTTGGCCAGATGGCTATCGTCGACCAGCGCCTGCACCTCGACCAGCATTGGCCgg                              <  8:26752/148‑1 (MQ=255)                                                                                                                    ggTCATGGGTCGGGATACCGCCATGGCGGTGGAGCACCGCCAACAACATGGCCAGACGGTTCTGGTCCAGGCCCAAGGTTACCCGCCGCGGGTTGGCCAGATGGCTATCGTCGACCAGCGCCTGc                                                   >  3:289018/1‑125 (MQ=255)                                                                                                                    ggTCATGGGTCGGGATACCGCCATGGCGGTGGAGCACCGCCAACAACATGGCCAGACGGTTCTGGTCCAGGCCCAAGGTTACCCGCCGCGGGTTGGCCAGATGGCTATCGTCGACCAGCGCCTGc                                                   <  4:289018/125‑1 (MQ=255)                                                                                                                     gTCATGTGTCGGGATAGGGGCATGGCGGTGGAGCACCGCCAACTACATGGCCTGACTGTTCTGGTCCAGGCCCTAGGTTTGACGCCGCGGGTTGGCCAGATGGGTATCGTGGACCAGCGCCTGCAACTCGTGCAGCATTGGGGGGGTa                           <  2:342001/148‑1 (MQ=255)                                                                                                                     gTCATGGGTCGGGATACCGCCATGGCGGTGGAGCACCGCCAACAACATGGCCAGACGGTTCTGGTCCAGGCCCAAGGTTACCCGCCGCGGGTTGGCCAGATGGCTATCGTCGACCAGCGCCTGCACCTCGACCAGCATTGGCCGGGt                            <  5:145015/147‑1 (MQ=255)                                                                                                                     gTCATGGGTCGGGATACCGCCATGGCGGTGGAGCACCGCCAACAACATGGCCAGACGGTTCTGGTCCAGGCCCAAGGTTACCCGCCGCGGGTTGGCCAGATGGCTATCGTCGACCAGCGCCTGCACCTCGACCAGCATTGGCCGGGTAc                          >  6:62812/1‑149 (MQ=255)                                                                                                                        aTGGGTCGGGATACCGCCATGGCGGTGGAGCACCGCCAACAACATGGCCAGACGGTTCTGGTCCAGGCCCCAGGTTACCCGCCGCGGGTTGGCCAGATGGCTATCGTCGACCAGCGCCTGCACCTCGACCAGCATTGGCCGGGTACCCt                       >  8:248557/1‑149 (MQ=255)                                                                                                                        aTGGGTCGGGATACCGCCATGGCGGTGGAGCACCGCCAACAACATGGCCAGACGGTTCTGGTCCAGGCCCAAGGTTACCCGCCGCGGGTTGGCCAGATGGCTATCGTCGACCAGCGCCTGCACCTCGACCAGCATTGGCCGGGTACCCt                       >  7:179873/1‑149 (MQ=255)                                                                                                                          gggTCGGGATACCGCCATGGCGGTGGAGCACCGCCAACAACATGGCCAGACGGTTCTGGTCCAGGCCCAAGGTTACCCGCCGCGGGTTGGCCAGATGGCTATCGTCGACCAGCGCCTGCACCTCGACCAGCATTGGCCGGGTACCCTcc                     >  2:197929/1‑149 (MQ=255)                                                                                                                            gTCGGGATACCGCCATGGCGGTGGAGCACCGCCAACAACATGGCCAGACGGTTCTGGTCCAGGCCCAAGGTTACCCGCCGCGGGTTGGCCAGATGGCTATCGTCGACCAGCGCCTGCACCTCGACCAGCATTGGCCGGGTACCCTcc                     >  7:110335/1‑147 (MQ=255)                                                                                                                            gTCGGGATACCGCCATGGCGGTGGAGCACCGCCAACAACATGGCCAGACGGTTCTGGTCCAGGCCCAAGGTTACCCGCCGCGGGTTGGCCAGATGGCTATCGTCGACCAGCGCCTGCACCTCGACCAGCATTGGCCGGGTACCCTcc                     >  4:159900/1‑147 (MQ=255)                                                                                                                                   tACCGCCATGGCGGTGGAGCACCGCCAACAACATGGCCAGACGGTTCTGGTCCAGGCCCAAGGTTACCCGCCGCGGGTTGGCCAGATGGCTATCGTCGACCAGCGCCTGCACCTCGACCAGCATTGGCCGGGTACCCTCCCAGGTTGc             >  8:149093/1‑148 (MQ=255)                                                                                                                                      cGCCATGGCGGTGGAGCACCGCCAACAACATGGCCAGACGGTTCTGGTCCAGGCCCAAGGTTACCCGCCGCGGGTTGGCCAGATGGCTATCGTCGACCAGCGCCTGCACCTCGACCAGCATTGGCCGGGTACCCTCCCAGGTTGCCAt          >  6:218250/1‑148 (MQ=255)                                                                                                                                        ccATGGCGGTGGAGCACCGCCAACAACATGGCCAGACGGTTCTGGTCCAGGCCCAAGGTTACCCGCCGCGGGTTGGCCAGATGGCTATCGTCGACCAGCGCCTGCACCTCGACCAGCATTGGCCGGGTACCCTCCCAGGTTGCCATcac       >  2:302987/1‑149 (MQ=255)                                                                                                                                          ttGGCGGTGCTCCACCGCCAACAACATGGCCAGACGGTTCTGGTCCAGGCCCAAGGTTACCCGCCGCGGGTTGGCCAGATGGCTATCGTCGACCAGCGCCTGCACCTCGACCAGCATTGGCCGGGTACCCTCCCAGGTTGCCATcacca     >  3:154031/2‑149 (MQ=255)                                                                                                                                          ttGGCGGTGCTCCACCGCCAACAACATGGCCAGACGGTTCTGGTCCAGGCCCAAGGTTACCCGCCGCGGGTTGGCCAGATGGCTATCGTCGACCAGCGCCTGCACCTCGACCAGCATTGGCCGGGTACCCTCCCAGGTTGCCATcacca     >  4:161683/2‑149 (MQ=255)                                                                                                                                          aTGGCGGTGGAGCACCGCCAACAACATGGCCAGACGGTTCTGGTCCAGGCCCAAGGTTACCCGCCGCGGGTTGGCCAGATGGCTATCGTCGACCAGCGCCTGCACCTCGACCAGCATTGGCCGGGTACCCTCCCAGGTTGCCATcacca     <  5:218250/149‑1 (MQ=255)                                                                                                                                           tGGCGGTGGAGCACCGCCAACAACATGGCCAGACGGTTCTGGTCCAGGCCCAAGGTTACCCGCCGCGGGTTGGCCAGATGGCTATCGTCGACCAGCGCCTGCACCTCGACCAGCATTGGCCGGGTACCCTCCCAGGTTGCCATcaccac    <  2:110803/149‑1 (MQ=255)                                                                                                                                             gCGGTGGTCCACCGCCAACAACATGGCCAGACGGTTCTGGTCCAGGCCCAAGGTTTCCCGCCGCGGGTTGGCCAGATGGCTATCGTCGCCCAGCGCCTGCACCTCGGCCAGCATTGGCCGGGTACGCTCCCAGGTTGCCATGACCACGc  >  6:41049/1‑149 (MQ=255)                                                                                                                                             gCGGTGGAGCACCGCCAACAACATGGCCAGACGGTTCTGGTCCAGGCCCAAGGTTACCCGCCGCGGGTTGGCCAGATGGCTATCGTCGACCAGCGCCTGCACCTCGACCAGCATTGGCCGGGTACCCTCCCAGGTTGCCATCACCACGc  <  8:179873/149‑1 (MQ=255)                                                                                                                                             gCGGTGGAGCACCGCCAACAACATGGCCAGACGGTTCTGGTCCAGGCCCAAGGTTACCCGCCGCGGGTTGGCCAGATGGCTATCGTCGACCAGCGCCTGCACCTCGACCAGCATTGGCCGGGTACCCTCCCAGGTTGCCATCACCACGc  <  6:290210/149‑1 (MQ=255)                                                                                                                                             gCGGTGGAGCACCGCCAACAACATGGCCAGACGGTTCTGGTCCAGGCCCAAGGTTACCCGCCGCGGGTTGGCCAGATGGCTATCGTCGACCAGCGCCTGCACCTCGACCAGCATTGGCCGGGTACCCTCCCAGGTTGCCATCACCACGc  <  4:47995/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                              AAACCATGGGCCAGTGGCCGGTTGCGCAAGCTCGACATCACCGCCGCCAGCAGGGCCAGGTCCGAAGCCGTTTCCAGCACCTTCACCCCGCCGACCACGTTGAGAAACACATCCTGGTCATGGGTCGGGATACCGCCATGGCGGTGGAGCACCGCCAACAACATGGCCAGACGGTTCTGGTCCAGGCCCAAGGTTACCCGCCGCGGGTTGGCCAGATGGCTATCGTCGACCAGCGCCTGCACCTCGACCAGCATTGGCCGGGTACCCTCCCAGGTTGCCATCACCACGC  >  NC_002947/5269167‑5269455

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