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	98.9%	327.7 / ‑4.4	90	A18V (GCG→GTG)	yjfJ	PspA/IM30 family protein
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (51/38);  total (51/39)
Fisher's exact test for biased strand distribution p-value = 4.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

TCCTGCTTGAGATAACCTCGCTGGTAGATAACGCATTGGATCTGGCTGAAATTACAGAAGAATATTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGCGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCTGAAAGATTTACG  >  NC_000913/4410461‑4410752                                                                                                                                                 |                                                                                                                                                    tCCTGCTTGAGATAACCTCGCTGGTAGATAACGCATTGGATCTGGCTGAAATTACAGAAGAATATTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGga                                                                                                                                                  >  1:117095/1‑148 (MQ=255)    tGCTTGAGATAACCTCGCTGGTAGATAACGCATTGGATCTGGCTGAAATTACAGAAGAATATTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATc                                                                                                                                            >  2:908998/1‑151 (MQ=255)        tGAGATAACCTCGCTGGTAGATAACGCATTGGATCTGGCTGAAATTACAGAAGAATATTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCAtt                                                                                                                                        >  1:446424/1‑151 (MQ=255)            aTAACCTCGCTGGTAGATAACGCATTGGATCTGGCTGAAATTACAGAAGAATATTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTgaag                                                                                                                                    <  1:1078496/151‑1 (MQ=255)            aTAACCTCGCTGGTAGATAACGCATTGGATCTGGCTGAAATTACAGAAGAATATTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTga                                                                                                                                      >  2:741693/1‑149 (MQ=255)              aaCCTCGCTGGTAGATAACGCATTGGATCTGGCTGAAATTACAGAAGAATATTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTgaagaa                                                                                                                                  <  1:248849/151‑1 (MQ=255)              aaCCTCGCTGGTAGATAACGCATTGGATCTGGCTGAAATTACAGAAGAATATTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGCGGAAGAATCCATTgaaga                                                                                                                                   <  2:996012/150‑1 (MQ=255)                ccTCGCTGGTAGATAACGCATTGGATCTGGCTGAAATTACAGAAGAATATTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAAc                                                                                                                                <  1:434720/151‑1 (MQ=255)                   cGCTGGTAGATAACGCATTGGATCTGGCTGAAATTACAGAAGAATATTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGTGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAAcc                                                                                                                               <  1:1211707/149‑1 (MQ=255)                   cGCTGGTAGATAACGCATTGGATCTGGCTGAAATTACAGAAGAATATTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAAccc                                                                                                                              <  2:245121/150‑1 (MQ=255)                      tGGTAGATAACGCATTGGATCTGGCTGAAATTACAGAAGAATATTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAgg                                                                                                                           >  1:339798/1‑150 (MQ=255)                         tAGATAACGCATTGGATCTGGCTGAAATTACAGAAGAATATTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGt                                                                                                                       >  2:305700/1‑151 (MQ=255)                           gATAACGCATTGGATCTGGCTGAAATTACAGAAGAATATTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTgc                                                                                                                     <  1:609470/151‑1 (MQ=255)                             tAACGCATTGGATCTGGCTGAAATTACAGAAGAATATTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTgcgc                                                                                                                   >  2:750615/1‑151 (MQ=255)                               aCGCATTGGATCTGGCTGAAATTACAGAAGAATATTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCAt                                                                                                                 <  1:949316/151‑1 (MQ=255)                                   aTTGGATCTGGCTGAAATTACAGAAGAATATTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTg                                                                                                             <  1:489368/151‑1 (MQ=255)                                     tGGATCTGGCTGAAATTACAGAAGAATATTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTgg                                                                                                            <  1:65041/150‑1 (MQ=255)                                      ggATCTGGCTGAAATTACAGAAGAATATTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGaa                                                                                                          <  1:844675/151‑1 (MQ=255)                                       gATCTGGCTGAAATTACAGAAGAATATTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAAc                                                                                                         >  1:38386/1‑151 (MQ=255)                                            ggCTGAAATTACAGAAGAATATTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAAc                                                                                                         <  2:468631/146‑1 (MQ=255)                                            ggCTGAAATTACAGAAGAATATTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAg                                                                                                       <  1:379735/148‑1 (MQ=255)                                            ggCTGAAATTACAGAAGAATATTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCat                                                                                                    <  1:99868/151‑1 (MQ=255)                                             gCTGAAATTACAGAAGAATATTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCat                                                                                                    >  2:1237209/1‑150 (MQ=255)                                              cTGAAATTACAGAAGAATATTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCatat                                                                                                  >  2:961934/1‑151 (MQ=255)                                                 aaaTTACAGAAGAATATTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCg                                                                                               >  1:195208/1‑151 (MQ=255)                                                         gaagaaTATTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCt                                                                                        >  1:1237654/1‑150 (MQ=255)                                                          aagaaTATTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTaa                                                                                      <  1:305694/151‑1 (MQ=255)                                                           agaaTATTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTaaa                                                                                     <  2:629524/151‑1 (MQ=255)                                                           agaaTATTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTaaa                                                                                     <  1:788239/151‑1 (MQ=255)                                                             aaTATTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGc                                                                                   >  1:399537/1‑151 (MQ=255)                                                              atatTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCt                                                                                  >  2:845553/1‑151 (MQ=255)                                                                   cacaCTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGaa                                                                               <  1:9708/149‑1 (MQ=255)                                                                   cacaCTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACt                                                                             <  2:396590/151‑1 (MQ=255)                                                                   cacaCTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACt                                                                             <  2:989748/151‑1 (MQ=255)                                                                    acacTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTc                                                                            >  2:1039484/1‑151 (MQ=255)                                                                      acTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGa                                                                          <  1:750585/151‑1 (MQ=255)                                                                         aaTTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGAt                                                                         <  1:121463/149‑1 (MQ=255)                                                                           ttACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATaa                                                                       >  1:464115/1‑149 (MQ=255)                                                                             aCAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGcc                                                                   >  2:1054782/1‑151 (MQ=255)                                                                              cAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCg                                                                  <  2:508419/151‑1 (MQ=255)                                                                              cAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCg                                                                  <  1:255613/151‑1 (MQ=255)                                                                                  gACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGaa                                                               >  2:298428/1‑150 (MQ=255)                                                                                   aCGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGaaa                                                              >  2:317921/1‑150 (MQ=255)                                                                                    cGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAAt                                                            >  1:1072689/1‑151 (MQ=255)                                                                                      ggTATTTATGGGAATTTTAAAAAGTTTATTTTCGCTGGTGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAAtct                                                          <  2:1103353/151‑1 (MQ=255)                                                                                       gTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAAtctc                                                         >  2:994839/1‑151 (MQ=255)                                                                                       gTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAAtctc                                                         >  1:327287/1‑151 (MQ=255)                                                                                        tATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTcg                                                        >  1:332126/1‑151 (MQ=255)                                                                                            tATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGtt                                                    >  1:1176915/1‑151 (MQ=255)                                                                                            tATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGtt                                                    >  2:93227/1‑151 (MQ=255)                                                                                             aTGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTg                                                   >  2:806425/1‑151 (MQ=255)                                                                                                 gAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATct                                               >  1:221896/1‑151 (MQ=255)                                                                                                 gAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATct                                               <  1:1100483/151‑1 (MQ=255)                                                                                                 gAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATct                                               <  1:550014/151‑1 (MQ=255)                                                                                                  aaTTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATctg                                              >  2:1169119/1‑151 (MQ=255)                                                                                                      ttAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTgg                                          >  2:1235711/1‑151 (MQ=255)                                                                                                       tAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTgg                                          >  1:1164193/1‑150 (MQ=255)                                                                                                         aaaaGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGgcgc                                       >  1:1243521/1‑151 (MQ=255)                                                                                                         aaaaGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGgcgc                                       <  2:712465/151‑1 (MQ=255)                                                                                                         aaaaGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGgcg                                        >  1:5323/1‑150 (MQ=255)                                                                                                         aaaaGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCGCGCGTTGATCTGCTGgcgc                                       <  1:526276/151‑1 (MQ=255)                                                                                                             gTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCggg                                    >  2:1134757/1‑150 (MQ=255)                                                                                                              tttatttaCGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTg                                  >  2:1187338/1‑151 (MQ=255)                                                                                                                 atttaCGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGaa                                <  2:660818/150‑1 (MQ=255)                                                                                                                 atttaCGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATAGTCGTTACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGaaa                               <  2:339806/151‑1 (MQ=255)                                                                                                                   ttaCGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAACTCTCGCGTCGATCTGCTGGCGCGGGTGAAAtt                             >  2:362765/1‑151 (MQ=255)                                                                                                                   ttaCGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAAtt                             >  2:493992/1‑151 (MQ=255)                                                                                                                   ttaCGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAAtt                             >  2:499530/1‑151 (MQ=255)                                                                                                                   ttaCGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAAtt                             >  2:402343/1‑151 (MQ=255)                                                                                                                     aCGCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTa                            >  2:712219/1‑150 (MQ=255)                                                                                                                       gCTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGt                         >  2:57278/1‑151 (MQ=255)                                                                                                                        cTGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTc                        >  1:492915/1‑151 (MQ=255)                                                                                                                         tGGGGAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTc                        >  2:1051815/1‑150 (MQ=255)                                                                                                                          ggggAAATCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCAc                      >  2:711974/1‑151 (MQ=255)                                                                                                                                aTCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATa                  >  2:1117956/1‑149 (MQ=255)                                                                                                                                aTCGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATa                  >  2:359194/1‑149 (MQ=255)                                                                                                                                  cGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGc               <  1:437335/150‑1 (MQ=255)                                                                                                                                  cGTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCt              <  2:1042165/151‑1 (MQ=255)                                                                                                                                   gTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCTg             >  2:701178/1‑151 (MQ=255)                                                                                                                                   gTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCTg             >  1:112121/1‑151 (MQ=255)                                                                                                                                   gTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCTg             <  1:671599/151‑1 (MQ=255)                                                                                                                                   gTTTATCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCTg             <  2:828320/151‑1 (MQ=255)                                                                                                                                       aTCTCCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCTGaa           >  2:508400/1‑149 (MQ=255)                                                                                                                                        tctcCCAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCTGaa           >  2:436035/1‑148 (MQ=255)                                                                                                                                           cccAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCTGAAAg         <  1:1233639/147‑1 (MQ=255)                                                                                                                                           cccAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCTGAAAg         <  1:817684/147‑1 (MQ=255)                                                                                                                                           cccAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCTGAAAGAt       <  2:33593/149‑1 (MQ=255)                                                                                                                                            ccAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCTGAAAGAt       <  1:104366/148‑1 (MQ=255)                                                                                                                                            ccAGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTAACGATAAGCTGAAAGAttt     >  2:550797/1‑150 (MQ=255)                                                                                                                                              aGGTGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCTGAAAGATTTACg  >  2:50417/1‑151 (MQ=255)                                                                                                                                                 |                                                                                                                                                    TCCTGCTTGAGATAACCTCGCTGGTAGATAACGCATTGGATCTGGCTGAAATTACAGAAGAATATTCACACTAATTACAAGGACGGGTATTTATGGGAATTTTAAAAAGTTTATTTACGCTGGGGAAATCGTTTATCTCCCAGGCGGAAGAATCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCTGAAAGATTTACG  >  NC_000913/4410461‑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: ‑