Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	472,849	A→G	G84G (GGA→GGG)	glnK →	nitrogen regulatory protein PII‑2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	472,849	0	A	G	82.6%	103.9 / 17.3	46	G84G (GGA→GGG)	glnK	nitrogen regulatory protein PII‑2
Reads supporting (aligned to +/- strand):  ref base A (3/5);  new base G (20/18);  total (23/23)
Fisher's exact test for biased strand distribution p-value = 6.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.19e-01

GCTGGTGACCGTGATAATCAAACCATTCAAGCTGGAAGACGTTCGTGAAGCGTTATCTTCCATTGGTATTCAGGGCCTGACCGTCACCGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGAAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCCGACGAAGCGGCGCTGTAATCTCTGGCACACAGCAACAGGAACGAAAAATGAAGATAGCGACGATAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTG  >  NC_000913/472603‑473091                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   gCTGGTGACCGTGATAATCAAACCATTCAAGCTGGAAGACGTTCGTGAAGCGTTAGCTTCCATTGGTATTCAGGGCCTGACCGTCACCGAAGTGAAAGGTTTCGGGCGTTAGAAAGTTCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGAGGAAGTGATCGATATCGTCAGTAAGGCGGTTTACACCGGgaaaa                                                                                                                                                                                                                                                <  2:310556/251‑1 (MQ=255)   tGGTGACCGTGATAATCAAACCATTCAAGCTGGAAGACGTTCGTGAAGCGTTATCTTCCATTGGTATTCAGGGCCTGACCGTCACCGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGGAAAAtt                                                                                                                                                                                                                                              >  1:213954/1‑251 (MQ=255)           gTGATAATCAAACCATTCAAGCTGGAAGACGTTCGTGAAGCGTTATCTTCCATTGGTATTCAGGGCCTGACCGTCACCGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCTTCAGTAATGCGGCTTACTCCGGGAAAATTGGCGACgg                                                                                                                                                                                                                                      <  1:171347/251‑1 (MQ=255)                     aaCCATTCAAGCTGGAAGACGTTCGTGAAGCGTTATCTTCCATTGGTATTCAGGGCCTGACCGTCACCGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTGCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGGAAAATTGGCGACGGCAAAATCTtc                                                                                                                                                                                                                            <  2:410418/251‑1 (MQ=255)                              aGCTGGAAGACGTTCGTGAAGCGTTATCTTCCATTGGTATTCAGGGCCTGACCGTCACCGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGGAAAATTGGCGACGGCAAAATCTTCGTCGCTGa                                                                                                                                                                                                                    >  2:20926/1‑250 (MQ=255)                                 gggAAGACGTTCGTGAAGCGTTATCTTCCAGTGGTATTCAGGGCCTGACCGTCACCGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGGAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTg                                                                                                                                                                                                                <  2:213954/250‑1 (MQ=255)                                           tCGTGAAGCGTTATCTTCCATTGGTATTCAGGGCCTGACCGTCACCGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGGAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTc                                                                                                                                                                                                       <  2:389527/250‑1 (MQ=255)                                            cGTGAAGCGTTATCTACCATTGGTATTCAGGGCCTGACCGTCACCGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGGAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATGGCAACGCGTCa                                                                                                                                                                                                      <  1:171015/250‑1 (MQ=255)                                                   cGTTATCTTCCATTGGTATTCAGGGCCTGACCGTCACCGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGGAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCattcgtatt                                                                                                                                                                                              >  1:261994/1‑251 (MQ=255)                                                            ccATTGGTATTCAGGGCCTGACCGTCACCGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGAAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGc                                                                                                                                                                                     <  1:238444/251‑1 (MQ=255)                                                                tGGTATTCAGGGCCTGACCGTCACCGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGAAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGaa                                                                                                                                                                                  <  1:1876/250‑1 (MQ=255)                                                                  ggaTTCAGGGCCTGACCGTCACCTAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAGAATTAAAAATTGATGTAGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCATCAGTAAGGCGGCTTACACCGGAAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCAATCGTATTCGTACCGGCGAAGcc                                                                                                                                                                               <  2:366083/249‑1 (MQ=255)                                                                    aTTCAGGGCCTGACCGTCACCGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGGAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCcg                                                                                                                                                                              <  2:368697/250‑1 (MQ=255)                                                                        aGGGCCTGACCGTCACCGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGAAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCcgacga                                                                                                                                                                          >  2:239635/1‑250 (MQ=255)                                                                                     cACCGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGGAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCCGACGAAGCGGCGCTGTAAt                                                                                                                                                            <  2:289743/251‑1 (MQ=255)                                                                                                  gttttCGGGCGTCAGAAAGGGCATGCCGAGCTGGACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGGAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCTTCAGTCGTATTCGTCCCGGCGAAGCCGACGAAGCGGCGCTTTAATATCTGGCACACAg                                                                                                                                               <  2:303087/249‑1 (MQ=255)                                                                                                            gTCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGAGATCGATATCGTCAGTAAGGCGGCTTACACCGGAAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCCGACGAAGCGGCGCTGTAATCTCTGGCACACAGCAACAGGAAc                                                                                                                                     >  1:44908/1‑251 (MQ=255)                                                                                                              cAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGGAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCCGACGAAGCGGCGCTGTAATCTCTGGCACACAGCAACAGGAACg                                                                                                                                    >  1:227817/1‑250 (MQ=255)                                                                                                                   aGGGCATGACGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGAAAAAATTGATGTGGCGATTGCTGATGACCAAATCGATGAAGTGATCGATATCGTCAGGAAGGCGGCTTACACCGGGAAAAATGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCAGGCGAAGCGGCGCTGTAATCTCTGGCAAACAGCAACAGGAACGAAAAAt                                                                                                                              >  2:487866/1‑251 (MQ=255)                                                                                                                            cGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGGAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCCGACGAAGCGGCGCTGTAATCTCTGGCACACAGCAACAGGAACGAAAAATGAAGATAGc                                                                                                                     >  1:423895/1‑251 (MQ=255)                                                                                                                                 tGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGGAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCCGACGAAGCGGCGCTGTAATCTCTGGCACACAGCAACAGGAACGAAAAATGAAGATAGcgacg                                                                                                                 >  1:464340/1‑250 (MQ=255)                                                                                                                                       ggggggCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGGAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCCGACGAAGCGGCGCTGTAATCTCTGGCACACAGCAACAGGAACGAAAAATGAAGATAGCGACGATcaaa                                                                                                           >  1:74036/1‑250 (MQ=255)                                                                                                                                        gggggCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGAGTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGAAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCCGACGAAGCGGCGCTGTAATCTCTGGCACACAGCAACAGGAACGAAAAATGAAGATAGCGACGATAAAAAc                                                                                                         <  1:214603/251‑1 (MQ=255)                                                                                                                                            gCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGGAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCCGACGAAGCGGCGCTGTAATCTCTGGCACACAGCAACAGGAACGAAAAATGAAGATAGCGACGATAAAAACTggg                                                                                                     <  1:53374/251‑1 (MQ=255)                                                                                                                                                  tACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGGAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCCGACGAAGCGGCGCTGTAATCTCTGGCACACAGCAACAGGAACGAAAAATGAAGATAGCGACGATAAAAACTGGgcttgct                                                                                               >  1:178896/1‑251 (MQ=255)                                                                                                                                                    cAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACGCCGGGAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCCGACGAAGCGTCGCTGTAATCTCTGGCACACAGCAACAGGAACGAAAAATGAAGATAGCGACGATAAAAACTGGGCTTGCTTc                                                                                             <  2:464340/251‑1 (MQ=255)                                                                                                                                                                       aaaGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGGAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTc                                                                                                                                                                                                       <  1:241550/126‑1 (MQ=255)                                                                                                                                                                       aaaGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGGAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTc                                                                                                                                                                                                       >  2:241550/1‑126 (MQ=255)                                                                                                                                                                                aTTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGGAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCCGACGAAGCGGCGCTGTAATCTCTGGCACACAGCAACAGGAACGAAAAATGAAGATAGCGACGATAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATg                                                                 >  1:475289/1‑251 (MQ=255)                                                                                                                                                                                          cGATTGCTGATGACCAACTCGCTGAAGTGATCGATCTCGTCAGTAAGGCGGCTTACACCGGGAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCCGACGAAGCGGCGCTGTAATCTCTGGCACACAGCAACAGGAACGAAAAATGAAGATAGCGACGATAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTg                                                       >  2:241047/1‑251 (MQ=255)                                                                                                                                                                                          cGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGAAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCCGACGAAGCGGCGCTGTAATCTCTGGCACACAGCAACAGGAACGAAAAATGAAGATAGCGACGATAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTg                                                       >  2:105915/1‑251 (MQ=255)                                                                                                                                                                                                 tgatgaCCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGGAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCCGACGAAGCGGCGCTGTAATCTCTGGCACACAGCAACAGGAACGAAAAATGAAGATAGCGACGATAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGc                                                >  1:146380/1‑251 (MQ=255)                                                                                                                                                                                                  gatgaCCAACTAGATGAAGTGATCGATATCGGCAGTAAGGCGGCTTACACCGGGAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCCGACGAAGCGGCGCTGTAATCTCTGGCACACAGCAACAGGAACGAAAAATGAAGATAGCGACGATAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGcc                                               <  1:126600/251‑1 (MQ=255)                                                                                                                                                                                                        cAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGGAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCCGACGAAGCGGCGCTGTAATCTCTGGCACACAGCAACAGGAACGAAAAATGAAGATAGCGACGATAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATaaa                                         <  2:146380/251‑1 (MQ=255)                                                                                                                                                                                                          aCTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGGAAAAGTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCCGACGAAGCGGCGCTGTAATCTCTGGCACACAGCAACAGGAACGAAAAATGAAGATAGCGACGATAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGc                                       <  2:346048/251‑1 (MQ=255)                                                                                                                                                                                                            tCGATGAAGTGATAGATATCGTCAGTAAGGCGGCTTACACCGGGAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGAGTCATTCGTATTCGTACCGGCGAAGCCGACGAAGCGGCGCTGTAATCTCTGGCACACAGCAACAGGAACGAAAAATGAAGATAACGACGCTAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTCCGGTGGCCGATAAAGCCg                                     >  1:348704/1‑251 (MQ=255)                                                                                                                                                                                                             cGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGGAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCCGACGAAGCGGCGCTGTAATCTCTGGCACACAGCAACAGGAACGaa                                                                                                                                  >  1:405892/1‑157 (MQ=255)                                                                                                                                                                                                             cGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGGAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCCGACGAAGCGGCGCTGTAATCTCTGGCACACAGCAACAGGAACGaa                                                                                                                                  <  2:405892/157‑1 (MQ=255)                                                                                                                                                                                                                    gTGATCGATATCGTCAGTAAGGCGTCTTACACCGGGAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCTGCATTCGTATTCGTACCGGCGAAGCCGACGAAGCGGCGCTGTAATCTCTGGCACACAGCAACAGGAACGAAAAATGAAGATAGCGACGATAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCGGATAAAGCCGACAATGCg                             <  2:219583/251‑1 (MQ=255)                                                                                                                                                                                                                        tCGATATCGTCAGTAAGGCGGCTTACACCGGGAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCCGACGAAGCGGCGCTGTAATCTCTGGCACACAGCAACAGGAACGAAAAATGAAGATAGCGACGATAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGttt                          >  1:34809/1‑250 (MQ=255)                                                                                                                                                                                                                         cGATATCGTCAGTAAGGCGGCTTACACCGGGAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGGCATTCGTATTCGTACCGGCGAAGCCGACGAAGCGGCTCTGTAATCTCTGGCCCACAGCAACAGGAACGAAAAATGAAGATAGCGACGATAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCGCCTGCGGTGGCCGATAAAGCCGACAATGCGTTTat                        <  2:69297/251‑1 (MQ=255)                                                                                                                                                                                                                           atatCGTCAGTAAGGCGGCTTACACCGGAAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCCGACGAAGCGGCGCTGTAATCTCTGGCACACAGCAACAGGAACGAAAAATGAAGATAGCGACGATAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTatga                      <  1:239635/251‑1 (MQ=255)                                                                                                                                                                                                                                 tCAGTAAGGCGGCTTACACCGGGAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCCGACGAAGCGGCGCTGTAATCTCTGGCACACAGCAACAGGAACGAAAAATGAAGATAGCGACGATAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGAttt                >  1:382980/1‑251 (MQ=255)                                                                                                                                                                                                                                    gTAAGGCGGCTTACACCGGGAAAATCGGCGACGTCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCCGACGAAGCGGCGCTGTAATCTCTGGCACACAGCAACAGGAACGAAAAATGAAGATAGCGACGATAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTa             <  1:422773/251‑1 (MQ=255)                                                                                                                                                                                                                                      aaGGCGGCTTACACCGGGAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCCGACGAAGCGGCGCTGTAATCTATGGCACACAGCCACAGGAACGACAAATGAAGATAGCGACGATAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTACTGGCTGCACCTGCGGTGGCCGCTAAAGCCGACAATGCGTTTATGATAATTTGtcct           >  2:285106/1‑248 (MQ=255)                                                                                                                                                                                                                                       aGGCGGCTTACACCGGGAAAATTGGCGACGGCAACATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCCGACGAAGCGGCGCTGTAATCTCTGGCACACAGCAACAGGAACGAAAAATGAAGATAGCGACGATAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTg          >  1:184570/1‑251 (MQ=255)                                                                                                                                                                                                                                               tACACCGGGAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCCGACGAAGCGGCGCTGTAATCTCTGGCACACAGCAACAGGAACGAAAAATGAAGATAGCGACGATAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTg  >  1:433004/1‑251 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   GCTGGTGACCGTGATAATCAAACCATTCAAGCTGGAAGACGTTCGTGAAGCGTTATCTTCCATTGGTATTCAGGGCCTGACCGTCACCGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGAAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCCGACGAAGCGGCGCTGTAATCTCTGGCACACAGCAACAGGAACGAAAAATGAAGATAGCGACGATAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTG  >  NC_000913/472603‑473091

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