Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,315,811	A→G	11.9%	F91S (TTC→TCC)	phnN ←	ribose 1,5‑bisphosphokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,315,811	0	A	G	11.9%	191.0 / 15.7	67	F91S (TTC→TCC)	phnN	ribose 1,5‑bisphosphokinase
Reads supporting (aligned to +/- strand):  ref base A (31/28);  new base G (4/4);  total (35/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 = 7.40e-01

TGGCGTTTTCACGGCCACGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGATATGGTTTTCACTTCCGGCGCTGGC  >  NC_000913/4315671‑4315959                                                                                                                                             |                                                                                                                                                     tGGCGTTTTCACGGCCACGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCg                                                                                                                                              >  2:1077146/1‑149 (MQ=255)              gCCACGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGAt                                                                                                                                 >  2:1212718/1‑149 (MQ=255)                cACGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCa                                                                                                                               >  2:277280/1‑149 (MQ=255)                          cAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGCCGGATTCGACGTGc                                                                                                                                           >  2:578919/1‑127 (MQ=255)                          cAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGCCGGATTCGACGTGc                                                                                                                                           <  1:578919/127‑1 (MQ=255)                           aGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTcgacgcc                                                                                                                    >  1:980424/1‑149 (MQ=255)                            ggCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTcgacgccg                                                                                                                   >  3:673089/1‑149 (MQ=255)                             gcgcTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAAtc                                                                                                                            <  2:540918/139‑1 (MQ=255)                             gcgcTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAAtc                                                                                                                            >  1:540918/1‑139 (MQ=255)                              cgcTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTcgacgccgac                                                                                                                 <  1:1077146/149‑1 (MQ=255)                                     ggaggaTCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGc                                                                                                                                               >  3:1248218/1‑112 (MQ=255)                                     ggaggaTCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCataa                                                                                                          >  4:467790/1‑149 (MQ=255)                                     ggaggaTCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGAAAGCGCGCCCGCGCCTGCGGCAGATGGGCGAGTGAGCCGTTGACCAGCACGTCGAATCCGGc                                                                                                                                               <  4:1248218/112‑1 (MQ=255)                                          aTCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCgg                                                                                                                                                <  1:245608/106‑1 (MQ=255)                                          aTCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCgg                                                                                                                                                >  2:245608/1‑106 (MQ=255)                                                 gCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTa                                                                                              <  4:673089/149‑1 (MQ=255)                                                   gAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGc                                                                                            <  1:1212718/149‑1 (MQ=255)                                                         tGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAAt                                                                                                                                                    >  1:638679/1‑87 (MQ=255)                                                         tGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAAt                                                                                                                                                    >  3:1127260/1‑87 (MQ=255)                                                         tGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAAt                                                                                                                                                    <  4:1127260/87‑1 (MQ=255)                                                         tGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAAt                                                                                                                                                    <  2:638679/87‑1 (MQ=255)                                                         tGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCa                                                                                      >  4:1418304/1‑149 (MQ=255)                                                         tGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCa                                                                                      >  3:1054567/1‑149 (MQ=255)                                                         tGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGCCGGATTCGACGTGc                                                                                                                                           >  4:413351/1‑96 (MQ=255)                                                         tGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGCCGGATTCGACGTGc                                                                                                                                           <  3:413351/96‑1 (MQ=255)                                                          gTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCg                                                                                           >  4:654385/1‑143 (MQ=255)                                                          gTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCg                                                                                           <  3:654385/143‑1 (MQ=255)                                                            aaaCAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTcgacg                                                                                                                      <  4:409601/114‑1 (MQ=255)                                                            aaaCAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTcgacg                                                                                                                      >  3:409601/1‑114 (MQ=255)                                                                aGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAg                                                                               <  3:1418304/149‑1 (MQ=255)                                                                          cagcGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCgacgccgacgcc                                                                                                              >  2:1072429/1‑108 (MQ=255)                                                                          cagcGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCgacgccgacgcc                                                                                                              <  1:1072429/108‑1 (MQ=255)                                                                                cGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGcc                                                               <  1:277280/149‑1 (MQ=255)                                                                                  aTTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGcccc                                                             >  3:216823/1‑149 (MQ=255)                                                                                      aTAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCa                                                                                                             >  1:1178394/1‑97 (MQ=255)                                                                                      aTAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCa                                                                                                             <  2:1178394/97‑1 (MQ=255)                                                                                           gcgcccgcgccTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGAttt                                                                   >  3:890019/1‑134 (MQ=255)                                                                                           gcgcccgcgccTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGAttt                                                                   <  4:890019/134‑1 (MQ=255)                                                                                               ccgcgccTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTaaaaa                                                >  4:511406/1‑149 (MQ=255)                                                                                                  cgccTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACt                                             <  3:511406/149‑1 (MQ=255)                                                                                                          gCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTgctcg                                     >  1:1209593/1‑149 (MQ=255)                                                                                                            aGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTc                                            <  3:172956/140‑1 (MQ=255)                                                                                                            aGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTc                                            >  4:172956/1‑140 (MQ=255)                                                                                                             gATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTgctcgctc                                  >  3:1415503/1‑149 (MQ=255)                                                                                                                  gcgcGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGc                             <  3:467790/149‑1 (MQ=255)                                                                                                                      gTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGc                                                                                    <  2:564000/90‑1 (MQ=255)                                                                                                                      gTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGc                                                                                    >  1:564000/1‑90 (MQ=255)                                                                                                                         aGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGATATgg                      >  4:1314646/1‑149 (MQ=255)                                                                                                                          gCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGATATGGt                     <  4:216823/149‑1 (MQ=255)                                                                                                                             gTTGACCGGCACGTCGACTGCGGCGTGCAGCCAGACACCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAACGGTAACAGATTTTGCCCCGCGCGGGAAAAAAACTCCTGCTCGCTCAGGGCGATATGGttgt                  <  4:1054567/149‑3 (MQ=255)                                                                                                                                   cAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCg                                                                                           <  4:652868/70‑1 (MQ=255)                                                                                                                                   cAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCg                                                                                           >  3:652868/1‑70 (MQ=255)                                                                                                                                    tgCACGCCGGATTCGACGTGCTGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTgc                                        <  1:872597/119‑1 (MQ=255)                                                                                                                                    tgCACGCCGGATTCGACGTGCTGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTgc                                        >  4:1418454/2‑120 (MQ=255)                                                                                                                                    tgCACGCCGGATTCGACGTGCTGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTgc                                        >  2:872597/2‑120 (MQ=255)                                                                                                                                    tgCACGCCGGATTCGACGTGCTGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTgc                                        <  3:1418454/119‑1 (MQ=255)                                                                                                                                    aGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGATATGGTTTTCACTTcc           <  3:1314646/149‑1 (MQ=255)                                                                                                                                      cACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGc                                                                >  4:1147134/1‑94 (MQ=255)                                                                                                                                      cACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGc                                                                <  3:1147134/94‑1 (MQ=255)                                                                                                                                         gTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTaaaaa                                                >  1:814069/1‑107 (MQ=255)                                                                                                                                         gTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTaaaaa                                                <  2:814069/107‑1 (MQ=255)                                                                                                                                          tCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTaaaa                                                 <  2:1325618/105‑1 (MQ=255)                                                                                                                                          tCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTaaaa                                                 >  1:1325618/1‑105 (MQ=255)                                                                                                                                          tCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGATATGGTTTTCACTTCCGGCGCt     <  4:1415503/149‑1 (MQ=255)                                                                                                                                           cGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGATATGGTTTTCACTTCCGGCGCTg    >  3:483407/1‑149 (MQ=255)                                                                                                                                            gAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAgg                               <  2:1091699/121‑1 (MQ=255)                                                                                                                                            gAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAgg                               >  1:1091699/1‑121 (MQ=255)                                                                                                                                             aaTCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGATATGGTTTTCa                <  2:127163/135‑1 (MQ=255)                                                                                                                                             aaTCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGATATGGTTTTCa                >  1:127163/1‑135 (MQ=255)                                                                                                                                             aaTCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGATATGGTTTTCACTTCCGGCGCTGGc  <  2:1209593/149‑1 (MQ=255)                                                                                                                                             |                                                                                                                                                     TGGCGTTTTCACGGCCACGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGATATGGTTTTCACTTCCGGCGCTGGC  >  NC_000913/4315671‑4315959

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