Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	953,720	A→T	5.9%	N249K (AAT→AAA)	focA ←	formate channel

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	953,720	0	A	T	5.9%	235.1 / 5.2	68	N249K (AAT→AAA)	focA	formate channel
Reads supporting (aligned to +/- strand):  ref base A (34/30);  new base T (2/2);  total (36/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.61e-01

AAGTGGATTTTTTATTTACTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATC  >  NC_000913/953572‑953867                                                                                                                                                     |                                                                                                                                                    aaGTGGATTTTTTATTTACTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGaaa                                                                                                                                                     <  3:194959/149‑1 (MQ=255)    tGGATTTTTTATTTACTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATtc                                                                                                                                                  >  2:90203/1‑149 (MQ=255)       aTTTTTTATTTACTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATtcatc                                                                                                                                               >  3:320895/1‑149 (MQ=255)         tttttATTTACTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCAc                                                                                                                                             >  2:462389/1‑149 (MQ=255)          ttttATTTACTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACg                                                                                                                                            >  3:96506/1‑149 (MQ=255)           tttatttaCTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACgg                                                                                                                                           <  2:518407/149‑1 (MQ=255)                   cTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTga                                                                                                                                   <  1:76970/149‑1 (MQ=255)                   cTGCGTACTTCGACAACCATTAATCGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTga                                                                                                                                   <  2:406875/149‑1 (MQ=255)                    tGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTgag                                                                                                                                  >  1:600714/1‑149 (MQ=255)                      cGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGaa                                                                                                                                <  2:343852/149‑1 (MQ=255)                       gTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGaaa                                                                                                                               <  2:600714/149‑1 (MQ=255)                         aCTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGaaaaa                                                                                                                             <  2:504023/149‑1 (MQ=255)                         aCTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGaaaaa                                                                                                                             <  1:348508/149‑1 (MQ=255)                             cGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTCAGAAAAAtttt                                                                                                                         <  1:462389/149‑1 (MQ=255)                              gACAACCATCAATGATGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTc                                                                                                                        <  4:257046/149‑1 (MQ=255)                               aCAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTcc                                                                                                                       <  1:414945/149‑1 (MQ=255)                               aCAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTcc                                                                                                                       >  4:386533/1‑149 (MQ=255)                               aCAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTAGCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTcc                                                                                                                       >  4:506726/1‑149 (MQ=255)                                 aaCCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCgg                                                                                                                     >  2:13869/1‑149 (MQ=255)                                   ccATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTg                                                                                                                   >  4:80620/1‑149 (MQ=255)                                               gtCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTg                                                                                                       >  2:304294/1‑149 (MQ=255)                                                 cGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCg                                                                                                     >  4:165094/1‑149 (MQ=255)                                                  gTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACGGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCATTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCgg                                                                                                    >  4:287626/1‑149 (MQ=255)                                                  gTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCgg                                                                                                    <  3:203552/149‑1 (MQ=255)                                                               tAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCgg                                                                                       <  4:52133/149‑1 (MQ=255)                                                               tAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCgg                                                                                       >  2:285638/1‑149 (MQ=255)                                                                 aaTGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCgggg                                                                                     <  1:13869/149‑1 (MQ=255)                                                                   tGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGAt                                                                                   <  3:88907/149‑1 (MQ=255)                                                                     aCCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGc                                                                                 >  3:183839/1‑149 (MQ=255)                                                                      cccAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACAGACTGCGGTCCAAAATTCCGGGGATGCg                                                                                >  2:403929/1‑149 (MQ=255)                                                                        cAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGaa                                                                              <  3:80620/149‑1 (MQ=255)                                                                        cAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGaa                                                                              >  4:162972/1‑149 (MQ=255)                                                                              gTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGAGCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTcgcg                                                                        >  4:359030/1‑149 (MQ=255)                                                                               tCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTga                                                                                                                                   >  2:393042/1‑89 (MQ=255)                                                                               tCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTga                                                                                                                                   <  1:393042/89‑1 (MQ=255)                                                                               tCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCgg                                                                       >  4:250507/1‑149 (MQ=255)                                                                                  acccaaccAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGAtt                                                                    >  2:260330/1‑149 (MQ=255)                                                                                                ccaccaccGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCgg                                                      >  3:308819/1‑149 (MQ=255)                                                                                                 caccaccGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCggg                                                     >  1:225281/1‑149 (MQ=255)                                                                                                    caccGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCg                                                                                        <  3:100837/111‑1 (MQ=255)                                                                                                    caccGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCg                                                                                        >  4:100837/1‑111 (MQ=255)                                                                                                     accGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCa                                                 >  1:442931/1‑149 (MQ=255)                                                                                                     accGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCa                                                 >  1:39051/1‑149 (MQ=255)                                                                                                         ttAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCAt                                                >  4:512066/2‑146 (MQ=255)                                                                                                         aTAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATaaa                                             <  2:177448/149‑1 (MQ=255)                                                                                                         aTAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATaaa                                             <  3:386533/149‑1 (MQ=255)                                                                                                         aTAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAACTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCAt                                                <  3:512066/146‑1 (MQ=255)                                                                                                           aaTGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACa                                           <  4:96506/149‑1 (MQ=255)                                                                                                            aTGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACAt                                          <  4:320895/149‑1 (MQ=255)                                                                                                                tGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGttt                                      <  1:203669/149‑1 (MQ=255)                                                                                                                tGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGttt                                      <  1:462519/149‑1 (MQ=255)                                                                                                                 gCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTg                                     >  2:217441/1‑149 (MQ=255)                                                                                                                  ccGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGc                                    <  3:424718/149‑1 (MQ=255)                                                                                                                    gATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGa                                  >  1:131853/1‑149 (MQ=255)                                                                                                                    gATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGa                                  <  1:285638/149‑1 (MQ=255)                                                                                                                     aTCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGAt                                 >  3:130417/1‑149 (MQ=255)                                                                                                                          aaCCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACtgt                            >  1:592508/1‑149 (MQ=255)                                                                                                                           aCCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACtgtg                           >  3:242173/1‑149 (MQ=255)                                                                                                                             cGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCt                         <  2:225281/149‑1 (MQ=255)                                                                                                                                     ggTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCg                 <  3:359030/149‑1 (MQ=255)                                                                                                                                      gTTATCAGTGATGACATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCg                                                                                        <  4:393343/77‑1 (MQ=255)                                                                                                                                      gTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCg                                                                                        >  3:393343/1‑77 (MQ=255)                                                                                                                                      gTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGc                >  4:618644/1‑149 (MQ=255)                                                                                                                                         aTCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTgg             <  4:391642/149‑1 (MQ=255)                                                                                                                                            cgTGATGAATTTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACt                              >  2:391369/2‑129 (MQ=255)                                                                                                                                            cgTGATGAATTTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACt                              <  1:391369/128‑1 (MQ=255)                                                                                                                                               gatgaATTTCATCACTGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATAc                               >  4:322755/1‑125 (MQ=255)                                                                                                                                               gatgaATTTCATCACTGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATAc                               <  3:322755/125‑1 (MQ=255)                                                                                                                                                  gaAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCaacaaaca    >  3:109233/1‑149 (MQ=255)                                                                                                                                                   aaaTTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACAt   >  2:198584/1‑149 (MQ=255)                                                                                                                                                    aaTTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATc  <  4:308819/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    AAGTGGATTTTTTATTTACTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATC  >  NC_000913/953572‑953867

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