Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,843,898	T→C	Y325C (TAC→TGC)	ynjI ←	inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,843,898	0	T	C	98.4%	238.2 / ‑3.7	64	Y325C (TAC→TGC)	ynjI	inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (40/22);  total (40/24)
Fisher's exact test for biased strand distribution p-value = 3.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

ACAAATGGGCACAATTCATTGCAGTTACGCTCTAATGTAGGCCGGGCAAGCGCAGCGCCCCCGGCAAAATTTCAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGAACAACTCGTGCG  >  NC_000913/1843752‑1844045                                                                                                                                                   |                                                                                                                                                    aCAAATGGGCACAATTCATTGCAGTTACGCTCTAATGTAGGCCGGGCAAGCGCAGCGCCCCCGGCAAAATTTCAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTg                                                                                                                                                 >  1:1058988/1‑151 (MQ=255)            caATTCATTGCAGTTACGCTCTAATGTAGGCCGGGCAAGCGCAGCGCCCCCGGCAAAATTTCAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACt                                                                                                                                      <  1:710265/151‑1 (MQ=255)            caATTCATTGCAGTTACGCTCTAATGTAGGCCGGGCAAGCGCAGCGCCCCCGGCAAAATTTCAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACt                                                                                                                                      >  1:88423/1‑151 (MQ=255)             aaTTCATTGCAGTTACGCTCTAATGTAGGCCGGGCAAGCGCAGCGCCCCCGGCAAAATTTCAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTg                                                                                                                                     >  1:920364/1‑151 (MQ=255)                 cATTGCAGTTACGCTCTAATGTAGGCCGGGCAAGCGCAGCGCCCCCGGCAAAATTTCAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTg                                                                                                                                 >  2:243443/1‑151 (MQ=255)                    tGCAGTTACGCTCTAATGTAGGCCGGGCAAGCGCAGCGCCCCCGGCAAAATTTCAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCg                                                                                                                              >  2:26977/1‑151 (MQ=255)                    tGCAGTTACGCTCTAATGTAGGCCGGGCAAGCGCAGCGCCCCCGGCAAAATTTCAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCg                                                                                                                              >  2:399815/1‑151 (MQ=255)                       aGTTACGCTCTAATGTAGGCCGGGCAAGCGCAGCGCCCCCGGCAAAATTTCAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTaa                                                                                                                           <  2:873905/151‑1 (MQ=255)                       aGTTACGCTCTAATGTAGGCCGGGCAAGCGCAGCGCCCCCGGCAAAATTTCAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTa                                                                                                                            >  1:276604/1‑150 (MQ=255)                       aGTTACGCTCTAATGTAGGCCGGGCAAGCGCAGCGCCCCCGGCAAAATTTCAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTa                                                                                                                            >  2:656232/1‑150 (MQ=255)                             gCTCTAATGTAGGCCGGGCAAGCGCAGCGCCCCCGGCAAAATTTCAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATa                                                                                                                      >  1:201720/1‑150 (MQ=255)                                 tAATGTAGGCCGGGCAAGCGCAGCGCCCCCGGCAAAATTTCAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACg                                                                                                                 >  1:782519/1‑151 (MQ=255)                                    tGTAGGCCGGGCAAGCGCAGCGCCCCCGGCAAAATTTCAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCtt                                                                                                              >  1:477875/1‑151 (MQ=255)                                     gTAGGCCGGGCAAGCGCAGCGCCCCCGGCAAAATTTCAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTc                                                                                                             <  2:276594/151‑1 (MQ=255)                                      tAGGCCGGGCAAGCGCAGCGCCCCCGGCAAAATTTCAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTc                                                                                                             >  2:874550/1‑150 (MQ=255)                                           cGGGCAAGCGCAGCGCCCCCGGCAAAATTTCAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTc                                                                                                       >  1:1065704/1‑151 (MQ=255)                                             ggCAAGCGCAGCGCCCCCGGCAAAATTTCAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTc                                                                                                       >  1:418265/1‑149 (MQ=255)                                              gCAAGCGCAGCGCCCCCGGCAAAATTTCAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGc                                                                                                    >  1:1057919/1‑151 (MQ=255)                                                aaGCGCAGCGCCCCCGGCAAAATTTCAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCAt                                                                                                  >  2:163783/1‑151 (MQ=255)                                                aaGCGCAGCGCCCCCGGCAAAATTTCAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCAt                                                                                                  >  2:395880/1‑151 (MQ=255)                                                     cagcgcCCCCGGCAAAATTTCAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGtt                                                                                             <  1:199947/151‑1 (MQ=255)                                                     cagcgcCCCCGGCAAAATTTCAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGtt                                                                                             >  1:263371/1‑151 (MQ=255)                                                               ggCAAAATTTCAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGc                                                                                   <  2:925190/151‑1 (MQ=255)                                                                 cAAAATTTCAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCg                                                                                  >  2:305003/1‑150 (MQ=255)                                                                   aaaTTTCAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGAt                                                                                >  2:473718/1‑150 (MQ=255)                                                                       ttCAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATagc                                                                             >  2:717183/1‑149 (MQ=255)                                                                           ggCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGa                                                                       <  1:563120/151‑1 (MQ=255)                                                                            gCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATagcagc                                                                          >  1:922679/1‑147 (MQ=255)                                                                            gCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGat                                                                      <  1:115164/151‑1 (MQ=255)                                                                              gTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGatat                                                                    <  2:769767/151‑1 (MQ=255)                                                                                ttATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGc                                                                  >  1:325757/1‑151 (MQ=255)                                                                                ttATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGc                                                                  >  2:466574/1‑151 (MQ=255)                                                                                        aTTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTaaa                                                           <  1:319724/150‑1 (MQ=255)                                                                                            aaCGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGaaa                                                       <  2:370531/150‑1 (MQ=255)                                                                                             aCGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATc                                                     <  2:607359/151‑1 (MQ=255)                                                                                             aCGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATc                                                     <  2:994784/151‑1 (MQ=255)                                                                                              cGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATcc                                                    <  1:608394/151‑1 (MQ=255)                                                                                               ggATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATccc                                                   >  2:147728/1‑151 (MQ=255)                                                                                               ggATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATccc                                                   >  1:621887/1‑151 (MQ=255)                                                                                                 atgatgCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCAt                                                 <  1:342891/151‑1 (MQ=255)                                                                                                    atgCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTg                                              >  2:414027/1‑151 (MQ=255)                                                                                                         ccccACAGAACATTTCTTATGGGCCAACGGCATATCTTACTGAAGTGCAAACAAACCTGCTTGTCGTAAAGATAACACGCGGCAAGAACAGAATAAGGGAAATGTAGGAGAAAGCAGCAGATATGACAATAAAGaaaaaaaagttgaata                                          <  2:1064553/150‑14 (MQ=255)                                                                                                           ccACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTAttt                                       <  2:1010656/151‑1 (MQ=255)                                                                                                           ccACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTAttt                                       <  2:346657/151‑1 (MQ=255)                                                                                                             aCGGCACATGTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATACCACGCTTCGAGTGCAGCATCCGTTAACTTTCTGCCTTAGCGGCATATATGCCAGTAAAGAACTCCCATTTGATGAttttt                                     <  2:229054/151‑1 (MQ=255)                                                                                                               ggAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTAttttttt                                   >  2:554634/1‑151 (MQ=255)                                                                                                               ggAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTAttttttt                                   >  1:767444/1‑151 (MQ=255)                                                                                                               ggAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTAttttttt                                   <  2:875009/151‑1 (MQ=255)                                                                                                                 aaCATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTg                                  >  2:852848/1‑150 (MQ=255)                                                                                                                   cATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATa                               <  1:317790/151‑1 (MQ=255)                                                                                                                   cATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATa                               <  1:142864/151‑1 (MQ=255)                                                                                                                       tCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAAtc                            <  2:454925/150‑1 (MQ=255)                                                                                                                        cTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAAtctt                          >  2:20817/1‑151 (MQ=255)                                                                                                                        cTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAAtctt                          >  1:663298/1‑151 (MQ=255)                                                                                                                          tATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAAtcttct                        >  2:771836/1‑151 (MQ=255)                                                                                                                           aTGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAAtcttct                        <  2:477863/150‑1 (MQ=255)                                                                                                                             gggCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCg                     >  2:32283/1‑151 (MQ=255)                                                                                                                              ggCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGc                    >  1:331013/1‑151 (MQ=255)                                                                                                                               gCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCt                   <  1:945559/151‑1 (MQ=255)                                                                                                                                   aCGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCg               >  2:101061/1‑151 (MQ=255)                                                                                                                                    cGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCg               >  1:1052378/1‑150 (MQ=255)                                                                                                                                      gCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGaac            >  1:831051/1‑151 (MQ=255)                                                                                                                                              tACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGAACAACTCGTg    >  1:216433/1‑151 (MQ=255)                                                                                                                                                cTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGAACAACTCGTgcg  >  2:430356/1‑151 (MQ=255)                                                                                                                                                   |                                                                                                                                                    ACAAATGGGCACAATTCATTGCAGTTACGCTCTAATGTAGGCCGGGCAAGCGCAGCGCCCCCGGCAAAATTTCAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGAACAACTCGTGCG  >  NC_000913/1843752‑1844045

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