Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,410,605	C→T	A18V (GCG→GTG)	yjfJ →	PspA/IM30 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,410,605	0	C	T	99.0%	346.4 / NA	97	A18V (GCG→GTG)	yjfJ	PspA/IM30 family protein
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (52/43);  total (53/44)
Fisher's exact test for biased strand distribution p-value = 4.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATCCTGCTTGAGATAACCTCGCTGGTAGATAACGCATTGGATCTGGCTGAAATTACAGAAGAATATTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGCGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCTGAAAGATTTACG  >  NC_000913/4410459‑4410752                                                                                                                                                   |                                                                                                                                                    tatCCTGCTTGAGATAACCTCGCTGGTAGATAACGCATTGGATCTGGCTGAAATTACAGAAGAATATTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGga                                                                                                                                                  >  2:1176935/1‑150 (MQ=255)  atCCTGCTTGAGATAACCTCGCTGGTAGATAACGCATTGGATCTGGCTGAAATTACAGAAGAATATTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGga                                                                                                                                                  >  2:394052/1‑149 (MQ=255)   tCCTGCTTGAGATAACCTCGCTGGTAGATAACGCATTGGATCTGGCTGAAATTACAGAAGAATATTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGga                                                                                                                                                  >  1:771042/1‑148 (MQ=255)    ccTGCTTGAGATAACCTCGCTGGTAGATAACGCATTGGATCTGGCTGAAATTACAGAAGAATATTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGgaagaa                                                                                                                                              <  1:193975/151‑1 (MQ=255)    ccTGCTTGAGATAACCTCGCTGGTAGATAACGCATTGGATCTGGCTGAAATTACAGAAGAATATTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGgaagaa                                                                                                                                              <  2:497850/151‑1 (MQ=255)    ccTGCTTGAGATAACCTCGCTGGTAGATAACGCATTGGATCTGGCTGAAATTACAGAAGAATATTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGgaaga                                                                                                                                               >  2:43854/1‑150 (MQ=255)     cTGCTTGAGATAACCTCGCTGGTAGATAACGCATTGGATCTGGCTGAAATTACAGAAGAATATTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGAGGAAATCGTTTATCTCCCAGGTGGAAGAAt                                                                                                                                             >  2:700577/1‑151 (MQ=255)      tGCTTGAGATAACCTCGCTGGTAGATAACGCATTGGATCTGGCTGAAATTACAGAAGAATATTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATc                                                                                                                                            >  2:214771/1‑151 (MQ=255)       gCTTGAGATAACCTCGCTGGTAGATAACGCATTGGATCTGGCTGAAATTACAGAAGAATATTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATcc                                                                                                                                           <  1:902999/151‑1 (MQ=255)        cTTGAGATAACCTCGCTGGTAGATAACGCATTGGATCTGGCTGAAATTACAGAAGAATATTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATcc                                                                                                                                           >  1:629293/1‑150 (MQ=255)          tGAGATAACCTCGCTGGTAGATAACGCATTGGATCTGGCTGAAATTACAGAAGAATATTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCAt                                                                                                                                         <  1:1155542/150‑1 (MQ=255)           gagaTAACCTCGCTGGTAGATAACGCATTGGATCTGGCTGAAATTACAGAAGAATATTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTg                                                                                                                                       >  2:275931/1‑151 (MQ=255)           gagaTAACCTCGCTGGTAGATAACGCATTGGATCTGGCTGAAATTACAGAAGAATATTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTg                                                                                                                                       >  1:285218/1‑151 (MQ=255)           gagaTAACCTCGCTGGTAGATAACGCATTGGATCTGGCTGAAATTACAGAAGAATATTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTg                                                                                                                                       <  1:1212669/151‑1 (MQ=255)           gagaTAACCACGCTGGTAGATAACGCATTGGATCTGGCTGAAATTACAGAAGAATATTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCAtt                                                                                                                                        <  2:250381/150‑1 (MQ=255)            agaTAACCTCGCTGGTAGATAACGCATTGGATCTGGCTGAAATTACAGAAGAATATTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTga                                                                                                                                      <  2:868893/151‑1 (MQ=255)                aaCCTCGCTGGTAGATAACGCATTGGATCTGGCTGAAATTACAGAAGAATATTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTACAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCAAGGTGGAAGAATCCATTgaaga                                                                                                                                   >  1:184485/1‑150 (MQ=255)                aaCCTCGCTGGTAGATAACGCATTGGATCTGGCTGAAATTACAGAAGAATATTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTgaaga                                                                                                                                   >  1:122795/1‑150 (MQ=255)                  ccTCGCTGGTAGATAACGCATTGGATCTGGCTGAAATTACAGAAGAATATTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAAc                                                                                                                                <  2:987894/151‑1 (MQ=255)                  ccTCGCTGGTAGATAACGCATTGGATCTGGCTGAAATTACAGAAGAATATTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAAc                                                                                                                                >  2:1018634/1‑151 (MQ=255)                   cTCGCTGGTAGATAACGCATTGGATCTGGCTGAAATTACAGAAGAATATTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAAcc                                                                                                                               >  2:902105/1‑151 (MQ=255)                     cGCTGGTAGATAACGCATTGGATCTGGCTGAAATTACAGAAGAATATTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCa                                                                                                                             <  2:175831/151‑1 (MQ=255)                         ggTAGATAACGCATTGGATCTGGCTGAAATTACAGAAGAATATTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGc                                                                                                                         <  2:685479/151‑1 (MQ=255)                              aTAACGCATTGGATCTGGCTGAAATTACAGAAGAATATTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTgcg                                                                                                                    <  1:1141299/151‑1 (MQ=255)                              aTAACGCATTGGATCTGGCTGAAATTACAGAAGAATATTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTgc                                                                                                                     <  1:48330/150‑1 (MQ=255)                              aTAACGAATTGGATCTGGCTGAAATTACAGAAGAATATTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTgc                                                                                                                     >  1:925570/1‑150 (MQ=255)                               tAACGCATTGGATCTGGCTGAAATTACAGAAGAATATTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTgcgc                                                                                                                   >  1:6035/1‑151 (MQ=255)                                  cGCATTGGATCTGGCTGAAATTACAGAAGAATATTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATg                                                                                                                <  1:433150/151‑1 (MQ=255)                                     aTTGGATCTGGCTGAAATTACAGAAGAATATTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTg                                                                                                             <  1:82988/151‑1 (MQ=255)                                        ggATCTGGCTGAAATTACAGAAGAATATTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGa                                                                                                           <  2:68794/150‑1 (MQ=255)                                        ggATCTGGCTGAAATTACAGAAGAATATTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGa                                                                                                           >  2:1149673/1‑150 (MQ=255)                                           tCTGGCTGAAATTACAGAAGAATATTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAAc                                                                                                         >  1:828728/1‑149 (MQ=255)                                              ggCTGAAATTACAGAAGAATATTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCat                                                                                                    >  1:135529/1‑151 (MQ=255)                                               gCTGAAATTACAGAAGAATATTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCat                                                                                                    >  1:347827/1‑150 (MQ=255)                                                cTGAAATTACAGAAGAATATTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCatat                                                                                                  >  2:1186888/1‑151 (MQ=255)                                                   aaaTTACAGAAGAATATTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCg                                                                                               >  1:355995/1‑151 (MQ=255)                                                        aCAGAAGAATATTCGCACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACg                                                                                          <  1:576215/151‑1 (MQ=255)                                                        aCAGAAGAATATTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACg                                                                                          >  2:421820/1‑151 (MQ=255)                                                          agaagaATATTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCt                                                                                        >  1:220875/1‑151 (MQ=255)                                                          agaagaATATTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCt                                                                                        >  1:625101/1‑151 (MQ=255)                                                             agaaTATTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTaaa                                                                                     <  1:298845/151‑1 (MQ=255)                                                              gaaTATTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTaa                                                                                      >  2:160526/1‑149 (MQ=255)                                                                atatTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCt                                                                                  >  1:1085798/1‑151 (MQ=255)                                                                  atTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTg                                                                                 >  1:1042594/1‑150 (MQ=255)                                                                  atTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGa                                                                                <  1:305677/151‑1 (MQ=255)                                                                   ttCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGa                                                                                >  2:561377/1‑150 (MQ=255)                                                                    tCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAAc                                                                              >  1:1031812/1‑151 (MQ=255)                                                                     cacaCTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAAc                                                                              <  1:839165/150‑1 (MQ=255)                                                                     cacaCTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACt                                                                             >  1:150308/1‑151 (MQ=255)                                                                      acacTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTc                                                                            >  2:650442/1‑151 (MQ=255)                                                                       cacTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCg                                                                           <  2:6436/151‑1 (MQ=255)                                                                          tAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGAt                                                                         >  1:753841/1‑150 (MQ=255)                                                                                cAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCg                                                                  <  2:510609/151‑1 (MQ=255)                                                                                 aaGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCgg                                                                 <  1:701741/151‑1 (MQ=255)                                                                                  aGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGa                                                                <  1:598121/151‑1 (MQ=255)                                                                                      cGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAAt                                                            <  1:59484/151‑1 (MQ=255)                                                                                       gggTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAAtc                                                           >  2:11135/1‑151 (MQ=255)                                                                                         gTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAAtctc                                                         >  2:45875/1‑151 (MQ=255)                                                                                            tttATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTcgcg                                                      >  2:535348/1‑151 (MQ=255)                                                                                            tttATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTcgcg                                                      >  1:1034599/1‑151 (MQ=255)                                                                                               aTGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGtt                                                    <  2:217878/150‑1 (MQ=255)                                                                                                   gAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATct                                               <  2:377695/151‑1 (MQ=255)                                                                                                   gAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATct                                               <  1:596993/151‑1 (MQ=255)                                                                                                      ttttAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATctgct                                            >  2:497514/1‑151 (MQ=255)                                                                                                       tttAAAAAGTTTATTTGCGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATctgct                                            <  1:644835/150‑1 (MQ=255)                                                                                                       tttAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGCGGAAGAATCCATTGAAGAACCCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAATCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATctgct                                            <  1:1386/150‑1 (MQ=255)                                                                                                         tAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGgc                                         >  2:549238/1‑151 (MQ=255)                                                                                                          aaaaaGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGgc                                         <  2:127430/150‑1 (MQ=255)                                                                                                             aaGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCgg                                     <  1:146198/151‑1 (MQ=255)                                                                                                              aGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCggg                                    >  2:1155591/1‑151 (MQ=255)                                                                                                                tttatttaCGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTg                                  >  2:479808/1‑151 (MQ=255)                                                                                                                 ttatttaCGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTg                                  >  1:468001/1‑150 (MQ=255)                                                                                                                    tttaCGCTGGGGAAATCGTTTATCGCCCAGGTGGAAGAATCCATTGAAGAAGCTCAGGGCGTGCGCATGCTGGAACAGCATATTGGTGATGCTAAAGCGGATCTCGATAAAGCCGGAAAAGCTCGCGTTGGTGTGCTGGCGCGGGTGaaa                               <  1:1064758/150‑1 (MQ=255)                                                                                                                     ttaCGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAAtt                             >  1:638526/1‑151 (MQ=255)                                                                                                                           tGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTc                        >  1:885467/1‑150 (MQ=255)                                                                                                                            tgggAAATCGTTTATCCCCCAGGTGGAAGAATCCATTGAAGAACCACAGGGCGCGCGCATGCTGGAACAGCATATTCGTGACGGTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGTGTTGATCTGCTGGCGCGGGCGAACTTAAggc                        >  1:735350/2‑147 (MQ=255)                                                                                                                            ggggAAATCGTTTATCTCCCAGGGGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGt                         >  1:19879/1‑148 (MQ=255)                                                                                                                             gggAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACg                     >  2:141846/1‑151 (MQ=255)                                                                                                                             gggAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACg                     <  1:1156195/151‑1 (MQ=255)                                                                                                                              ggAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGa                    <  1:1051570/151‑1 (MQ=255)                                                                                                                              ggAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGa                    <  2:98149/151‑1 (MQ=255)                                                                                                                               gAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGAt                   <  1:1018574/151‑1 (MQ=255)                                                                                                                                 aaTCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATa                  >  1:931212/1‑150 (MQ=255)                                                                                                                                  aTCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATa                  >  1:1146091/1‑149 (MQ=255)                                                                                                                                  aTCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAg                <  2:215898/151‑1 (MQ=255)                                                                                                                                    cGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCt              <  1:289862/151‑1 (MQ=255)                                                                                                                                    cGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCt              >  1:804985/1‑151 (MQ=255)                                                                                                                                     gTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCTg             >  1:754243/1‑151 (MQ=255)                                                                                                                                     gTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTAGCGCGGGTGAAATTAAGTCACGATAAGCTg             >  2:667157/1‑151 (MQ=255)                                                                                                                                        tATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCTGaa           <  2:796437/150‑1 (MQ=255)                                                                                                                                             cccAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCTGAAAg         <  2:629311/147‑1 (MQ=255)                                                                                                                                             cccAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCTGAAAGAttt     >  2:151699/1‑151 (MQ=255)                                                                                                                                             cccAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCTGAAAGAttt     <  2:993167/151‑1 (MQ=255)                                                                                                                                             cccAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCTGAAAGAtt      <  1:549895/150‑1 (MQ=255)                                                                                                                                               cAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCTGAAAGATTTAc   <  2:678365/151‑1 (MQ=255)                                                                                                                                                aGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCTGAAAGATTTACg  <  2:359769/151‑1 (MQ=255)                                                                                                                                                aGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCTGAAAGATTTACg  >  2:721155/1‑151 (MQ=255)                                                                                                                                                   |                                                                                                                                                    TATCCTGCTTGAGATAACCTCGCTGGTAGATAACGCATTGGATCTGGCTGAAATTACAGAAGAATATTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGCGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCTGAAAGATTTACG  >  NC_000913/4410459‑4410752

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