Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,816,848	C→T	T58M (ACG→ATG)	yicC →	UPF0701 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,816,848	0	C	T	98.4%	239.5 / ‑3.6	62	T58M (ACG→ATG)	yicC	UPF0701 family protein
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (29/32);  total (30/32)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.43e-01

GGCGTGAAATCAAGGGTGAATGGGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAATGCAGAGTGACGAAGGGGAAATCAACCCG  >  NC_000913/3816704‑3816993                                                                                                                                                 |                                                                                                                                                  ggCGTGAAATCAAGGGTGAATGGGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCgg                                                                                                                                             <  2:1222500/151‑1 (MQ=255) ggCGTGAAATCAAGGGTGAATGGGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCgg                                                                                                                                             >  1:775318/1‑151 (MQ=255)    gTGAAATCAAGGGTGAATGGGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCATGATGCGCGGTaa                                                                                                                                          <  1:562533/151‑1 (MQ=255)     tGAAATCAAGGGTGAATGGGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGACGCGCGGTaa                                                                                                                                          >  2:414166/1‑150 (MQ=255)         aTCAAGGGTGAATGGGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTgg                                                                                                                                     <  1:677434/151‑1 (MQ=255)           cAAGGGTGAATGGGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGa                                                                                                                                    <  2:83099/150‑1 (MQ=255)            aaGGGTGAATGGGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGa                                                                                                                                    <  2:272230/149‑1 (MQ=255)            aaGGGTGAATGGGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAAt                                                                                                                                  <  2:555452/151‑1 (MQ=255)                   aaTGGGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCt                                                                                                                           >  2:241506/1‑151 (MQ=255)                        gggAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCt                                                                                                                      <  1:733508/151‑1 (MQ=255)                        gggAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCt                                                                                                                      <  1:633267/151‑1 (MQ=255)                           aGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATg                                                                                                                   >  1:165112/1‑151 (MQ=255)                            gcgcACCCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATg                                                                                                                   <  2:686787/150‑1 (MQ=255)                             cgcAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAg                                                                                                                 <  1:670145/151‑1 (MQ=255)                              gcAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGc                                                                                                                <  1:1053099/151‑1 (MQ=255)                                aaCCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGcc                                                                                                               >  2:35674/1‑150 (MQ=255)                                 aCCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAg                                                                                                             <  2:1101079/151‑1 (MQ=255)                                        aaaTGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAg                                                                                                      <  1:835543/151‑1 (MQ=255)                                          aTGCGCTCGGTAGACCAGCGTTATCTGGAAACTTACTTGCGACTGCCGTAGCAGGTCCGTAGCCTTGAGCCTATCGGTGGCTAGGGTATTCGTTCTCGCCTGATGCACGGTAAATGTGAATGTTCCCTGCGCTATGAGCCAGATGTTAgcg                                                                                                    <  2:1195209/151‑1 (MQ=255)                                          aTGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAgcg                                                                                                    >  1:1097429/1‑151 (MQ=255)                                          aTGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAgcg                                                                                                    >  1:174709/1‑151 (MQ=255)                                                 cagTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAgg                                                                                             <  2:601032/149‑1 (MQ=255)                                                 cGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCaa                                                                                               <  2:736037/149‑1 (MQ=255)                                                  ggTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGt                                                                                            >  1:493488/1‑151 (MQ=255)                                                     aaaCCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGa                                                                                          <  1:836560/150‑1 (MQ=255)                                                       aCCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCt                                                                                       <  2:153841/151‑1 (MQ=255)                                                        ccAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTg                                                                                      <  1:217253/151‑1 (MQ=255)                                                            cGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATcc                                                                                  >  1:802951/1‑151 (MQ=255)                                                                aTCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCa                                                                               <  2:383899/150‑1 (MQ=255)                                                                     gAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGaaa                                                                         <  2:258745/151‑1 (MQ=255)                                                                             cttTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGc                                                                  <  2:675017/150‑1 (MQ=255)                                                                             cttTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCt                                                                 >  2:729250/1‑151 (MQ=255)                                                                               ttCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTa                                                                >  1:455313/1‑150 (MQ=255)                                                                                   tCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAAc                                                             >  2:343883/1‑149 (MQ=255)                                                                                   tCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAAc                                                             >  1:1023023/1‑149 (MQ=255)                                                                                    cTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGc                                                          >  2:244223/1‑151 (MQ=255)                                                                                      gCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTg                                                        <  2:1175147/151‑1 (MQ=255)                                                                                           aGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAAc                                                   >  1:517157/1‑151 (MQ=255)                                                                                              aGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGc                                                >  1:1221701/1‑151 (MQ=255)                                                                                                 tCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCcg                                              <  1:380163/150‑1 (MQ=255)                                                                                                  ccGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCcgcg                                            >  2:97850/1‑151 (MQ=255)                                                                                                   cGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCcgcg                                            >  2:940216/1‑150 (MQ=255)                                                                                                      aGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACt                                        >  2:640664/1‑151 (MQ=255)                                                                                                            gAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTaa                                  <  2:740715/151‑1 (MQ=255)                                                                                                             aaCCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTaaa                                 <  2:1097399/151‑1 (MQ=255)                                                                                                             aaCCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTaa                                  >  1:700757/1‑150 (MQ=255)                                                                                                                 tGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTAGTAACTGCCGCGAACTGGGTAAAAATg                             >  2:982354/1‑151 (MQ=255)                                                                                                                   tCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGCAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAATGc                            >  1:58310/1‑150 (MQ=255)                                                                                                                     gTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAATGc                            >  2:158237/1‑148 (MQ=255)                                                                                                                     gTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAATGc                            >  1:1147149/1‑148 (MQ=255)                                                                                                                     gTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAATGc                            >  1:773061/1‑148 (MQ=255)                                                                                                                       tCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAATGCagag                        <  2:394855/150‑1 (MQ=255)                                                                                                                           gAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAATGCAGAGTGACg                   <  1:88313/151‑1 (MQ=255)                                                                                                                             gCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAATGCAGAGTGACGaa                 <  2:232683/151‑1 (MQ=255)                                                                                                                             gCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAATGCAGAGTGACGa                  <  2:666403/150‑1 (MQ=255)                                                                                                                               gTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAATGCAGAGTGACGAAgg               <  2:638031/151‑1 (MQ=255)                                                                                                                                 aTTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAATGCAGAGTGACGAAgggg             >  1:766283/1‑151 (MQ=255)                                                                                                                                  ttCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAATGCAGAGTGACGAAgggg             >  2:948902/1‑150 (MQ=255)                                                                                                                                     gTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAATGCAGAGTGACGAAGGGGAAAt         >  1:466700/1‑151 (MQ=255)                                                                                                                                          cGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAATGCAGAGTGACGAAGGGGAAATCAAcc    >  1:208719/1‑151 (MQ=255)                                                                                                                                            ccTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAATGCAGAGTGACGAAGGGGAAATCAACCCg  >  2:147010/1‑151 (MQ=255)                                                                                                                                             cTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAATGCAGAGTGACGAAGGGGAAATCAACCCg  <  1:341743/150‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                  GGCGTGAAATCAAGGGTGAATGGGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAATGCAGAGTGACGAAGGGGAAATCAACCCG  >  NC_000913/3816704‑3816993

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