Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,296,381	+GC	intergenic (+587/+55)	gltP → / ← yjcO	glutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,296,380	1	.	C	97.3%	243.2 / ‑1.2	73	intergenic (+586/+56)	gltP/yjcO	glutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base C (43/28);  total (44/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*	NC_000913	4,296,380	2	.	G	97.2%	240.9 / ‑1.4	73	intergenic (+586/+56)	gltP/yjcO	glutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base G (43/27);  total (44/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑GTTTA‑‑C‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCG  >  NC_000913/4296293‑4296515                                                                                          ||                                                                                                                                           tGCAACGCTTACCGCGCCTTGTCAGGCCCACGCCAGCCACCGCAATAGCCTGCTTTCGCGTGATTTTGTAGGTCGGATAAGGC‑GTTTACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAc                                                                             <  2:1140841/150‑1 (MQ=1)  gCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑GTTTACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAcc                                                                            >  3:8626/1‑150 (MQ=9)    gACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑GTTTACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGt                                                                          >  4:164576/1‑150 (MQ=12)     aCGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑GTTTACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGtt                                                                         <  3:35472/150‑1 (MQ=255)      cGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑GTTTACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTg                                                                        <  1:737475/150‑1 (MQ=255)       gCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑GTTTACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTgg                                                                       >  3:541576/1‑150 (MQ=255)       gCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑GTTTACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTgg                                                                       <  3:393963/150‑1 (MQ=255)       gCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑GTTTACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTgg                                                                       <  1:39521/150‑1 (MQ=255)         ttAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑GTTTACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGtt                                                                     >  3:133271/1‑150 (MQ=255)         ttAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑GTTTACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTCGGtt                                                                     >  4:256636/1‑150 (MQ=255)         ttAACGCGTCGTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑GTTTACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGtt                                                                     <  1:746378/150‑1 (MQ=255)          tAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑GTTTACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTa                                                                    >  1:742755/1‑150 (MQ=255)          tAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑GTTTACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATATTAACCGTTGGTTa                                                                    >  4:458349/1‑150 (MQ=255)            aCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑GTTTACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAt                                                                  >  2:162051/1‑150 (MQ=255)              gcgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑GTTTACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAttt                                                                >  2:824385/1‑150 (MQ=255)              gcgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑GTTTACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACAGTTGGTTAAttt                                                                >  3:499482/1‑150 (MQ=255)                gTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑GTTTACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAttttt                                                              <  3:292355/150‑1 (MQ=255)                 tCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑GTTTACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCATATCAGGCCTCTCTTAACCGTAGGTTAACTTTTc                                                             >  1:117346/1‑150 (MQ=11)                  cTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑GTTTACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCa                                                            >  3:185074/1‑150 (MQ=255)                  cTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑GTTTACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCa                                                            >  2:663042/1‑150 (MQ=255)                  cTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑GTTTACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCa                                                            <  2:395243/150‑1 (MQ=255)                   ttATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑GTTTACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaa                                                           >  1:472716/1‑150 (MQ=255)                     aTCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑GTTTACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAc                                                         <  4:408587/150‑1 (MQ=255)                      tCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑GTTTACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAct                                                        >  2:933643/1‑150 (MQ=255)                        aGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑GTTTACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCa                                                            >  2:1046145/1‑144 (MQ=255)                        aGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑GTTTACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCa                                                            <  1:1046145/144‑1 (MQ=255)                        aGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑GTTTACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAActct                                                      >  2:531157/1‑150 (MQ=255)                        aGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑GTTTACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACCTTAGCCGGTGGTTAATTTTTCCATctct                                                      >  2:563805/1‑150 (MQ=11)                         ggCCTTCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑GTTTACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCtt                                                     <  1:977018/150‑1 (MQ=255)                         ggCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑GTTTACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCtt                                                     <  1:858855/150‑1 (MQ=255)                         ggCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑GTTTACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCtt                                                     <  1:1263186/150‑1 (MQ=255)                          gCCTACGCCAGACAGCGCAATAGCATGCTTTAGCGTGAACTTGTAGCGCGGATAAAGC‑GTTTA‑‑CGACG‑CACCCGACATCAATACCAAATGCGACACTTACCGAGTCTGATAAGGCCAATCTTAACCGTCGGAGAATTGATCAAACTCGTc                                                    >  4:15541/1‑150 (MQ=2)                             tACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑GTTTACGC‑‑CG‑CATCCGACATCAATGCCTGATGTGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTcaca                                                 >  2:1008410/1‑150 (MQ=255)                             tACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑GTTTACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTcaca                                                 >  1:630231/1‑150 (MQ=255)                             tACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑GTTTACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTcaca                                                 >  2:98627/1‑150 (MQ=255)                             tACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑GTTTACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTcaca                                                 >  2:755876/1‑150 (MQ=255)                             tACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑GTTTACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTcaca                                                 >  2:1258391/1‑150 (MQ=255)                              aCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑GTTTACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAttt                                                                <  1:74710/134‑1 (MQ=255)                              aCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑GTTTACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAttt                                                                >  2:74710/1‑134 (MQ=255)                              aCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑GTTTACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCacac                                                >  1:145057/1‑150 (MQ=255)                                gCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGACGGCTAAGGC‑GTTTACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAccc                                              >  2:106045/1‑150 (MQ=17)                                 ccAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑GTTTACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCg                                             <  1:219174/150‑1 (MQ=255)                                  cAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑GTTTACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCgg                                            <  3:6631/150‑1 (MQ=255)                                  cAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑GTTTACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCgg                                            >  1:63951/1‑150 (MQ=255)                                   aGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑GTTTACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTa                                                                    >  2:449093/1‑125 (MQ=255)                                   aGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑GTTTACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTa                                                                    <  1:449093/125‑1 (MQ=255)                                   aGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑GTTTACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGt                                           <  1:607107/150‑1 (MQ=255)                                     aCAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑GTTTACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTAt                                         >  1:235431/1‑150 (MQ=255)                                      cAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑GTTTACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATc                                        <  1:437860/150‑1 (MQ=255)                                       aGCGCGATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑GTTTACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCa                                       <  1:20472/150‑1 (MQ=255)                                       aGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑GTTTACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCa                                       >  2:1213544/1‑150 (MQ=255)                                        gcgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑GTTTACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATcca                                      >  1:251303/1‑148 (MQ=255)                                          gcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑GTTTACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAc                                    >  2:996710/1‑150 (MQ=255)                                               aGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑GTTTACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAACCTTc                               >  2:739205/1‑150 (MQ=255)                                                 ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑GTTTACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCa                             >  1:286678/1‑150 (MQ=255)                                                 ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑GTTTACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCa                             >  4:399417/1‑150 (MQ=255)                                                 ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑GTTTACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCa                             <  1:742653/150‑1 (MQ=255)                                                 ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑GTTTACGC‑‑CG‑CATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGGTGGTTAATCTTTCTATCTCGTCACTCCCGGTGTCGGACCCtgcgt                             >  2:742653/1‑146 (MQ=11)                                                  cTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑GTTTACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCAt                            >  3:78596/1‑150 (MQ=255)                                                   tGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑GTTTACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCAtc                           >  2:1265398/1‑149 (MQ=255)                                                   tGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑GTTTACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATa                           >  3:313746/1‑150 (MQ=255)                                                   tGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑GTTTACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATa                           <  1:314093/150‑1 (MQ=255)                                                     aTTTAGCGTGATTTTGTAGGTCGGATAAGGC‑GTTTACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACa                         >  2:1261000/1‑150 (MQ=255)                                                           cGTGATTTTGTAGGTCGGATAAGGC‑GTTTACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAg                   <  1:563805/150‑1 (MQ=255)                                                            gtaaTTTTGTGGGTCGGTTAAGGCGGTTTC‑‑C‑‑CGCCATCCGGCAGCGATGCCTGGGGCGACTCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAgg                  <  1:912144/147‑1 (MQ=37)                                                            gTGATTTTGTAGGTCGGATAAGGC‑GTTTACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAgg                  <  3:488179/150‑1 (MQ=255)                                                            gTGATTTTGTAGGTCGGATAAGGC‑GTTTACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAgg                  <  2:646921/150‑1 (MQ=255)                                                             tGATTTTGTAGGTCGTATACGTC‑GTTTACCC‑‑CG‑CATCCGCCATCAATTCCTGAGCCGACGCTTTCCGCGTCTTATCAGCCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGt                 <  1:1163276/150‑1 (MQ=255)                                                                 tttGTAGGTCGGATAAGGC‑GTTAACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCaa             <  1:20745/150‑1 (MQ=255)                                                                    gTAGGTCGGATAAGGC‑GTTTACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGCATCAAACCCCTCCATACAGCTCACGTTCAACCa          >  2:509949/1‑150 (MQ=25)                                                                     tAGGTCGGATAAGGC‑GTTTACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAg         >  4:145185/1‑150 (MQ=255)                                                                       ggTCGGATAAGGC‑GTTTACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTg       <  3:76250/150‑1 (MQ=255)                                                                            gATAAGGC‑GTTTACGC‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAgcg  >  3:157408/1‑150 (MQ=255)                                                                                          ||                                                                                                                                           TGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑GTTTA‑‑C‑‑CG‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCG  >  NC_000913/4296293‑4296515

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