Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	991,569	A→T	9.2%	D63V (GAC→GTC)	der →	GTPase Der

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	991,569	0	A	T	9.2%	197.2 / 10.3	65	D63V (GAC→GTC)	der	GTPase Der
Reads supporting (aligned to +/- strand):  ref base A (29/30);  new base T (3/3);  total (32/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.63e-01

ATCCACCATGTTCAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGACCAGATGATCGCCGAGCACCTGCGCAAGCGGAACAAGGCG  >  NC_002947/991426‑991717                                                                                                                                                |                                                                                                                                                     aTCCACCATGTTCAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGTCACCGGTGATACCACCGGTGTCGatca                                                                                                                                                <  1:1287012/150‑4 (MQ=255)  tCCACCATGTTCAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTg                                                                                                                                               >  2:1158033/1‑150 (MQ=255)  tCCACCATGTTCAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTg                                                                                                                                               <  1:625067/150‑1 (MQ=255)   ccaccaTGTTCAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTgg                                                                                                                                              <  1:1117395/150‑1 (MQ=255)   ccaccaTGTTCAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTgg                                                                                                                                              <  2:1184862/150‑1 (MQ=255)      ccaTGTTCAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCa                                                                                                                                           >  2:950832/1‑150 (MQ=255)        aTGTTCAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATg                                                                                                                                         >  1:1313125/1‑150 (MQ=255)          gTTCAACCGCCTGACCAATACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGa                                                                                                                                       <  1:262827/150‑1 (MQ=255)             cAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGa                                                                                                                                    >  2:238037/1‑150 (MQ=255)               aCCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGaaa                                                                                                                                  <  1:673466/150‑1 (MQ=255)                 cGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGTCACCGGTGATACCACCGGTGTCGat                                                                                                                                                  >  1:769006/1‑131 (MQ=255)                 cGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGTCACCGGTGATACCACCGGTGTCGat                                                                                                                                                  <  2:769006/132‑2 (MQ=255)                  gCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGaaaaaa                                                                                                                               >  2:981645/1‑150 (MQ=255)                         aaaGATCCGGGATGCAATCGTCGGTGACCTGTCGGGGCTGACCCGTGACCGCAAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGAAAAAAAAAAAAAATGGCCGa                                                                                                                        <  2:440302/149‑1 (MQ=255)                          aaGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAg                                                                                                                       <  2:1313125/150‑1 (MQ=255)                               ccGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGt                                                                                                                                                <  2:250147/120‑1 (MQ=255)                               ccGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGt                                                                                                                                                >  1:250147/1‑120 (MQ=255)                               ccGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTc                                                                                                                  <  2:656000/150‑1 (MQ=255)                                    aTGCCTTCGTCGGTGTCCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTGCCTTCATCCTGATCGACGCCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCa                                                                                                             <  1:1229405/150‑1 (MQ=255)                                        cATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGc                                                                                                         >  2:964088/1‑150 (MQ=255)                                                gTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTgaag                                                                                                 >  1:654232/1‑150 (MQ=255)                                                  gACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTgaagaa                                                                                               <  2:654232/150‑1 (MQ=255)                                                     cTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTACTTTATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGcc                                                                                            >  1:532854/1‑150 (MQ=255)                                                           ggCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTAt                                                                                      >  2:740018/1‑150 (MQ=255)                                                             ccTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTAtgt                                                                                    <  1:950832/150‑1 (MQ=255)                                                              cTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTAtgtg                                                                                   <  2:1353203/150‑1 (MQ=255)                                                                   ccGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGtt                                                                              >  1:1095089/1‑150 (MQ=255)                                                                   ccGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGtt                                                                              <  1:740018/150‑1 (MQ=255)                                                                     gTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGaa                                                                                                                                                  >  1:838444/1‑80 (MQ=255)                                                                     gTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGaa                                                                                                                                                  <  2:838444/80‑1 (MQ=255)                                                                          cGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGtt                                                                       <  1:238037/150‑1 (MQ=255)                                                                              aGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATg                                                                   >  1:368860/1‑150 (MQ=255)                                                                               gTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGc                                                                  <  2:368860/150‑1 (MQ=255)                                                                                       aTGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCgg                                                          <  1:45382/150‑1 (MQ=255)                                                                                       aTGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCgg                                                          >  2:140800/1‑150 (MQ=255)                                                                                            aGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGt                                                                                                                                                >  2:740601/1‑59 (MQ=255)                                                                                            aGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGt                                                                                                                                                <  1:740601/59‑1 (MQ=255)                                                                                            aGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGa                                                                                          <  1:903977/113‑1 (MQ=255)                                                                                            aGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGa                                                                                          >  2:903977/1‑113 (MQ=255)                                                                                              cTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTAt                                                      <  2:1203799/147‑1 (MQ=255)                                                                                              cTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTAt                                                      >  1:1203799/1‑147 (MQ=255)                                                                                                  cAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCg                                               <  1:964088/150‑1 (MQ=255)                                                                                                   aGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCTTGGCTCATGGCCGTGCGCTTATTACtcctg                                              >  2:307111/1‑146 (MQ=255)                                                                                                       tCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGa                                          <  2:1101258/150‑1 (MQ=255)                                                                                                         gTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGAcc                                        >  1:951883/1‑150 (MQ=255)                                                                                                          ttCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTg                                                                          >  2:1217152/1‑115 (MQ=255)                                                                                                          ttCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTg                                                                          <  1:1217152/115‑1 (MQ=255)                                                                                                          ttCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGACCa                                       >  1:688779/1‑150 (MQ=255)                                                                                                                 aTCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGAt                                                                    >  2:170173/1‑114 (MQ=255)                                                                                                                 aTCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGAt                                                                    <  1:170173/114‑1 (MQ=255)                                                                                                                  tCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGACCAGATGATCg                               >  2:1032706/1‑150 (MQ=255)                                                                                                                     tGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGGCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGACCAGATGATGGCCg                            >  2:1319267/1‑150 (MQ=255)                                                                                                                     acttCGTCACCGGTGATACCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGACCAGATGATCGCCg                            >  1:255986/4‑150 (MQ=255)                                                                                                                       ttcgtcacCGGTGATACCACCGGTGTCGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCgg                                                          >  1:483405/6‑118 (MQ=39)                                                                                                                       ttcgtcacCGGTGATACCACCGGTGTCGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCgg                                                          <  2:483405/113‑1 (MQ=39)                                                                                                                       aTCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGACCAGATGATCGCCGAg                          >  2:1304910/1‑150 (MQ=255)                                                                                                                          gtcacCGGTGATACCACCGGTGTCGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGACCAGATGATCGCCGAGCAc                       >  2:129692/3‑150 (MQ=255)                                                                                                                            cacCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGACCAGATGATCGCCGAGCACCt                     <  1:1032706/150‑1 (MQ=255)                                                                                                                            cacCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGACCAGATGATCGCCGAGCACCt                     <  1:981645/150‑1 (MQ=255)                                                                                                                                gatggtATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGACCAGATGATCGCCGAGCACCTgcgc                 >  2:483560/3‑150 (MQ=255)                                                                                                                                 gtgataccACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGACCAGATGATCGCCGAGCACCTGCGCa                >  2:405744/8‑150 (MQ=255)                                                                                                                                    gtATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGACCAGATGATCGCCGAGCACCTGCGCAAGc             >  2:1277529/1‑150 (MQ=255)                                                                                                                                      aTCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGACCAGATGATCGCCGAGCACCTGCGCAAGCgg           <  2:951883/150‑1 (MQ=255)                                                                                                                                         aCCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGACCAGATGATCGCCGAGCACCTGCGCAAGCGGAAc        <  1:405744/150‑1 (MQ=255)                                                                                                                                               gACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCGGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGACCAGATGATCGCCGAGCACCTGCGCAAGCGGAACAAggcg  <  1:1158033/150‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                     ATCCACCATGTTCAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGACCAGATGATCGCCGAGCACCTGCGCAAGCGGAACAAGGCG  >  NC_002947/991426‑991717

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