Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,316,619	A→T	5.2%	S13C (AGC→TGC)	PP_1150 →	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,316,619	0	A	T	5.2%	237.8 / 2.0	77	S13C (AGC→TGC)	PP_1150	membrane protein
Reads supporting (aligned to +/- strand):  ref base A (38/35);  new base T (2/2);  total (40/37)
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.13e-01

GCCAATCGAATATACATGATTAGCTAGCTAATTAATCGTTTCCGCTTTATCCTGCTCATCA‑TTGTCTATCTTTTCAGCAAGACTTTTCCCCGCAAGGTCCGCCTCGTGCCCATCACCTTTCAGGCCCTGTTCGCGCCCAGCAGCCTCGCCCTCAAGTTTGCCATCAAGACCCTGCTCGGGGGCGGCCTGGCGTTGTGGCTGGCGATGCGCTGGGGGCTGGAGCAACCCTCCTGGGCATTGATGACTGCTTTCATCGTCGCCCAGCCGCTATCGGGCATGGTGGTGC  >  NC_002947/1316478‑1316763                                                                                                                                               |                                                                                                                                                 gCCAATCGAATATACATGATTAGCTAGCTAATTAATCGTTTCCGCTTTATCCTGCTCATCA‑TTGTCTATCTTTTCAGCAAGACTTTTCCCCGCAAGGTCCGCCTCGTGCCCATCACCTTTCAGGCCCTGTTCGCGCCCAGCAGCCTCGc                                                                                                                                           <  1:260040/149‑1 (MQ=255)  ccAATCGAATATACATGATTAGCTAGCTAATTAATCGTTTCCGCTTTATCCTGCTCATCA‑TTGTCTATCTTTTCAGCAAGACTTTTCCCCGCAAGGTCCGCCTCGTGCCCATCACCTTTCAGGCCCTGTTCGCGCCCAGCAGCCTCGcc                                                                                                                                          >  4:241954/1‑149 (MQ=255)      tCGAATATACATGATTAGCTAGCTAATTAATCGTTTCCGCTTTATCCTGCTCATCA‑TTGTCTATCTTTTCAGCAAGACTTTTCCCCGCAAGGTCCGCCTCGTGCCCATCACCTTTCAGGCCCTGTTCGCGCCCAGCAGCCTCGCCCTCa                                                                                                                                      <  8:44494/149‑1 (MQ=255)        gAATATTAATGATTAGCTAGCTAATTAATCGTTTACGCTTTATCCTGCTCATCA‑TTTTCTATCTTTTCAGCAAGACTTTTCCCCGCAAGGTCCGCCTCGTGCCCATCACCTTTCAGGCCCTGTTCGCGCCCAGCAGCCTCGCCCTCaa                                                                                                                                     <  8:173029/148‑1 (MQ=255)         aaTATACATGATTAGCTAGCTAATTAATCGTTTCCGCTTTATCCTGCTCATCA‑TTGTCTATCTTTTCAGCAAGACTTTTCCCCGCAAGGTCCGCCTCGTGCCCATCACCTTTCAGGCCCTGTTCGCGCCCAGCAGCCTCGCCCTCAAGt                                                                                                                                   >  7:200181/1‑149 (MQ=255)           tataCATGATTAGCTAGCTAATTAATCGTTTCCGCTTTATCCTGCTCATCA‑TTGTCTATCTTTTCAGCAAGACTTTTCCCCGCAAGGTCCGCCTCGTGCCCATCACCTTTCAGGCCCTGTTCGCGCCCAGCAGCCTCGCCCTCAAGttt                                                                                                                                 <  3:241954/149‑1 (MQ=255)            ataCATGATTAGCTAGCTAATTAATCGTTTCCGCTTTATCCTGCTCATCA‑TTGTCTATCTTTTCAGCAAGACTTTTCCCCGCAAGGTCCGCCTCGTGCCCATCACCTTTCAGGCCCTGTTCGCGCCCAGCAGCCTCGCCCTCAAGTTTg                                                                                                                                >  7:141796/1‑149 (MQ=255)              aCATGATTAGCTAGCTAATTAATCGTTTCCGCTTTATCCTGCTCATCA‑TTGTCTATCTTTTCAGCAAGACTTTTCCCCGCAAGGTCCGCCTCGTGCCCATCACCTTTCAGGCCCTGTTCGCGCCCAGCAGCCTCGCCCTCAAGTTTGc                                                                                                                               >  4:314111/1‑148 (MQ=255)               cATGATTAGCTATCTAATTAATCGTTTCCGCTTTATCCTGCTCATAA‑TTGTCTATCTTTTCAGCAAGACTTTTCCCCGCAAGGTCCGCCTCGTGCCCATCACCTTTCAGGCCCTGTTCGCGCCCAGCAGCCTCGCCCTCAAGTTTGcc                                                                                                                              <  5:176620/148‑1 (MQ=255)                 tGATTAGCTAGCTAATTAATCGTTTCCGCTTTATCCTGCTCATCA‑TTGTCTATCTTTTCAGCAAGACTTTTCCCCGCAAGGTCCGCCTCGTGCCCATCACCTTTCAGGCCCTGTTCGCGCCCAGCAGCCTCGCCCTCAAGTTTGCCATc                                                                                                                           <  8:200181/149‑1 (MQ=255)                          agctaATTAATCGTTTCCGCTTTATCCTGCTCATCA‑ATGTCTATCTTTTCAGCAACAATTTAACCAGCAAGGTCCGCCTCGTGCCCATCACCTTTCAGGCCCTGTTCCAGAAAAGCAGCCTCGCCCTCAACTTTGCCATCAAGAccat                                                                                                                   >  7:286657/1‑146 (MQ=255)                            ctaATTAATCGTTTCCGCTTTATCCTGCTCATCA‑TTGTCTATCTTTTCAGCAAGACTTTTCCCCGCAAGGTCCGCCTCGTGCCCATCACCTTTCAGGCCCTGTTCGCGCCCAGCAGCCTCGCCCTCAAGTTTGCCATCAAGACCCTGCt                                                                                                                <  4:268012/149‑1 (MQ=255)                              aattaatCGTTTCCGCTTTATCCTGCTCATCAGTTGTCTATCTTTTCAGCAAGACTTTTCCCCGCAAGGTCCGCCTCGTGCCCATCACCTTTCAGGCCCTGTTCGCGCCCAGCAGCCTCGCCCTCAAGTTTGCCATCAAGACCCTGCTc                                                                                                               >  1:45722/1‑149 (MQ=255)                              aattaatCGTTTCCGCTTTATCCTGCTCATCA‑TTGTCTATCTTTTCAGCAAGACTTTTCCCCGCAAGGTCCGCCTCGTGCCCATCACCTTTCAGGCCCTGTTCGCGCCCAGCTGCCTCGCCCTCAAGTTTGCCATCAAGACCCTGCTCg                                                                                                              >  2:9618/1‑149 (MQ=255)                              aattaatCGTTTCCGCTTTATCCTGCTCATCA‑TTGTCTATCATTTCAGCAAGACTTTTCCCCGCAAGGTCCGCCTCGTGCCCATCACCTTTCAGGCCCTGTTCGCGCCCAGCTGCCTCGCCCTCAAGTTTGCCATCAAGACCCTGCTCg                                                                                                              >  2:15664/1‑149 (MQ=255)                                  aatCGTTTCCGCTTTATCCTGCTCATCA‑TTGTCTATCTTTTCAGCAAGACTTTTCCCCGCAAGGTCCGCCTCGTGCCCATCACCTTTCAGGCCCTGTTCGCGCCCAGCAGCCTCGCCCTCAAGTTTGCCATCAAGACCCTGCTCggggg                                                                                                          >  2:41943/1‑149 (MQ=255)                                  aatCGTTTCCGCTTTATCCTGCTCATCA‑TTGTCTATCTTTTCAGCAAGACTTTTCCCCGCAAGGTCCGCCTCGTGCCCATCACCTTTCAGGCCATGTTCGCGCCCAGCAGCCTCGCCCTCAAGTTTGCCATCAAGACCCTGCTCggggg                                                                                                          >  3:42410/1‑149 (MQ=255)                                           atgtttATCCTGATCATCA‑TTGTCTATCTTTTCAGCAAGACTTTTCGCCGCAAGGTCCGCCTCGTGCCTATCACCTTTCAGTCCCTGTTCGCGCCCAGCAGCGTCTCCCTCATGTTTGCGATCAAGACCCTGCTCTGGTGCTGCCTGGc                                                                                                 <  4:196285/146‑1 (MQ=255)                                                tATCCTGCTCATCA‑TTGTCTATCTTTTCAGCAAGACTTTTCCCCGCAAGGTCCGCCTCGTGCCCATCACCTTTCAGGCCCTGTTCGCGCCCAGCAGCCTCGCCCTCAAGTTTGCCATCAAGACCCTGCTCGGGGGCGGCCTGGCGtt                                                                                              <  4:42410/147‑1 (MQ=255)                                                  tCCTGCTCATCA‑TTGTCTATCTTTTCAGCAAGACTTTTCCCCGCAAGGTCCGCCTCGTGCCCATCACCTTTCAGGCCCTGTTCGCGCCCAGCAGCCTCGCCCTCAAGTTTGCCATCAAGACCCTGCTCGGGGGCGGCCTGGCGTtgtg                                                                                           <  1:41943/148‑1 (MQ=255)                                                      gCTCATCA‑TTGTCTATCTTTTCAGCAAGACTTTTCCCCGCAAGGTCCGCCTCGTGCCCATCACCTTTCAGGCCCTGTTCGCGCCCAGCAGCCTCGCCCTCAAGTTTGCCATCAAGACCCTGCTCGGGGGCGGCCTGGCGTTGtggctgg                                                                                      >  4:239061/1‑149 (MQ=255)                                                               ttGTCTATCTTTTCAGCAAGACTTTTCCCCGCAAGGTCCGCCTCGTGCCCATCACCTTTCAGGCCCTGTTCGCGCCCAGCAGCCTCGCCCTCAAGTTTGCCATCAAGACCCTGCTCGGGGGCGGCCTGGCGTTGTGGCTGGCGATgcgc                                                                              >  8:22538/1‑149 (MQ=255)                                                               ttGTCTATCTTTTCAGCAAGACTTTTCCCCGCAAGGTCCGCCTCGTGCCCATCACCTTTCAGGCCCTGTTCGCGCCCAGCAGCCTCGCCCTCAAGTTTGCCATCAAGACCCTGCTCGGGGGCGGCCTGGCGTTGTGGCTGGCGATgcgc                                                                              <  7:22538/149‑1 (MQ=255)                                                               ttGTCTATCTTTTCAGCAAGACTTTTCCCCGCAAGGTCCGCCTCGTGCCCATCACCTTTCAGGCCCTGTTCGCGCCCAGCAGCCTCGCCCTCAAGTTTGCCATCAAGACCCTGCTCGGGGGCGGCCTGGCGTTGTGGCTGGCGATgcgc                                                                              >  4:269461/1‑149 (MQ=255)                                                               ttGTCTATCTTTTCAGCAAGACTTTTCCCCGCAAGGTCCGCCTCGTGCCCATCACCTTTCAGGCCCTGTTAGCGCCCAGCAGCCTCGCCCTCAAGTTTGCCATCAAGACCCTGCTCGGGGGCGGCCTGGCGTTGTGGCTGGCGATgcgc                                                                              >  8:325217/1‑149 (MQ=255)                                                                 gTCTATCTTTTCAGCAAGACTTTTCCCCGCAAGGTCCGCCTCGTGCCCATCACCTTTCAGGCTCTGTTCGCGCCCAGCAGCCTCGCCCTCAAGTTTGCCATCAAGACCCTGCTCGGGGGCGGCCTGGCGTTGTGGCTGGCGATGCGCTg                                                                            >  6:236147/1‑149 (MQ=255)                                                                  tCTATCTTTTCAGCAAGACTTTTCCCCGCAAGGTCCGCCTCGTGCCCATCACCTTTCAGGCCCTGTTCGCGCCCAGCAGCCTCGCCCTCAAGTTTGCCATCAAGACCCTGCTCGGGGGCGGCCTGGCGTTGTGGCTGGCGATGCGCTgg                                                                           >  7:208597/1‑149 (MQ=255)                                                                    tATCTTTTCAGCAAGACTTTTCCCCGCAAGGTCCGCCTCGTGCCCATCACCTTTCAGGCCCTGTTCGCGCCCAGCAGCCTCGCCCTCAAGTTTGCCATCAAGACCCTGCTCGGGGGCGGCCTGGCGTTGTGGCTGGCGATGCGCTgggg                                                                         >  6:330511/1‑149 (MQ=255)                                                                     aTCTTTTCAGCAAGACTTTTCCCCGCAAGGTCCGCCTCGTGCCCATCACCTTTCAGGCCCTGTTCGCGCCCAGCAGCCTCGCCCTCAAGTTTGCCATCAAGACCCTGCTCGGGGGCGGCCTGGCGTTGTGGCTGGCGATGCGCTggggg                                                                        >  6:243189/1‑149 (MQ=255)                                                                     aTCTTTTCAGCAAGACTTTTCCCCGCAAGGTCCGCCTCGTGCCCATCACCTTTCAGGCCCTGTTCGCGCCCAGCAGCCTCGCCCTCAAGTTTGCCATCAAGACCCTGCTCGGGGGCGGCCTGGCGTTGTGGCTGGCGATGCGCTggggg                                                                        >  4:150323/1‑149 (MQ=255)                                                                     aTCTTTTCAGCAAGACTTTTCCCCGCAAGGTCCGCCTCGTGCCCATCACCTTTCAGGCCCTGTTCGCGACCAGCAGCCTCGCCCTCAAGTTTGCCATCAAGACCCTGCTCGGGGGCGGCCTGGCGTTGTGGCTGGCGATGCGCTgg                                                                           <  2:222412/146‑1 (MQ=255)                                                                     aTCTTTTCAGCAAGACTTTTCCCCGCAAGGTCCGCCTCGTGCCCATCACCTTTCAGGCCCTGTTCGCGACCAGCAGCCTCGCCCTCAAGTTTGCCATCAAGACCCTGCTCGGGGGCGGCCTGGCGTTGTGGCTGGCGATGCGCTgg                                                                           >  1:222412/1‑146 (MQ=255)                                                                      tCTTTTCAGGAAGACTATTCCCGGCAAGGTGCTACTCGTGACAATCCCCTTTCATGCCCTGTTCGCGCCCATCAGACTCGCCCTCAAGTTTGGCATCAATACCCTGCTCGGGGGCTGCATGGCGTTGTGGCTGGCGATGCTCTgggt                                                                         <  8:286657/147‑2 (MQ=255)                                                                          ttCAGCAAGACTTTTCCCCGCAAGGTCCGCCTCGTGCCCATCACCTTTCAGGCCCTGTTCGCGCCCAGCAGCCTCGCCCTCAAGTTTGCCATCAAGACCCTGCTCGGGGGCGGCCTGGCGTTGTGGCTGGCGATGCGCTGGGGGCTgg                                                                    <  5:243189/148‑1 (MQ=255)                                                                              gCAAGACTTTTCCCCGCAAGGTCCGCCTCGTGCCCATCACCTTTCAGGCCCTGTTCGCGCCCAGCAGCCTCGCCCTCAAGTTTGCCATCAAGACCCTGCTCGGGGGCGGCCTGGCGTTGTGGCTGGCGATGCGCTGGGGGCTGGAGCaa                                                               <  2:219444/149‑1 (MQ=255)                                                                                aaGACTTTTCCCCGCAAGGTCCGCCTCGTGCCCATCACCTTTCAGGCCCTGTTCGCGCCCAGCAGCCTCGCCCTCAAGTTTGCCATCAAGACCCTGCTCTGGGGCGGCCTGGCGTTGTGGCTGGCGATGCGCTGGGGGCTGGAGCAAcc                                                             >  2:51341/1‑149 (MQ=255)                                                                                aaGACTTTTCCCCGCAAGGTCCGCCTCGTGCCCATCACCTTTCAGGCCCTGTTCGCGCCCAGCAGCCTCGCCCTCAAGTTTGCCATCAAGACCCTGCTCGGGGGCGGCCTGGCGTTGTGGCTGGCGATGCGCTGGGGGCTGGAGCAAcc                                                             <  3:269461/149‑1 (MQ=255)                                                                                 aGACTTTTCCCCGCAAGGTCCGCCTCGTGCCCATCACCTTTCAGGCCCTGTTCGCGCCCAGCTGCCTCGCCCTCAAGTTTGCCATCAAGACCCTGCTCGGGGGCGGCCTGGCGTTGTGGCTGGCGATGCGCTGGGGGCTGGAGCAAccc                                                            <  1:15664/149‑1 (MQ=255)                                                                                 aGACTTTTCCCCGCAAGGTCCGCCTCGTGCCCATCACCTTTCAGGCCCTGTTCGCGCCCAGCTGCCTCGCCCTCAAGTTTGCCATCAAGACCCTGCTCGGGGGCGGCCTGGCGTTGTGGCTGGCGATGCGCTGGGGGCTGGAGCAAccc                                                            <  1:9618/149‑1 (MQ=255)                                                                                 aGACTTTTCCCCGCAAGGTCCGCCTCGTGCCCATCACCTTTCAGGCCCTGTTCGCGCCCAGCAGCCTCGCCCTCAAGTTTGCCATCAAGACCCTGCTCGGGGGCGGCCTGGCGTTGTGGCTGGCGATGCGCTGCGGGCTGGAGCAAcc                                                             >  4:235505/1‑148 (MQ=255)                                                                                    cTTTTCCCCGCAAGGTCCGCCTCGTGCCCATCACCTTTCAGGCCCTGTTCGCGCCCAGCAGCCTCGCCCTCAAGTTTGCCATCAAGACCCTGCTCGGGGGCGGCCTGGCGTTGTGGCTGGCGCTGCGCAGGGGGATGGAGCAACcctcc                                                         >  4:305662/1‑149 (MQ=255)                                                                                     ttttCCCCGCAAGGTCCGCCTCGTGCCCATCACCTTTCAGGCCCTGTTCGCGCCCAGCAGCCTCGCCCTCAAGTTTGCCATCAAGAcc                                                                                                                     >  7:27897/1‑88 (MQ=255)                                                                                     ttttCCCCGCAAGGTCCGCCTCGTGCCCATCACCTTTCAGGCCCTGTTCGCGCCCAGCAGCCTCGCCCTCAAGTTTGCCATCAAGAcc                                                                                                                     <  8:27897/88‑1 (MQ=255)                                                                                       ttCCCCGCAAGGTCCGCCTCGTGCCCATCACCTTTCAGGCCCTGTTCGCGCCCAGCAGCCTCGCCCTCAAGTTTGCCATCAAGACCCTGCTCGGGGGCGGCCTGGCGTTGTGGCTGGCGATGCGCTGGGGGCTGGAGCAACCCTCCTgg                                                      >  4:356243/1‑149 (MQ=255)                                                                                         ccccGCAAGGTCCGCCTCGTGCCCATCACCTTTCAGGCTCTGTTCGCGCCCAGCAGCCTCGCCCTCAAGTTTGCCATCAAGACCCTGCTCGGGGGCGGCCTGGCGTTGTGGCTGGCGATGCGCTGGGGGCTGGAGCAACCCTCCTggg                                                     <  5:236147/148‑1 (MQ=255)                                                                                           ccGCAAGGTCCGCCTCGTGCCCATCACCTTTCAGGCCCTGTTCGCGCCCAGCAGCCTCGCCCTCAAGTTTGCCATCAAGACCCTGCTCGGGGGCGGCCTGGCGTTGTGGCTGGCGATGCGCTGGGGGCTGGAGCAACCCTCCTGGGCAt                                                  <  7:325217/149‑1 (MQ=255)                                                                                           ccGCAAGGTCCGCCTCGTGCCCATCACCTTTAAGGCCCTGTTCGTGCCCAGCAGCCTCGCCCTCAAGTTTGCCATCAAGACCCTGCTCGGGGGCGGCCTGGCGTTGTGGCTGGCGATGCGCTGGGGCCTGGAGCAACCCTCCTGGGCAt                                                  >  7:136386/1‑149 (MQ=255)                                                                                                   tCCGCCTCGTGCCCATCACCTTTCAGGCCCTGTTCGCGCCCAGCAGCCTCGCCCTCAAGTTTGCCATCAAGACCCTGCTCGGGGGCGGCCTGGCGTTGTGGCTGGCGATGCGCTGGGGGCTGGAGCAACCCTCCTGGGCATTGATGACt                                          <  5:330511/149‑1 (MQ=255)                                                                                                     cGCCTCGTGCCCATCACCTTTCAGGCCCTGTTCGCGCCCAGCAGCCTCGCCCTCAAGTTTGCCATCAAGACCCTGCTCGGGGGCGGCCTGGCGTTGTGGCTGGCGATGCGCTGGGGGCTGGAGCAACCCTCCTGGGCATTGACGACTGc                                        >  2:33184/1‑149 (MQ=255)                                                                                                        cTCGTGCCCATCACCTTTCAGGCCCTGTTCGCGCCCAGCAGCCTCGCCCTCAAGTTTGCCATCAAGACCCTGCTCGGGGGCGGCCTGGCGTTGTGGCTGGCGATGCGCTGGGGGCTGGAGCAACCCTCCTGGGCATTGATGACTGCttt                                     <  3:239061/149‑1 (MQ=255)                                                                                                        cTCGTGCCCATCACCTTTCAGGCCCTGTTCGCGCCCAGCAGCCTCGCCCTCAAGTTTGCCATCAAGACCCTGCTCGGGGGCGGCCTGGCGTTGTGGCTGGCGATGCGCTGGGGGCTGGAGCAACCCTCCTGGGCATTGATGACTGCttt                                     >  6:323519/1‑149 (MQ=255)                                                                                                           gTGCCCATCACCTTTCAGGCCCTGTTCGCGCCCAGCAGCCTCGCCCTCAAGTTTGCCATCAAGACCCTGCTCGGGGGCGGCCTGGCGTTGtggctgg                                                                                      >  5:321196/1‑97 (MQ=255)                                                                                                           gTGCCCATCACCTTTCAGGCCCTGTTCGCGCCCAGCAGCCTCGCCCTCAAGTTTGCCATCAAGACCCTGCTCGGGGGCGGCCTGGCGTTGtggctgg                                                                                      <  6:321196/97‑1 (MQ=255)                                                                                                              cccATCACCTTTCAGGCCCTGTTCGCGCCCAGCAGCCTCGCCCTCAAGTTTGCCATCAAGACCCTGCTCGGGGGCGGCCTGGCGTTGTGGCTGGCGATGCGCTGGGGGCTGGAGCAAc                                                              >  4:269830/1‑118 (MQ=255)                                                                                                              cccATCACCTTTCAGGCCCTGTTCGCGCCCAGCAGCCTCGCCCTCAAGTTTGCCATCAAGACCCTGCTCGGGGGCGGCCTGGCGTTGTGGCTGGCGATGCGCTGGGGGCTGGAGCAAc                                                              <  3:269830/118‑1 (MQ=255)                                                                                                                 aTCACCTTTCAGGCCCTGTTCGCGCCCAGCAGCCTCGCCCTCAAGTTTGCCATCAAGACCCTGCTCGGGGGCGGCCTGGCGTTGTGGCTGGCGATGCGCTGGGGGCTGGAGCAACCCTCCTGGGCATTGATGACTGCTTTCATCGTCGc                            <  2:45722/149‑1 (MQ=255)                                                                                                                    aCCTTTCAGGCCCTGTTCGCGCCCAGCAGCCTCGCCCTCAAGTTTGCCATCAAGACCCTGCTCGGGGGCGGCCTGGCGTTGTGGCTGGCGATGCGCTGGGGGCTGGAGCAACCCTCCTGGGCATTGATGACTGCTTTCATCGTCGcc                           >  2:226256/1‑147 (MQ=255)                                                                                                                     ccTTTCAGGCCCTGTTCGCGCCCAGCAGCCTCGCCCTCAAGTTTGCCATCAAGACCCTGCTCGGGGGCGGCCTGGCGTTGTGGCTGGCGATGCGCTGGGGGCTGGAGCAACCCTCCTGGGCATTGATGACTGCTTTCATCGTCGCCCAg                        >  7:220033/1‑149 (MQ=255)                                                                                                                     ccTTTCAGGCCCTGTTCGCGCCCAGCAGCCTCGCCCTCAAGTTTGCCATCAAGACCCTGCTCGGGGGCGGCCTGGCGTTGTGGCTGGCGATGCGCTGGGGGCTGGAGCAACCCTCCTGGGCATTGATGACTGCTTTCATCGTCGCCCAg                        >  3:36886/1‑149 (MQ=255)                                                                                                                      gtttGAGGGCCTGTTCGGGCCCAGCAGCCTCGCGCTCAAGTTTGGGATCAAGTCCTTGCTCGGGGGCGGGCTGGCGTTGTGGGTGGGGATGCGCTGGGGGCTGGAGcta                                                               <  4:332124/108‑3 (MQ=255)                                                                                                                      cTTTCAGGCCCTGTTCGCGCCCAGCAGCCTCGCCCTCAAGTTTGCCATCAAGACCCTGCTCGGGGGCGGCCTGGCGTTGTGGCTGGCGATGCGCTGGGGGCTGGAGCAACCCTCCTGGGCATTGATGACTGCTTTCATCGTCGc                            <  6:320516/144‑1 (MQ=255)                                                                                                                      cTTTCAGGCCCTGTTCGCGCCCAGCAGCCTCGCCCTCAAGTTTGCCATCAAGACCCTGCTCGGGGGCGGCCTGGCGTTGTGGCTGGCGATGCGCTGGGGGCTGGAGCAACCCTCCTGGGCATTGATGACTGCTTTCATCGTCGc                            >  5:320516/1‑144 (MQ=255)                                                                                                                      cTTTCAGGCCCTGTTCGCGCCCAGCAGCCTCGCCCTCAAGTTTGCCATCAAGACCCTGCTCGGGGGCGGCCTGGCGTTGTGGCTGGCGATGCGCTGGGGGCTGGAGCAACCCTCCTGGGCATTGATGACTGCTTTCATCGTCGCCCAGc                       >  8:235070/1‑149 (MQ=255)                                                                                                                               ccTGTTCGCGCCCAGCAGCCTCGCCCTCAAGTTTGCCATCAAGACCCTGCTCGGGGGCGGCCTGGCGTTGTGGCTGGCGATGCGCTGGGGGCTGGAGCAACCCTCCTGGGCATTGATGACTg                                         <  4:198888/122‑1 (MQ=255)                                                                                                                               ccTGTTCGCGCCCAGCAGCCTCGCCCTCAAGTTTGCCATCAAGACCCTGCTCGGGGGCGGCCTGGCGTTGTGGCTGGCGATGCGCTGGGGGCTGGAGCAACCCTCCTGGGCATTGATGACTg                                         >  3:198888/1‑122 (MQ=255)                                                                                                                               ccTGTTCGCGCCCAGCAGCCTCGCCCTCAAGTTTGCCATCAAGACCCTGCTCGGGGGCGGCCTGGCGTTGTGGCTGGCGATGCGCTGGGGGCTGGAGCAACCCTCCTGGGCATTGATGACTGCTTTCATCGTCGCCCAGCCGCTATCgg              <  3:314111/149‑1 (MQ=255)                                                                                                                                   ttCGCGCCCAGCAGCCTCGCCCTCAAGTTTGCCATCAAGACCCTGCTCGGGGGCGGCCTGGCGTTGTGGCTGGCGATGCGCTGGGGGCTGGAGCAACCCTCCTGGGCATTGATGACTGCTTTCATCGTCGCCCAGCCGCTATCGGGCAt          >  5:92776/1‑149 (MQ=255)                                                                                                                                   ttCGCGCCCAGCAGCCTCGCCCTCAAGTTTGCCATCAAGACCCTGCTCGGGGGCGGCCTGGCGTTGTGGCTGGCGATGCGCTGGGGGCTGGAGCAACCCTCCTGGGCATTGATGACTGCTTTCATCGTCGCCCAGCCGCTATCGGGCAt          >  2:58901/1‑149 (MQ=255)                                                                                                                                      gcgcCCAGCAGCCTCGCCCTCAAGTTTGCCATCAAGACCCTGCTCGGGGGCGGCCTGGCGTTGTGGCTGGCGATGCGCTGGGGGCTGGAGCAACCCTCCTGGGCATTGATGACTGCTTTCATCGTCGCCCAGCCGCTATCGGGCAtggt       >  6:353501/1‑149 (MQ=255)                                                                                                                                       cgcCCAGCAGCCTCGCCCTCAAGTTTGCCATCAAGACCCTGCTCGGGGGCGGCCTGGCGTTGTGGCTGGCGATGCGCTGGGGGCTGGAGCAACCCTCCTGGGCATTGATGACTGCTTTCATCGTCGCCCAGCCGCTATCGGGCAtggtg      <  8:208597/149‑1 (MQ=255)                                                                                                                                       cgcCCAGCAGCCTCGCCCTCAAGTTTGCCATCAAGACCCTGCTCGGGGGCGGCCTGGCGTTGTGGCTGGCGATGCGCTGGGGGCTGGAGCAACCCTCCTGGGCATTGATGACTGCTTTCATCGTCGCCCAGCCGCTATCGGGCAtggtg      <  6:92776/149‑1 (MQ=255)                                                                                                                                       cgcCCAGCAGCCTCGCCCTCAAGTTTGCCATCAAGACCCTGCTCGGGGGCGGCCTGGCGTTGTGGCTGGCGATGCGCTGGGGGCTGGAGCAACCCTCCTGGGCATTGATGACTGCTTTCATCGTCGCCCAGCCGCTATCGGGCAtggtg      <  2:169861/149‑1 (MQ=255)                                                                                                                                       cgcCCAGCAGCCTCGCCCTCAAGTTTGCCATCAAGACCCTGCTCGGGGGCGGCCTGGCGTTGTGGCTGGCGATGCGCTGGGGGCTGGAGCAACCCTCCTGGGCATTGATGACTGCTTTCATCGTCGCCCAGCCGCTATCGGGCAtggtg      >  1:169861/1‑149 (MQ=255)                                                                                                                                        gcCCAGCAGCCTCGCCCTCAAGTTTGCCATCAAGACCCTGCTCGGGGGCGGCCTGGCGTTGTGGCTGGCGATGCGCTGGGGGCTGGAGCAACCCTCCTGGGCATTGATGACTGCTTTCATCGTCGCCCAGCCGCTATCGGGCAtggtgg     <  3:150323/149‑1 (MQ=255)                                                                                                                                         cccAGCAGCCTCGCCCTCAAGTTTGCCATCAAGACCCTGCTCGGGGGCGGCCTGGCGTTGTGGCTGGCGATGCGCTGGGGGCTGGAGCAACCCTCCTGGGCATTGATGACTGCTTTCATCGTCGCCCAGCCGCTATCGGGCATggtggt    <  7:235070/149‑1 (MQ=255)                                                                                                                                         cccAGCAGCCTCGCCCTCAAGTTTGCCATCAAGACCCTGCTCGGGGGCGGCCTGGCGTTGTGGCTGGCGATGCGCTGGGGGCTGGAGCAACCCTCCTGGGCATTGATGACTGCTTTCATCGTCGCCCAGCCGCTATCGGGCATggtggt    <  4:36886/149‑1 (MQ=255)                                                                                                                                          ccAGCAGCCTCGCCCTCAAGTTTGCCATCAAGACCCTGCTCGGGGGCGGCCTGGCGTTGTGGCTGGCGATGCGCTGGGGGCTGGAGCAACCCTCCTGGGCATTGATGACTGCTTTCATCGTCGCCCAGCCGCTATCGGGCATGgtggtg   <  3:356243/149‑1 (MQ=255)                                                                                                                                           cagcagCCTCGCCCTCAAGTTTGCCATCAAGCACCTGCTGGGGGGCGGCCTGGGGTTGTTGCTGGCGATGCGCTGGGTGCTGGAGCAACCCTCCTGGGCATTGATGACTGCTTTCATCGTCGCCCAGCCGCTATCGGGCATGgtggtgg  <  8:136386/149‑2 (MQ=255)                                                                                                                                           cagcagCCTCGCCCTCAAGTTTGCCATCAAGACCCTGCTAGGGGGCGGCCTGGCGTTGTGGCTGGCGATGCGCTGGGGGCTGGAGCAACCCTCCTGGGCATTGATGACTGCTTTCATCGTCGCCCAGCCGCTATCGGGCATGgtggtg   >  3:41922/1‑148 (MQ=255)                                                                                                                                               |                                                                                                                                                 GCCAATCGAATATACATGATTAGCTAGCTAATTAATCGTTTCCGCTTTATCCTGCTCATCA‑TTGTCTATCTTTTCAGCAAGACTTTTCCCCGCAAGGTCCGCCTCGTGCCCATCACCTTTCAGGCCCTGTTCGCGCCCAGCAGCCTCGCCCTCAAGTTTGCCATCAAGACCCTGCTCGGGGGCGGCCTGGCGTTGTGGCTGGCGATGCGCTGGGGGCTGGAGCAACCCTCCTGGGCATTGATGACTGCTTTCATCGTCGCCCAGCCGCTATCGGGCATGGTGGTGC  >  NC_002947/1316478‑1316763

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