Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,388,041	T→G	T50P (ACG→CCG)	aaeB ←	aromatic carboxylic acid efflux pump subunit AaeB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,388,041	0	T	G	97.8%	148.0 / NA	45	T50P (ACG→CCG)	aaeB	aromatic carboxylic acid efflux pump subunit AaeB
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (23/21);  total (23/22)
Fisher's exact test for biased strand distribution p-value = 4.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACAATGATCAGCGCGGTATAACCGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGT‑ACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGC  >  NC_000913/3387794‑3388286                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       acaATGACCAGCGCAGTGTAACCGGCCAGCCCCCACGCATACGAGCTTTCTCTTCCTACCAGCGAGGATCTCCAGGTACACAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGCGCGCGGATCATCGCAATGATGATCACCAGTCCGGCACTACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGACTCCGGTTCACCTCCCGCAGCAACGGCCGG‑Acc                                                                                                                                                                                                                                                    >  2:65489/1‑251 (MQ=255)  caATGATCAGCGCGGTATAACCGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGG‑ACcg                                                                                                                                                                                                                                                   >  1:419916/1‑251 (MQ=255)     tGATCAGCGCGGTATAACCGGCCAGCCCCCACGCATACGAGATTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAAAGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCACCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGG‑ACcggcg                                                                                                                                                                                                                                                >  1:267277/1‑251 (MQ=255)                              cccccACGCATACGAGGTTTCTAGTCGTACCAGCGAGGATATCCAGGTACAAAATCCTGCCCAGCTACAGCACCCCAGAATCGTCAATAATGGGGCGCGGATCATCTCAATGATGATCACCAGGCCGGCAAGACAGCCAATAAATGTGCGGATGATGCGTAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGG‑ACCGGCGGCAACAATCGCCGCTGTCAGTACCg                                                                                                                                                                                                                       <  1:480764/251‑1 (MQ=255)                               ccccACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGG‑ACCGGCGGCAACAATCGCCGCTGTCAGTACCGc                                                                                                                                                                                                                      <  1:136970/251‑1 (MQ=255)                                                ttCTATTCGTACCAGCGAGGCTGTCCAGGTACAAAAACCGGCCCAGAGACAGCCCACCAGAATCATCAATAATGGTGCGCGGTTCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGATGAGCGCCAGAAGACGGTTCACCTCCCGCAGCAAAGGCCGG‑ACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTcc                                                                                                                                                                                                      <  2:416392/250‑1 (MQ=255)                                                 tCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGG‑ACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCa                                                                                                                                                                                                     <  1:39230/250‑1 (MQ=255)                                                           ccAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGG‑ACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTg                                                                                                                                                                                           >  1:260040/1‑250 (MQ=255)                                                            cagtgagtgtatCCAGGTCCAAACCCCGGCCCAGATACAGCACCCCAGTATAATCAATAATGGTGCGCGGATCATCGCAATGATGATCCCCAGTCCTGCAATACAGCCAATAAATGTGCCGATGATGCGCAATAAGCCACTATAGCGAATAGCGCCAGAATTCGGTTCCCCTCCCGCAGCAAAGGCCGG‑ACCTGCGGCACCAATCGCCGCTGTCAGTACCGCCCAGCGTGGCTGTTCCAGCGGGAAGTgca                                                                                                                                                                                         <  2:389852/242‑3 (MQ=255)                                                              gCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGG‑ACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAg                                                                                                                                                                                       <  2:172461/251‑1 (MQ=255)                                                                             aagactttaccagccAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATCCGGTTCACCTCCCGCAGCAAAGGCCGG‑ACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGcc                                                                                                                                                                        >  2:229266/14‑251 (MQ=255)                                                                                aaaaaCCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGG‑ACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGAGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGt                                                                                                                                                                     >  1:100719/1‑251 (MQ=255)                                                                                   aaCCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGG‑ACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCg                                                                                                                                                                                                           <  2:392459/210‑1 (MQ=255)                                                                                   aaCCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGG‑ACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCg                                                                                                                                                                                                           >  1:392459/1‑210 (MQ=255)                                                                                           cAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACTATTCCGGCAATACAGCCAATAAATGTGCCTATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATGCGGTTCACCTCCCGCAGCAAAGGCCGG‑ACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCACCCTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGt                                                                                                                                                          <  1:19028/251‑1 (MQ=255)                                                                                           cAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGG‑ACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGt                                                                                                                                                          <  2:9696/251‑1 (MQ=255)                                                                                                        ccAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGG‑ACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTAcc                                                                                                                                             >  1:12562/1‑251 (MQ=255)                                                                                                               catcaATAATGGTGCGCGGATCGTCGCAATGATGATCACCCGGCCGGCAATACAGCCAATAAATGTGCCGATGATCCACAAAAAGCCCCGCCAGCGAAGAGCGCCAGAATACGGTTCACCTCACCCAGCAAAGGCCGG‑ACCGGCGGCAACAATCGCCGCTGTCAGTGCCGCCCAGCGTGGCGTCTCCAGCTGCAAGCGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCaaat                                                                                                                                       <  1:5123/250‑2 (MQ=255)                                                                                                               catcaATAATGGTGCGCGGAGCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGG‑ACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAg                                                                                                                                       <  1:438119/250‑1 (MQ=255)                                                                                                                                ggATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGG‑ACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCa                                                                                                                     <  1:310909/251‑1 (MQ=255)                                                                                                                                       cGCAATGATGATCACCAGTCCGTCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGG‑ACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGgag                                                                                                               <  2:462745/250‑1 (MQ=255)                                                                                                                                            tgatgaTCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGG‑ACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAAt                                                                                                          >  1:168171/1‑250 (MQ=255)                                                                                                                                                  tCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGG‑ACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTAGCCAGTTTTACCGCCAAGAGAATATGGTGGTTTGCAATGGAGAAAATACCCATc                                                                                                   >  1:133491/1‑251 (MQ=255)                                                                                                                                                      cAGTACGGCAATACAGCCAATAAAGGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATGCGGTGCACCTCCCGCAGCAAAGGCCGG‑ACCGGCGGCAACAATCGCCGCTGTCAGTACCGCGCAGGGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGGCAGTACAATGGCGGTCGCGAGTTTTAGCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGa                                                                                               <  1:65489/251‑1 (MQ=255)                                                                                                                                                          ccGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGG‑ACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTaa                                                                                           <  2:260040/251‑1 (MQ=255)                                                                                                                                                           cGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGG‑ACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGGGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGAACAATGGCGGTCGGCAGGTTTACCGCCAAACGGAAAAGATGGGTAGCAAAGGAGAAAAAAACCATCgcgaataa                                                                                           >  2:339959/1‑243 (MQ=255)                                                                                                                                                                   cAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGG‑ACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCCGTTTTTCCGCAAAGCGAATATGTTGGTTAGCAATTGAGAAAATACCCCTCGTGAtttaacaaaatcc                                                                                  >  2:382198/1‑240 (MQ=255)                                                                                                                                                                    aGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGG‑ACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTTCCGCAAAGCGCATATGTTGGGTAGCAATGGAGAAAAAAACCATCGTGATTAACCAAACttc                                                                                  >  1:346217/1‑248 (MQ=255)                                                                                                                                                                         aTAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGG‑ACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGGTTTACCGCAAAGCGAATATGTTGGTTAGGAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAgg                                                                            >  1:147333/1‑251 (MQ=255)                                                                                                                                                                                 gCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGG‑ACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGc                                                                    <  1:229266/251‑1 (MQ=255)                                                                                                                                                                                    gatgatGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGG‑ACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGACGCCAGTTTTACCCCAAAACGAATATGTTGGGTAGCAAAGGAGGAAAAACCCATACTGATTAACCCAACTCCCCCCGGCGAAGGGccct                                                                 >  2:397913/1‑249 (MQ=255)                                                                                                                                                                                    gatgatGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCAGCAGCAAAGGCCGG‑ACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTTCCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAAATCACGCAGGCGATGGGCCAt                                                                 >  1:48400/1‑251 (MQ=255)                                                                                                                                                                                         tgCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGG‑ACCGGTGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCCCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCCTCGTGATTAACCCAACCCACGCAGGCGATGGGCCCTTTTAc                                                            >  1:330035/1‑251 (MQ=255)                                                                                                                                                                                                 aaGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGG‑ACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAAc                                                    <  1:474658/251‑1 (MQ=255)                                                                                                                                                                                                    ccACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGG‑ACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGa                                                  >  1:50764/1‑250 (MQ=255)                                                                                                                                                                                                     cACGATAGCGAATAGCGCCAGCATACGGTTCACCTCCCGCAGCAAAGGCCGG‑ACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTAAAATGGCGGGTGCCAGGTTTACCGCAAAGCGAATATGCTGGGTAGCAATGGAGAAAATACCCATCGTGGTTAACCAAACTCCAGCAGGCGAGGGGCCATTTTACGGGAGAACGAAt                                                >  2:20588/1‑251 (MQ=255)                                                                                                                                                                                                           aGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGG‑ACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTgg                                           <  1:188347/250‑1 (MQ=255)                                                                                                                                                                                                            gCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGG‑ACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCt                                         <  2:12562/251‑1 (MQ=255)                                                                                                                                                                                                                  aGCGCCAGAAGACGGTTCACCTCCCGCAGCAAAGGCCGG‑ACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTtcgtc                                   <  2:419916/251‑1 (MQ=255)                                                                                                                                                                                                                        aGAATACGCGTCACCGCCCGCAGCAAAGGC‑GTGACCGTCGGCAACAATGGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGTCGGTCTCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAAGCCTCGCTTTCGTCTCGATc                             <  2:61492/251‑1 (MQ=255)                                                                                                                                                                                                                                          cGCAGCAAAGGCCGG‑ACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCCGTTTTACCGCCAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGaccac           >  1:8887/1‑251 (MQ=255)                                                                                                                                                                                                                                          cGCAGCAAAGGCCGG‑ACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCa                                                                                  >  2:426360/1‑180 (MQ=255)                                                                                                                                                                                                                                          cGCAGCAAAGGCCGG‑ACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCa                                                                                  <  1:426360/180‑1 (MQ=255)                                                                                                                                                                                                                                           gcagcaAAGGCCGG‑ACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCAAGCCGGCGATGGGCCATTTTACGGGAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGaccacc          >  1:406058/1‑251 (MQ=255)                                                                                                                                                                                                                                                   ggCCGG‑ACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAATGACCACCACTGTAGc  >  1:359865/1‑251 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       ACAATGATCAGCGCGGTATAACCGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGT‑ACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGC  >  NC_000913/3387794‑3388286

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