Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,010,979	T→G	8.8%	W14G (TGG→GGG)	bcsG →	endoglucanase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,010,979	0	T	G	8.8%	198.0 / 8.4	68	W14G (TGG→GGG)	bcsG	endoglucanase
Reads supporting (aligned to +/- strand):  ref base T (32/30);  new base G (2/4);  total (34/34)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01

GCCTGTGTTCGGCATTACGTGGCGTGTTGCGCCGTCGCCTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATG‑CTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAAGTGGCTGCGTATT  >  NC_002947/3010838‑3011122                                                                                                                                              |                                                                                                                                                 gcCTGTGTTCGGCATTACGTGGCGTGTTGCGCCGTCGCCTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTgg                                                                                                                                           >  4:28970/1‑149 (MQ=255)    tgtgTTCGGCATTACGTGGCGTGTTGCGCCGTCGCCTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCAGTCCAGGCCAGCGGGccagtacc                                                                                                                                        >  3:40247/1‑143 (MQ=255)        ttCGGCATTACGTGGCGTGTTGCGCCGTCGCCTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTgg                                                                                                                                      >  3:180023/1‑147 (MQ=255)        ttCGGCATTACGTGGCGTGTTGCGCCGTCGCCTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCAGTCCAGGCCAGCGGGCCAGTACCGGGa                                                                                                                                    <  6:211970/149‑1 (MQ=38)        ttCGGCATTACGTGGCGTGTTGCGCCGTCGCCTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCAGTCCAGGCCAGCGGGCCAGTACCGGGa                                                                                                                                    <  4:40247/149‑1 (MQ=38)                 aCGTGGCGTGTTGCGCCGTCGCCTTCCCCCCCGTTAGCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAAc                                                                                                                           >  1:246489/1‑149 (MQ=255)                    tGGCGTGTTGCGCCGTCGCCTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACctc                                                                                                                        >  8:168752/1‑149 (MQ=255)                           ttGCGCCGTCGCCTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTAttttt                                                                                                                 <  1:48559/149‑1 (MQ=255)                               gcCGTCGCCTT‑CCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGc                                                                                                             >  6:193640/1‑148 (MQ=255)                                  gTCGCCTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCaaa                                                                                                          <  3:28970/149‑1 (MQ=255)                                     gCCTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGtt                                                                                                       >  7:326527/1‑149 (MQ=255)                                     gCCTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGtt                                                                                                       <  4:353472/149‑1 (MQ=255)                                     gCCTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGtt                                                                                                       >  2:245216/1‑149 (MQ=255)                                        ttCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTct                                                                                                    >  4:333638/1‑149 (MQ=255)                                           ccccccGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTctgc                                                                                                  <  1:245216/148‑1 (MQ=255)                                            cccccGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTctgctg                                                                                                <  7:49640/149‑1 (MQ=255)                                             ccccGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGc                                                                                                                                               <  6:104519/101‑1 (MQ=255)                                             ccccGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGc                                                                                                                                               >  5:104519/1‑101 (MQ=255)                                             ccccGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTggg                                                                                                                                     >  8:287398/1‑111 (MQ=255)                                             ccccGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTggg                                                                                                                                     <  7:287398/111‑1 (MQ=255)                                             ccccGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTgg                                                                                               <  3:333638/149‑1 (MQ=255)                                             ccccGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTgg                                                                                               >  5:102246/1‑149 (MQ=255)                                              cccGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGt                                                                                              >  1:125511/1‑149 (MQ=255)                                                      cTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTgg                                                                                      <  4:180023/149‑1 (MQ=255)                                                      cTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTgg                                                                                      >  1:257436/1‑149 (MQ=255)                                                       tGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTgg                                                                                                                                      >  3:287204/1‑100 (MQ=255)                                                       tGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTgg                                                                                                                                      <  4:287204/100‑1 (MQ=255)                                                       tGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTggg                                                                                     >  8:192244/1‑149 (MQ=255)                                                       tGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTggg                                                                                     >  8:187053/1‑149 (MQ=255)                                                        gAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTggg                                                                                     >  3:131365/1‑148 (MQ=255)                                                        gAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTggg                                                                                     >  7:29113/1‑148 (MQ=255)                                                            ccGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATG‑c                                                                               >  4:185998/1‑149 (MQ=255)                                                            ccGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATG‑c                                                                               <  7:192244/149‑1 (MQ=255)                                                            ccGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATG‑c                                                                               <  7:187053/149‑1 (MQ=255)                                                              gCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATG‑CTc                                                                             <  7:168752/149‑1 (MQ=255)                                                                   tcGTACGCCTCCGCTCTTCACGCTATGAACCCCATAATTTCTAATACCATCTTCTGCCTCCCTGTACTGGCTCTCTGGCCTGTACTTGTATGGTTGAACCTCAATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATG‑CTCTACt                                                                         <  2:332215/147‑1 (MQ=255)                                                                        cgccgccgcGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATG‑CTCGACTTCCAgg                                                                   <  3:185998/149‑1 (MQ=255)                                                                        cgccgccgcGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATG‑CTCGACTTCCAg                                                                    <  6:102246/148‑1 (MQ=255)                                                                              cgcgCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATG‑CTCGACTTCCAGGCGCTGc                                                             <  2:125511/149‑1 (MQ=255)                                                                                      aCGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATG‑CTCGACTTCCAGGCGCTGCCGAACCtg                                                     <  2:257436/149‑1 (MQ=255)                                                                                                 cccAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTg                                                                                                                                       >  3:183617/1‑57 (MQ=255)                                                                                                 cccAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTg                                                                                                                                       <  4:183617/57‑1 (MQ=255)                                                                                                          aCTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGc                                  <  1:263684/149‑1 (MQ=255)                                                                                                           cTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATG‑CTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCtt                                <  2:246489/149‑1 (MQ=255)                                                                                                            tAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATG‑CTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCt                                    <  8:225815/144‑1 (MQ=255)                                                                                                            tAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATG‑CTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCt                                    >  7:225815/1‑144 (MQ=255)                                                                                                             aaTACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATG‑CTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGc                                          <  2:7250/137‑1 (MQ=255)                                                                                                             aaTACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATG‑CTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGc                                          >  1:7250/1‑137 (MQ=255)                                                                                                              aTACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATG‑CTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGt                              <  8:29113/148‑1 (MQ=255)                                                                                                                  caccaTCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATG‑CTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGc                                  <  8:45706/140‑1 (MQ=255)                                                                                                                  caccaTCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATG‑CTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGc                                  >  7:45706/1‑140 (MQ=255)                                                                                                                  caccaTCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATG‑CTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCtt                         >  8:327755/1‑149 (MQ=255)                                                                                                                   accaTCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTc                                                                                                     >  6:271292/1‑73 (MQ=255)                                                                                                                   accaTCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTc                                                                                                     <  5:271292/73‑1 (MQ=255)                                                                                                                   accaTCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATG‑CTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGttt                                      >  4:319768/1‑135 (MQ=255)                                                                                                                   accaTCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATG‑CTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCttt                        >  3:24325/1‑149 (MQ=255)                                                                                                                      aTCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAATGTTTCTGCTGGTCGCGCTGGGCATG‑CTCGACTTCCAGGCGCTGCCGAACCTGTGGTTCGGTGCGTTTCTGCTTGTCCCTTTGcc                     >  4:231939/1‑149 (MQ=255)                                                                                                                       tCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATG‑CTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCg                    <  8:326527/149‑1 (MQ=255)                                                                                                                        cACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATG‑CTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCgg                   >  2:109273/1‑149 (MQ=255)                                                                                                                         aCCCTCCCAGTCCAGGCCAGCGGGCCAGTACCGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATG‑CTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGc                                  <  6:17893/133‑1 (MQ=37)                                                                                                                           ccTCCCAGTCCAGGCCAGCGGGCCAGTACCGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATG‑CTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCt                                            >  2:249269/1‑121 (MQ=37)                                                                                                                           ccTCCCAGTCCAGGCCAGCGGGCCAGTACCGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATG‑CTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCt                                            <  1:249269/121‑1 (MQ=37)                                                                                                                                     cTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATG‑CTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAAGTGGCTGc       <  4:24325/148‑1 (MQ=255)                                                                                                                                        gCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATG‑CTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTg                                   <  8:325656/117‑1 (MQ=255)                                                                                                                                        gCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATG‑CTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTg                                   >  7:325656/1‑117 (MQ=255)                                                                                                                                        gCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATG‑CTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCgg                   >  5:363030/1‑133 (MQ=255)                                                                                                                                        gCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATG‑CTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCgg                   <  6:363030/133‑1 (MQ=255)                                                                                                                                         cccGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATG‑CTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAAGTGGCTGCGTAtt  >  6:366824/1‑149 (MQ=255)                                                                                                                                              |                                                                                                                                                 GCCTGTGTTCGGCATTACGTGGCGTGTTGCGCCGTCGCCTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATG‑CTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAAGTGGCTGCGTATT  >  NC_002947/3010838‑3011122

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