Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_pae_pgi	4,295,641	+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_3_pae_pgi	4,295,640	1	.	C	95.2%	171.1 / 0.6	63	intergenic (+586/+56)	gltP/yjcO	glutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base . (0/3);  new base C (27/33);  total (27/36)
Fisher's exact test for biased strand distribution p-value = 2.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*	NC_000913_3_pae_pgi	4,295,640	2	.	G	95.2%	169.1 / 0.7	63	intergenic (+586/+56)	gltP/yjcO	glutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base . (0/3);  new base G (27/33);  total (27/36)
Fisher's exact test for biased strand distribution p-value = 2.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑A‑‑CCGCA‑‑TCC‑GACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGC  >  NC_000913_3_pae_pgi/4295554‑4295774                                                                                          ||                                                                                                                                          gCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCA‑‑TCC‑GACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAc                                                                            <  2:498487/149‑1 (MQ=1) gCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCA‑‑TCC‑GACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAc                                                                            <  2:962030/149‑1 (MQ=2) gCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCA‑‑TCC‑GACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAc                                                                            >  1:561507/1‑149 (MQ=2)       cTTAACGCGTCTTATCAGGCCTATGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCTCA‑‑TCC‑GACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTgg                                                                      <  1:1582769/149‑1 (MQ=11)        ttAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCACA‑‑TCC‑GACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGt                                                                     >  2:1293369/1‑149 (MQ=255)             gcgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCA‑‑TCC‑GACATCAATGCCTGATGCGACGCTTGCCGCGTCTTTTCAGGCCTTTCTTAATCGTTGGTTAAtt                                                                >  1:652962/1‑149 (MQ=11)              cgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCA‑‑TCC‑GACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAttt                                                               <  1:946031/149‑1 (MQ=255)                tCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCA‑‑TCC‑GACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTCTCTTAACCGTTGGTTACttttt                                                             >  2:335154/1‑149 (MQ=11)                   tATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCA‑‑TCC‑GACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaa                                                          >  1:921183/1‑149 (MQ=255)                    aTCAGGCCTACGCCAGACAGCGCAAAAGACTGCTAAATCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCA‑‑TCC‑GACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTccac                                                         >  2:349866/1‑146 (MQ=14)                        ggCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCA‑‑TCC‑GACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAActct                                                     <  2:561507/149‑1 (MQ=255)                        ggCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCA‑‑TCC‑GACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAActct                                                     >  1:1534380/1‑149 (MQ=255)                         gCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCA‑‑TCC‑GACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCtt                                                    <  1:151663/149‑1 (MQ=255)                         gCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGAACAGTCGGTT‑‑ACGCCGCA‑‑TCC‑GACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCtt                                                    >  1:721391/1‑149 (MQ=17)                           cTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCA‑‑TCC‑GACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTca                                                  <  1:1212609/149‑1 (MQ=255)                                ccAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCA‑‑TCC‑GACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAccc                                             <  1:77294/149‑1 (MQ=255)                                 cAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCA‑‑TCC‑GACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAtt                                                                >  2:1160703/1‑129 (MQ=255)                                 cAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCA‑‑TCC‑GACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAtt                                                                <  1:1160703/129‑1 (MQ=255)                                 cAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCA‑‑TCC‑GACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCg                                            <  1:919382/149‑1 (MQ=255)                                  aGACAGCGCAATAGCCTGATTTAGCTTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCA‑‑TCC‑GACATCAATGCCTGATGCGCGGCATGCCTCGTTCTGAAAGGCTTTGCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCgg                                           <  1:1584077/149‑1 (MQ=14)                                  aGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCA‑‑TCC‑GACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTCTTCAAACTCTTCACACCCgg                                           <  1:1598043/149‑1 (MQ=255)                                    aCAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGGTAAGGCGTTT‑‑ACGCCGCA‑‑TCC‑GACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCCAACTCTTCTCCCCCggga                                         >  1:107423/1‑147 (MQ=17)                                          cAATAGCCTGATTTAGCGTGATTTTGTAAGACGGATAAGGCGTTT‑‑ACGCCGCA‑‑TCC‑GACATCAACACCCGATGCGACGCTGGGCGCGTCTTATCAAGCCTATCTTAACCGTTAGTTAATTTTTCAAACCCTTTACACCCGGTATCAAAc                                   >  2:260937/1‑149 (MQ=17)                                           aaTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCA‑‑TCC‑GACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAcc                                  <  2:250976/149‑1 (MQ=255)                                              aGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCA‑‑TCC‑GACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCtt                               <  1:533092/149‑1 (MQ=255)                                               tgcggATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCA‑‑TCC‑GACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTc                              <  2:842472/145‑1 (MQ=255)                                               gCCTGATTTAGCGTTATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCA‑‑TCC‑GACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTc                              <  2:1372879/149‑1 (MQ=255)                                               gCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCA‑‑TCC‑GACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCATATCTTACCCGTTGGTTAATTTTTCACACTCTTCTACCCCGGTATCCAAcacatc                              >  2:463333/1‑144 (MQ=17)                                               gCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGAG‑TT‑‑A‑‑CCGGAGTTCCGGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTc                              <  1:530884/149‑1 (MQ=21)                                                ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCA‑‑TCC‑GACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTcc                             >  1:256879/1‑149 (MQ=255)                                                ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCA‑‑TCC‑GACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTcc                             <  1:64709/149‑1 (MQ=255)                                                 gtGATTTAGCGTGATTTTGTAGGTCGGATAAGGGGTTT‑‑ACGCCGCA‑‑TAC‑GACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCa                            >  2:872772/2‑149 (MQ=21)                                                 cTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCA‑‑TCC‑GACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAACCCCTTCCa                            >  1:985718/1‑149 (MQ=255)                                                 cTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCA‑‑TCC‑GACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTCCCa                            >  1:1492758/1‑149 (MQ=255)                                                  ggATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTTCC‑‑CCGAA‑‑TCC‑GACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTccct                           <  2:1488143/148‑3 (MQ=21)                                                   gATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCA‑‑TCC‑GACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATa                          <  2:1167496/149‑1 (MQ=255)                                                     tttAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCA‑‑TCC‑GACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACa                        <  1:940202/149‑1 (MQ=255)                                                      tgaGCGGGATTTTGTAGGCCGGAGAAGGCGTTT‑‑ACGCCGCA‑‑TCC‑GACATCAATGCCTGATGCGACGCTGGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAg                       <  1:551014/147‑1 (MQ=21)                                                      taggCGTGATTTTGTAGGTCTGATAAGGCGTTT‑‑ACGCCGCA‑‑TCC‑GACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAg                       <  2:399117/146‑1 (MQ=21)                                                           gTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCA‑‑TCC‑GACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAGACTCTTCACACCCGGTATCAAACCCATCCATAAAACTCAg                  >  1:84698/1‑149 (MQ=21)                                                           gTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCA‑‑TCC‑GACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAACCCCTTCCATACCGCCCac                  >  2:1369854/1‑148 (MQ=21)                                                           atGATTTTGTAGGGCGGATAAGGCGTTT‑‑ACGCCGCA‑‑TCC‑GACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAg                  <  2:595398/148‑1 (MQ=255)                                                            tGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCA‑‑TCC‑GACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAgg                 <  1:1094173/149‑1 (MQ=255)                                                              aTTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCA‑‑TCC‑GACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCaaa                                    <  1:559860/128‑1 (MQ=255)                                                              aTTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCA‑‑TCC‑GACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCaaa                                    >  2:559860/1‑128 (MQ=255)                                                                 ttGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCA‑‑TCC‑GACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCaa            <  1:801860/149‑1 (MQ=255)                                                                 ttGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCA‑‑TCC‑GACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCaa            <  1:718947/149‑1 (MQ=255)                                                                  tGTCGGTCGGATAAGGCGTTTCCG‑‑CCGCA‑‑TCC‑GACATCAATGCCTGATGCGACGTTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAAc           <  1:1310827/149‑1 (MQ=25)                                                                      ggTCGGATAAGGGGTTT‑‑ACGCCGCA‑‑TCC‑GACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGt       <  2:438776/149‑1 (MQ=255)                                                                      ggTCGGATAAGGCGTTT‑‑ACGCCGCA‑‑TCC‑GACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCtt                                                    <  2:1578664/104‑1 (MQ=255)                                                                      ggTCGGATAAGGCGTTT‑‑ACGCCGCA‑‑TCC‑GACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGt       <  1:1367809/149‑1 (MQ=255)                                                                      ggTCGGATAAGGCGTTT‑‑ACGCCGCA‑‑TCC‑GACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATGTTTCAAACTCtt                                                    >  1:1578664/1‑104 (MQ=17)                                                                       gTCGGATAAGGCGTTT‑‑ACGCCGCA‑‑TCC‑GACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATa                          <  1:517480/129‑1 (MQ=255)                                                                       gTCGGATAAGGCGTTT‑‑ACGCCGCA‑‑TCC‑GACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATa                          >  2:517480/1‑129 (MQ=255)                                                                       gTCGGATAAGGCGTTT‑‑ACGCCGCA‑‑TCC‑GACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTg      >  1:234542/1‑149 (MQ=255)                                                                        tCGGATAAGGCGTTT‑‑ACGCCGCA‑‑TCC‑GACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGc     >  1:1301188/1‑149 (MQ=255)                                                                          ggATAAGGCGTTT‑‑ACGCCGCA‑‑TCC‑GACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAg   <  2:691855/149‑1 (MQ=255)                                                                          ggATAAGGCGTTT‑‑ACGCCGCA‑‑TCC‑GACATCAATAAGTGTTGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCCGGTTCACCCAGTGCAg   >  1:793484/1‑149 (MQ=25)                                                                           gATAAGGCGTTT‑‑ACGCCGCA‑‑TCC‑GACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGtt                                                                        >  1:1401324/1‑79 (MQ=11)                                                                           gATAAGGCGTTT‑‑ACGCCGCA‑‑TCC‑GACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGtt                                                                        <  2:1401324/79‑1 (MQ=11)                                                                           gATAAGGCGTTT‑‑ACGCCGCA‑‑TCC‑GACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAgc  <  1:9450/149‑1 (MQ=255)                                                                           gATAAGGCGTTT‑‑ACGCCGCA‑‑TCC‑GACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGGATCAAACCCTTCCATACAGCTCAGGTTAAACCAGTGCAgc  >  2:129417/1‑149 (MQ=25)                                                                           gATAAGGCGTTT‑‑ACGCCGCA‑‑TCC‑GACATCAATGACTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAgc  >  1:651250/1‑149 (MQ=255)                                                                                          ||                                                                                                                                          GCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑A‑‑CCGCA‑‑TCC‑GACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGC  >  NC_000913_3_pae_pgi/4295554‑4295774

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