Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	362,541	T→A	A213A (GCA→GCT)	lacY ←	lactose permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	362,541	0	T	A	97.4%	117.0 / NA	38	A213A (GCA→GCT)	lacY	lactose permease
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base A (21/16);  total (21/17)
Fisher's exact test for biased strand distribution p-value = 4.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01

GGGCGTTTTTCCCACCGATGCGATTAATGATCAGTGGCGCAAAGAACATAATCGAGGCGTTAAGTAATTCGCCCATTGTCGTTACGTAGCCAAATACCCGCGTACCCTGTTCACCGGTAGCAAAGAACGAAGTAAAGAAATTAGCAAACTGTTGGTCAAAAACATCGTAGGTGCAGGAAACGCCAATAACATACAGTGACAAAAACCACAGTTTTGGCTGTCTG‑AACAGTTCCAGTGCCAGCTTAAGGCTAAATGCCGAATGGTTGGCACCTACCGCATTGGCAACCGTGGCAGAAGAGGGCGCATCCGTTTTGGCGAAAAAGAGTAAAACGGCGAGGATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGC  >  NC_000913/362306‑362776                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                            gggCGTTTTTCCCACCGCTGCGATTAATGATCAGTGGCGCAAAGAACATAATCGAGGCGTTAAGTAATTCGCCCATTGTCGTTACGTAGACAAATACCCGCGTACCCTGTTCACCGGTAGCCAAGAACGAAGTAAAGAAATTAGCAAACTGTTGGTCAAAAACATCGTAGGTGCAGGAAACGCCAATAACATACAGTGACAAAAACCACAGTTTTGGCTGTCTG‑AACCGTTCCAGAGCCAGCTTAAGGCt                                                                                                                                                                                                                               >  1:215146/1‑250 (MQ=255)                            tGATCAGTGGCGCAAAGAACATAATCGAGGCGTTAAGTAATTCGCCCATTGTCGTTACGTAGCCAAATACCCGCGTACCCTGTTCACCGGTAGCAAAGAACGAAGTAAAGAAATTAGCAAACTGTTGGTCAAAAACATCGTAGGTGCAGGAAACGCCAATAACATACAGTGACAAAAACCACAGTTTTGGCTGTCTG‑AACAGTTCCAGAGCCAGCTTAAGGCTAAATGCCGAATGGTTGGCACCTACCGc                                                                                                                                                                                                    >  1:100505/1‑250 (MQ=255)                             gATCAGTGGCGCAAAGAACATAATCGAGGCGTTAAGTAATTCGCCCATTGTCGTTACGTAGCCAAATACCCGCGTACCCTGTTCACCGGTAGCAAAGAACGAAGTAAAGAAATTAGCAAACTGTTGGTCAAAAACATCGTAGGTGCAGGAAACGCCAATAACATACAGTGACAAAAACCACAGTTTTGGCTGTCTG‑AACAGTTCCAGAGCCAGCTTAAGGCTAAATGCCGAATGGTTGGCACCTACCGCAt                                                                                                                                                                                                  <  1:292929/251‑1 (MQ=255)                               tCAGTGGCGCAAAGAACATAATCGAGGCGTTAAGTAATTCGCCCATTGTCGTTACGTAGCCAAATACCCGCGTACCCTGTTCACCGGTAGCAAAGAACGAAGTAAAGAAATTAGCAAACTGTTGGTCAAAAACATCGTAGGTGCAGGAAACGCCAATAACATACAGTGACAAAAACCACAGTTTTGGCTGTCTG‑AACAGTTCCAGAGCCAGCTTAAGGCTAAATGCCGAATGGTTGGCACCTACCGCATTg                                                                                                                                                                                                >  1:189551/1‑251 (MQ=255)                                             aaCATAATCGAGGCTTTAAGTAATTCGCCCATTGTCTTTACGTATCCAAATACCCGCGTACCCTGTTCACCGGTAGCAAAGAACGAAGTAAAGAAATTAGCAAACTGTTGGTCAAAAACATCGGAGGTGCAGGAAACGCCAATAACATACAGTGACAAAAACCACAGTTTTGGCTGTCTG‑AACAGTTCCATAGCCAGCTTAAGGGTAAATGCCGAATGGTGGGCACCTACCGCAGTGGCAACCGTGGCaga                                                                                                                                                                                  <  1:152201/251‑1 (MQ=255)                                              aCATAATCGAGGCGTTAAGAAATTCGCCCATTGTCGTTACGTAGCCAAATACCCGCGTACCCTGTTCACCGGTAGCAAAGAACGAAGTAAAGAAATTAGCAAACTGTTGGTCAAAAACATCGTAGGTGCAGGAAACGCCAATAACATACAGTGACAAAAACCACAGTTTTGGCTGTCTG‑AACAGTTCCAGAGCCAGCTTAAGGCTAAATGCCGAATGGTTGGCACCTACCGCATTGGCAACCGTGGCaga                                                                                                                                                                                  >  1:319572/1‑250 (MQ=255)                                               cATAATCGAGGCGTTAAGTAATTCGCCCCTTGTCGTTACGTAGCCAAATACCCGCGTACCCTGTTCCCCGGTAGCAACGAACGAAGTAAAGAAATTAGCAAACTGTTGGTCAAAAACATCGTAGGGGCAGGAAACGCCAATCACATACAGCGACACAAACCACAGTTTTGGCTGTCTG‑AACAGTTCCAGAGCCAGCTTAAGGCTAAATGCCGAATGGTTGGCACCTACCGCATTGGCACCCGTGGCaga                                                                                                                                                                                  >  1:282016/1‑249 (MQ=255)                                                                        cccATTGTCGTTACGTAGCCAAATACCCGCGTACCCTGTTCACCGGTAGCAAAGAACGAAGGAAAGAAATTAGCAAACTGTTGGACAAAAACATCGGAGGTGCAGGAAACGCTAATAACATACAGTGACAAAAACCACAGTTTTGGCTGTCTG‑AACAGTTCCAGAGCCAGCTTAAGGCTAAAGGCCGAATGGTTGGCACCTACCGCATTGGCAACCGTGGCAGAAGAGGGCGCATCCGTTTTGGCGaaaa                                                                                                                                                        <  2:268494/250‑1 (MQ=255)                                                                        cccATTGTCGTTACGTAGCCAAATACCCGCGTACACTGTTCACCGGTAGCAAAGAACGAAGTAAAGAAATTAGCAAACTGTTGGTCAAAAACATCGTAGGTGCAGGAAACGCCAATAACATACAGTGACAAAAACCACAGTTTTGGCTGTCTG‑AACAGTTCCAGAGCCAGCTTAAGGCTAAATGCCGAATGGTTGGCACCTACCGCATTGGCAACCGTGGCAGAAGAGGGCGCATCCGTTTTGGCGaaaa                                                                                                                                                        >  1:268494/1‑250 (MQ=255)                                                                             tGTCGTTAAGTAGCCAAATACCCGCGTCCCCTGTTCACCGGTAGCAAAGAACGAAGAAAAGAAATTAGCAAACTGTTGGTCAAAAACATCGTAGGAGCAGGAAACGCAAATAACATACAGTGACAAAAACCACAGTTTTGGCTGTCTG‑AACAGTTCCAGAGCCAGCTTAAGGCTAAATGCCGAATGGTTGGCACCTACCCCATTGGCAACCGTGGCAGAAGAGGGCGCATCCGTTTTGGCGAAAAAGAGt                                                                                                                                                   >  2:258036/1‑250 (MQ=255)                                                                                 gtttcGTAGCCAAATACCCGCGTACCCTGTTCGCCGGTAGCAAAGAACGAAGTAAAGAAATTAGCAAACTGTTGGTCAAAAACAGCGTAGGTGCAGGAAACGCCAATAACATACAGTGACAAAAACCACAGTTTTGGGTGGCTG‑AACAGTTCCAGAGCCAGCTTAAGGCGAAATGCCGAAGTGTTGTCGCGTATCGCATTGGCAACGGTGGCAGAATAGGGCGCATCCGTTTTGGCGAAATAGAGTAAAAc                                                                                                                                              <  2:194713/247‑1 (MQ=255)                                                                                  ttACGTAGCCAAATACCCGCGTACCCTGTTCACCGGTAGCAAAGAACGAAGTAAAGAAATTAGCAAACTGTTGGTCAAAAACATCGTAGGTGCAGGAAACGCCAATAACATACAGTGACAAAAACCACAGTTTTGGCTGTCTG‑AACAGTTCCAGAGCCAGCTTAAGGCTAAATGCCGAATGGTTGGCACCTACCGCATTGGCAACCGTGGCAGAAGAGGGCGCATCCGTTTTGGCGAACAAGAGTAAAACg                                                                                                                                             >  2:186735/1‑251 (MQ=255)                                                                                    aCGTAGCCAAATACCCGCGTACCCTGTTCACCGGTAGCAAAGAACGAAGTAAAGAAATTAGCAAACTGTTGGTCAAAAACATCGTAGGTGCAGGAAACGCCAATAACATACAGTGACAAAAACCACAGTTTTGGCTGTCTG‑AACAGTTCCAGAGCCAGCTTAAGGCTAAATGCCGAATGGTTGGCACCTACCGCATTGGCAACCGTGGCAGAAGAGGGCGCATCCGTTTTGGCGAAAAAGAGTAAAACgg                                                                                                                                            <  1:186735/250‑1 (MQ=255)                                                                                      gTAGCCAAATACCCGCGTACCCTGTTCACCGGTAGCAAAGAACGAAGTAAAGAAATTAGCAAACTGTTGGTCAAAAACATCGTAGGGGCAGGAAACGCCAATAACATACAGTGACAAAAACCACAGTTTTGGCTGTCTG‑AACAGTTCCAGAGCCAGCTTAAGGCTAAATGCCGAATGGTTGGCACCTACAGCATTGGCAACCGTGGCAGAAGAGGGCGCATCCGTTTTGGCGAAAAAGCATAAAACGGCg                                                                                                                                          >  2:389889/1‑250 (MQ=255)                                                                                       tAGCCAAATACCCGCGTACCCTGTTCACCGGTAGCAAAGAACGACGTACAGACATTAGCAAACTGTTGGCCAAACACCTCGTAGGTGCAGGAAACGCCCATAACATACAGTGACAAAAACCCCAGTTTTGGCTGTCTG‑AACAGTTCCAGAGCCAGCTTAAGGCTAAATGCCGAATGGTTGGCACCTACCGCATTGGCAACCGTGGCAGAAGAGGGCGCATCCGTTTTGGCGAAACAGAGTAAAACGGCg                                                                                                                                          >  1:389546/1‑249 (MQ=255)                                                                                                     cgTACCCTGTTAACCGGTAGCAAAGCACGAAGTAAAGAAATTAGCAAACTGGTGGTCAAAAACATCGTAGGTGCAGGAAACGCCAAGAACATACAGTGACAAAAACCACAGTTTTGGCTGTCTG‑ACCAGTTCCAGTGACAGCTTAAGGCTAAATGCCGAATTGTTGGCACCTACCGCATTGGCAACCGTGGCAGAAGAGGGCGCATCCGTTTGGGCGAAAAAGAGTAAAACGGCGATGATGAGTGCACAGc                                                                                                                          <  2:51707/251‑1 (MQ=255)                                                                                                         cccTGTTCACCGGTAGCAAAGAACGAAGTAAAGAAATTAGCAAACTGTTGGTCAAAAACATCGTAGGTGCAGGAAACGCCAATAACATACAGTGACAAAAACCACAGTTTTGGCTGTCTG‑AACAGTTCCAGAGCCAGCTTAAGGCTAAATGCCGAATGGTTGGCACCTACCGCATTGGCAACCGTGGCAGAAGAGGGCGCATCCGTTTTGGCGAAAAAGAGTAAAACGGCGAGGATGAGTGCACAGCCaga                                                                                                                      <  1:297816/251‑1 (MQ=255)                                                                                                         cccTGTTCACCGGTAGCAAAGAACGAAGTAAAGAAATGAGCAAACTGTTGGTCAAAAACATCGTAGGTGCAGGAAACGCCAATAACAGACAGTGACAAAAACCACAGTTTTGGCTGTCTG‑AACAGTTCCAGAGCCAGCTTAAGGCTAAATGCCGAATGGTTGGCACCTACCGCATTGGCAACCGTGGCAGAAGAGGGCGCATCCGTTTTGGCGAAAAAGAGTAAAACGGCGAGGATGAGTGCACAGCCag                                                                                                                       <  1:256108/250‑1 (MQ=255)                                                                                                               tCACCGGTAGCAAAGAACGAAGTAAAGAAATTAGCAAACTGTTGGTCAAAAACATCGTAGGTGCAGGAAACGCCAATAACATACAGTGACAAAAACCACAGTTTTGGCTGTCTG‑AACAGTTCCAGAGCCAGCTTAAGGCTAAATGCCGAATGGTTGGCACCTACCGCATTGGCAACCGTGGCAGAAGAGGGCGCATCCGTTTTGGCGAAAAAGAGTAAAACGGCGAGGATGAGTGCACAGCCAGAGCcca                                                                                                                 <  2:369936/250‑1 (MQ=255)                                                                                                                       aGCAAAGAAAGAAGGAAAGAAATTAGCAAACTGTTGGTCAAAAACATCGGATGTGCAGGAAACGCCAATAACATACAGTGACAAAAACCGCAGGTTTGGCGGT‑TGTAACAGTTCCAGAGCCAGCTTAAGGGTAAATGCCGAATGGTTGGCACCTACCGCATTGGCAACCGTGGCAGAAGAGGGCGCATCCGTTTTGGCGAAAAAGAGTAAAACTGCGAGGATGATTGCACAGCCAGAGCCTAGCCAGaaa                                                                                                         <  2:214088/250‑1 (MQ=255)                                                                                                                         cAAAGAACGAAGTAAAGAAATTAGCAAACTGTTGGTCAAAAACATCGTAGGTGCAGGAAACGCCAATAACATACAGTGACAAAAACCACAGTTTTGGCTGTCTG‑AACAGTTCCAGAGCCAGCTTAAGGCTAAATGCCGAATGGTTGGCACCTACCGCATTGGCAACCGTGGCAGAAGAGGGCGCATCCGTTTTGGCGAAAAAGAGTAAAACGGCGAGGATGAGTGCACAGCCAGAGCCCAGCCAGaaaaa                                                                                                       >  2:28666/1‑249 (MQ=255)                                                                                                                                 gAAGTAAAGAAATTAGCAAACTGTTGGTCAAAAACATCGTAGGTGCAGGAAACGCCAATAACATACAGTGACAAAAACCACAGTTTTGGCTGTCTG‑AACAGTTCCAGAGCCAGCTTAAGGCTAAATGCCGAATGGTTGGCACCTACCGCATTGGCAACCGTGGCAGAAGAGGGCGCATCCGTTTTGGCGAAAAAGAGTAAAACGGCGAGGATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGatt                                                                                              >  1:61004/1‑251 (MQ=255)                                                                                                                                                    aCTGTTGGTCAAAAACATCGTAGGTGCAGGAAACGCCAATAACATACAGTGACAAAAACCACAGTTTTGGCTGTCTG‑AACAGTTCCAGAGCCAGCTTAAGGCTAAATGCCGAATGGTTGGCACCTACCGCATTGGCAACCGTGGCAGAAGAGGGCGCATCCGTTTTGGCGAAAAAGAGTAAAACGGCGAGGATGAGTGCACAGCCAGAGCCCAGCCAGAAAACACACCGATTATTGATGGTGAACatgatg                                                                           >  1:19287/1‑251 (MQ=255)                                                                                                                                                      tGTTGGTCAAAAACATCGTAGGTGCAGGAAACGCCAATAACATACAGCGACAAAAACCCAAGTTTTGGCTGTCTG‑AACAGTTCCAGAGCCCGCTTAAGGCTAAATGCCGAATGGTTGGCAACTACCGCATTGGCAACCGAGGCAGAAGAGGGCGCATCCGTTTTGGCGAAAAAGCGTAAAACGGCGAGGATGAGTGCAAAGCCAGAGCCCAGCCAGCAAACAAACTGATTATTGATGGAGAACATGACGcc                                                                         >  2:365927/1‑251 (MQ=255)                                                                                                                                                                 aaCATCGTAGGTGCAGGAAACGCCAATAACATACAGTGACAAAAACCACAGTTTTGGCTGTCTG‑AACAGTTCCAGAGCCAGCTTAAGGCTAAATGCCGAATGGTTGGCACCTACCGCATTGGCAACCGTGGCAGAAGAGGGCGCATCCGTTTTGGCGAAAAAGAGTAAAACGGCGAGGATGAGTGCACAGCCAGAGCCCAGCCAGaaaa                                                                                                        <  2:321406/209‑1 (MQ=255)                                                                                                                                                                 aaCATCGTAGGTGCAGGAAACGCCAATAACATACAGTGACAAAAACCACAGTTTTGGCTGTCTG‑AACAGTTCCAGAGCCAGCTTAAGGCTAAATGCCGAATGGTTGGCACCTACCGCATTGGCAACCGTGGCAGAAGAGGGCGCATCCGTTTTGGCGAAAAAGAGTAAAACGGCGAGGATGAGTGCACAGCCAGAGCCCAGCCAGaaaa                                                                                                        >  1:321406/1‑209 (MQ=255)                                                                                                                                                                    aTCGTAGGTGCAGGAAACGCCAATAACATACAGTGACAAAAACCACAGTTTTGGCTGTCTG‑AACAGTTCCAGAGCCAGCTTAAGGCTAAATGCCGAATGGTTGGCACCTACCGCATTGGCAACCGTGGCAGAAGAGGGCGCATCCGTTTTGGCGAAAAAGAGTAAAACGGCGAGGATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGcac                                                           >  1:73630/1‑251 (MQ=255)                                                                                                                                                                          ggTGCAGGAAACGCCAATAACATACAGTGACAAAAACCACAGTTTTGGCTGTCTG‑AACAGTTCCAGAGCCAGCTTAAGGCTAAATGCCGAATGGTTGGCACCTACCGCATTGGCAACCGTGGCAGAAGAGGGCGCATCCGTTTTGGCGAAAAAGAGTAAAACGGCGAGGATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAgcg                                                     >  1:17008/1‑251 (MQ=255)                                                                                                                                                                          ggTGCAGGAAACGCCAATAACATACAGTGACAAAAACCACAGTTTTGGCTGTCTG‑AACAGTTCCAGAGCCAGCTTAAGGCTAAATGCCGAATGGTTGGCACCTACCGCATTGGCAACCGTGGCAGAAGAGGGCGCATCCGTTTTGGCGAAAAAGAGTAAAACGGCGAGGATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAgcg                                                     <  2:61004/251‑1 (MQ=255)                                                                                                                                                                               aGGAAACGCCAATAACATACAGTGACAAAAACCACAGTTTTGGCTGTCTG‑AACAGTTCCAGAGCCAGCTTAAGGCTAAATGCCGAATGGTTGGCACCTACCGCATTGGCAACCGTGGCAGAAGAGGGCGCATCCGTTTTGGCGAAAAAGAGTAAAACGGCGCGGATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGccc                                                  >  1:3860/1‑249 (MQ=255)                                                                                                                                                                                        cAATAACATGCAGTAACAAAATCCCCAGTGTGGGCTGTCTG‑AACAGTGCCAGAGCCAGCGTAAGGCAAAATGCCGAATGGTTGGCTCCTACCGCATTGGCAACCGTGGCAGAAGAGGGCGCATCCGTTTTGGCGAAAAAGAGTAAAACGGCGAGGATGAGTGCACAGCCAGTGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGATCACAGCGCCCAGCCAacaccg                                       <  2:270205/251‑3 (MQ=255)                                                                                                                                                                                            aaCATACAGTGACAAAAACCACAGTTTTGGCTGTCTG‑AACAGTTCCAGAGCCAGCTTAAGGCTAAATGCCGAATGGTTGGCACCTACCGCATTGGCAACCGTGGCAGAAGAGGGCGCATCCGTTTTGGCGAAAAAGAGTAAAACGGCGAGGATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAACACAGCCaa                                   <  2:383171/251‑1 (MQ=255)                                                                                                                                                                                              cATACAGTGACAAAAACCACAGTTTTGGCTGTCTG‑AACAGTTCCAGAGCCAGCTTAAGGCTAAATGCCTAATGGTTGGCGCCTACCGCATTGGCAACCGTGGCAGAAGAGGGCGCATCCGTTTTGGCGAAAAAGAGTAAAACGGCGAGGATGAGTGCACAGCCAGAGcc                                                                                                                   <  1:237844/169‑1 (MQ=255)                                                                                                                                                                                                     tGACAAAAACCACAGTTTTGGCTGTCTG‑AACAGTTCCAGAGCCAGCTTAAGGCTAAATGCCGAATGGTTGGCACCTACCGCATTGGCAACCGTGGCAGAAGAGGGCGCATCCGTTTTGGCGAAAAAGAGTAAACCGGCGAGGATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAACACAGCCAAACATCcgcg                          >  1:372897/1‑251 (MQ=255)                                                                                                                                                                                                                ccaGTTTTGGCTTGCTG‑AACAGTTCCAGAGCCAGCTTAAGGCTAAATGCCGAATGGTTGGCACGTACCGCATGGGCAACCGTGGCAGAAGAGGGCGCAACCGTTTTGGCGAAAAAGAGTAAAACGGCGAGGATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCTAGGCACACAGCGTCCAGCGAACACAGCGGAGCATCCGCGGGCTACCaaag               <  2:282016/250‑2 (MQ=255)                                                                                                                                                                                                                     tttgggtggatG‑AACAGTTCCAGAGAAAGCTTAAGGCTAAATGCCGAATGGTTGGCACCTCCCGCGTTGACAACCGTAGCAGAAGAGAGCGCATCCGTTTTGGCGAAAAAGAGTAAAACGGCGAGGATGAGTGCACAGCCAGAGCCCAGCCAGAAAACCAACTGATTATTGATGGTGAACATGATGCCGCAAATCGAGGCACACACCGCCCACCCAACACAGCCAACCATCCGCGCGCGACCAAATTCGa           >  2:246985/11‑250 (MQ=255)                                                                                                                                                                                                                        ggCTGTCTG‑AACAGTTCCAGAGCCAGCTTAAGGCTAAATGCCGAATGGTTGCCACCTACCGCATTGGCAACCGTGGCAGAAGAGGGCGCATCCGTTTTGGCGAAAAAGAGTAAAACGGCGAGGATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAAtt       >  2:32089/1‑251 (MQ=255)                                                                                                                                                                                                                                g‑AACAGTTCCAGAGCCAGCTTAAGGCTAAATGCCGAATGGTTGGCACCTACCGCATTGGCAACCGTGGCAGAAGAGGGCGCATCCGTTTTGGCGAAAAAGAGTAAAACGGCGAGGATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGc  <  2:341954/248‑1 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                            GGGCGTTTTTCCCACCGATGCGATTAATGATCAGTGGCGCAAAGAACATAATCGAGGCGTTAAGTAATTCGCCCATTGTCGTTACGTAGCCAAATACCCGCGTACCCTGTTCACCGGTAGCAAAGAACGAAGTAAAGAAATTAGCAAACTGTTGGTCAAAAACATCGTAGGTGCAGGAAACGCCAATAACATACAGTGACAAAAACCACAGTTTTGGCTGTCTG‑AACAGTTCCAGTGCCAGCTTAAGGCTAAATGCCGAATGGTTGGCACCTACCGCATTGGCAACCGTGGCAGAAGAGGGCGCATCCGTTTTGGCGAAAAAGAGTAAAACGGCGAGGATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGC  >  NC_000913/362306‑362776

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