Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_pae_pgi	4,399,663	T→G	F39C (TTT→TGT)	hfq →	global sRNA chaperone; HF‑I, host factor for RNA phage Q beta replication

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_pae_pgi	4,399,663	0	T	G	100.0%	294.8 / NA	83	F39C (TTT→TGT)	hfq	global sRNA chaperone; HF‑I, host factor for RNA phage Q beta replication
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (40/43);  total (40/43)

GTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATG  >  NC_000913_3_pae_pgi/4399518‑4399803                                                                                                                                                   |                                                                                                                                             gTACAAATAAGCAGATAAGGAAACGAGAGAATGGCCACGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGGATTAAGCTGCAAGGGCAAATCGAGTCTTGTGAt                                                                                                                                          >  2:219036/1‑150 (MQ=255)      aaTAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTGTGATCAGtt                                                                                                                                     <  4:29528/150‑1 (MQ=255)       aTAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTGTGATCAGTTc                                                                                                                                    <  2:471638/150‑1 (MQ=255)          aGCATATAAGGAAAAGAGATACTGGCTAAGTGGCCATCTATACAAGTTCCGTACCTGGACGAACTGGG‑TCGGGTACGTGTTCCAGTATCTATTTGTTTGGTGAGTGGTATTAAGCTGCAAGGGCAATTCGAGGCTTGTGATCAGTCCGTg                                                                                                                                 <  4:459535/150‑1 (MQ=255)            cATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTGTGATCAGTTCGTGAt                                                                                                                               >  3:491906/1‑150 (MQ=255)              tataAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTGTGATCAGTTCGTGATcc                                                                                                                             >  2:508546/1‑150 (MQ=255)               ataAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTGTGATCAGTTCGTGATCCt                                                                                                                            >  2:505368/1‑150 (MQ=255)                 aaGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTGTGATCAGTTCGTGATCCTGt                                                                                                                          >  2:61106/1‑150 (MQ=255)                  aGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTGTGATCAGTTCGTGATCCTGtt                                                                                                                         >  4:647805/1‑150 (MQ=255)                  aGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTGTGATCAGTTCGTGATCCTGtt                                                                                                                         <  1:180469/150‑1 (MQ=255)                   ggAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTGTGATCAGTTCGTGATCCTGTTg                                                                                                                        >  1:487475/1‑150 (MQ=255)                        agagagAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTGTGATCAGTTCGTGATCCTGTTGaaaaa                                                                                                                   <  4:202716/150‑1 (MQ=255)                        agagagAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTGTGATCAGTTCGTGATCCTGTTGaaaaa                                                                                                                   <  1:1011057/150‑1 (MQ=255)                            agaATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTGTGATCAGTTCGTGATCCTGTTGAAAAACACg                                                                                                               <  1:1070090/150‑1 (MQ=255)                              aaTGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTGTGATCAGTTCGTGATCCTGTTGAAAAACACGGt                                                                                                             <  2:487475/150‑1 (MQ=255)                                    tAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTGTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCa                                                                                                       >  1:271669/1‑150 (MQ=255)                                    tAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTGTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCa                                                                                                       >  2:627314/1‑150 (MQ=255)                                    tAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGATGCAGTGCCGGTTTATCAGTCGTAGTGAGTTAGGCTCTGGGGGCTAGTTGGGTCGTGTGATCATTTCGTGTTCCGGTTGAAAAGCTCGGTCAGCCa                                                                                                       >  3:635960/1‑150 (MQ=255)                                       ggggCAAGCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTGTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGAt                                                                                                    <  2:476069/150‑1 (MQ=255)                                        gggCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTGTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATg                                                                                                   >  2:691270/1‑150 (MQ=255)                                            aaTCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTGTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGttt                                                                                               >  4:606682/1‑150 (MQ=255)                                              tCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTGTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTAc                                                                                             <  2:30970/150‑1 (MQ=255)                                                tttACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTGTGATCAGTTCGTGATCCTGTTGAAAAAAACGGCCAGCCAGATGGTTTACaa                                                                                           >  3:343243/1‑150 (MQ=255)                                                   aCAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTGTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTAc                                                                                             >  1:397165/1‑145 (MQ=255)                                                   aCAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTGTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTAc                                                                                             <  2:397165/145‑1 (MQ=255)                                                     aaGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTGTGATCAGTTCGTGATCCTGtt                                                                                                                         <  3:383460/115‑1 (MQ=255)                                                     aaGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTGTGATCAGTTCGTGATCCTGtt                                                                                                                         >  4:383460/1‑115 (MQ=255)                                                     aaGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTGTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCAcg                                                                                      >  3:575464/1‑150 (MQ=255)                                                     aaGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTGTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCAcg                                                                                      <  4:348531/150‑1 (MQ=255)                                                       gATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTGTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCAcgcg                                                                                    >  2:482648/1‑150 (MQ=255)                                                             ttCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTGTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCt                                                                              >  2:397468/1‑150 (MQ=255)                                                                cccgaCGCAGTGCGTTCAGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTGTGATCAGTTCGTGATCCTGtt                                                                                                                         >  1:390298/5‑105 (MQ=255)                                                                cccgaCGCAGTGCGTTCAGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTGTGATCAGTTCGTGATCCTGtt                                                                                                                         <  2:390298/101‑1 (MQ=255)                                                                cTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTGTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACt                                                                           >  1:1152214/1‑150 (MQ=255)                                                                 tGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTGTGATCAGTTCGTGATCCTg                                                                                                                           <  2:241266/101‑1 (MQ=255)                                                                 tGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTGTGATCAGTTCGTGATCCTg                                                                                                                           >  1:241266/1‑101 (MQ=255)                                                                     cGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTGTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACtgttgt                                                                      >  2:305161/1‑150 (MQ=255)                                                                      gCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTGTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTc                                                                     <  3:216207/150‑1 (MQ=255)                                                                      gCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTGTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTc                                                                     <  4:118075/150‑1 (MQ=255)                                                                        aCTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTGTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTccc                                                                   >  1:606524/1‑150 (MQ=255)                                                                             g‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTGTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGtctc                                                              <  4:306310/150‑1 (MQ=255)                                                                             g‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTGTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGtctc                                                              >  4:117600/1‑150 (MQ=255)                                                                               tCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTGTGATCAGTTCGTGATc                                                                                                                              <  2:491105/85‑1 (MQ=255)                                                                               tCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTGTGATCAGTTCGTGATc                                                                                                                              >  1:491105/1‑85 (MQ=255)                                                                                 gggAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTGTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGGTTACAAGCACGCGATTTCTACTGGTGTCCCGGCTCGcc                                                           >  4:298602/1‑150 (MQ=255)                                                                                  ggAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTGTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGccc                                                          >  4:532526/1‑150 (MQ=255)                                                                                   gAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTGTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCg                                                         <  1:994002/150‑1 (MQ=255)                                                                                   gAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTGTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCg                                                         <  1:1124542/150‑1 (MQ=255)                                                                                    aaCGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTGTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCgg                                                        <  3:117600/150‑1 (MQ=255)                                                                                     aCGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTGTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGt                                                       <  2:106647/150‑1 (MQ=255)                                                                                     aCGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTGTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGt                                                       >  3:325253/1‑150 (MQ=255)                                                                                        tgtTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTGTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTtc                                                    <  1:505368/150‑1 (MQ=255)                                                                                        ggtTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTGTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTtc                                                    <  1:1168667/149‑1 (MQ=255)                                                                                          ttCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTGTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCa                                                                                                       <  4:348056/97‑1 (MQ=255)                                                                                          ttCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTGTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCa                                                                                                       >  3:348056/1‑97 (MQ=255)                                                                                             cAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTGTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCtcatc                                               >  3:140430/1‑150 (MQ=255)                                                                                              aGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTGTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCtcatca                                              <  1:1036297/150‑1 (MQ=255)                                                                                              aGTTTCTATTTATTTGGTGAAGGGTATTAAGCTGCAAGGGCAAATCGAGTCTTGTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCtcatca                                              <  4:509648/150‑1 (MQ=255)                                                                                                tttCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTGTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATcaca                                            <  2:116184/150‑1 (MQ=255)                                                                                                tttCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTGTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATcaca                                            <  2:1041241/150‑1 (MQ=255)                                                                                                tttCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTGTGATCAGTTAATAAAAATGTAAAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATcaca                                            <  2:1016977/150‑1 (MQ=255)                                                                                                 ttCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTGTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAg                                           >  1:782720/1‑150 (MQ=255)                                                                                                    tatttattTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTGTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAg                                                                                          <  2:779125/100‑1 (MQ=255)                                                                                                    tatttattTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTGTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAg                                                                                          >  1:779125/1‑100 (MQ=255)                                                                                                    tatttattTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTGTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTaa                                        <  1:508546/150‑1 (MQ=255)                                                                                                       ttatttGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTGTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTaacaa                                     >  4:723669/1‑150 (MQ=255)                                                                                                        tatttGGTGAATGCTATTAAGCTGCAAGGGCAAATCGAGTCTTGTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTaacaac                                    <  2:1036735/150‑1 (MQ=255)                                                                                                                 aaTGGTATTAAGCTGCAAGGGCAAATCGAGTCTTGTGATCAGTTCGTGATACTGAAGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGc                           <  2:660342/150‑1 (MQ=255)                                                                                                                    ggTATTAAGCTGCAAGGGCAAATCGAGTCTTGTGATCAGTTCGTGATCCTGTTGAAAAaca                                                                                                                 <  3:574381/61‑1 (MQ=255)                                                                                                                    ggTATTAAGCTGCAAGGGCAAATCGAGTCTTGTGATCAGTTCGTGATCCTGTTGAAAAaca                                                                                                                 >  4:574381/1‑61 (MQ=255)                                                                                                                    ggTATTAAGCTGCAAGGGCAAATCGAGTCTTGTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGc                                  >  2:166967/1‑140 (MQ=255)                                                                                                                    ggTATTAAGCTGCAAGGGCAAATCGAGTCTTGTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGc                                  <  1:166967/140‑1 (MQ=255)                                                                                                                        ttAAGCTGCAAGGGCAAATCGAGTCTTGTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACca                    <  1:482648/150‑1 (MQ=255)                                                                                                                           aGCTGCAAGGGCAAATCGAGTCTTGTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACcagca                 <  4:655335/150‑1 (MQ=255)                                                                                                                             cTGCAAGGGCAAATCGAGTCTTGTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGt               <  1:691270/150‑1 (MQ=255)                                                                                                                               gCAAGGGCAAATCGAGTCTTGTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTaa             >  4:165719/1‑150 (MQ=255)                                                                                                                                   gggCAAATCGAGTCTTGTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTAc         <  2:140197/150‑1 (MQ=255)                                                                                                                                      cAAATCGAGTCTTGTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACcat      >  2:1113312/1‑150 (MQ=255)                                                                                                                                      cAAATCGAGTCTTGTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACcat      >  3:84379/1‑150 (MQ=255)                                                                                                                                        aaTCGAGTCTTGTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACGAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACcatca    >  4:194486/1‑150 (MQ=255)                                                                                                                                          tCGAGTCTTGTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATg  <  1:627314/150‑1 (MQ=255)                                                                                                                                            gAGTCTTGTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCtca                                                 >  4:530971/1‑101 (MQ=255)                                                                                                                                            gAGTCTTGTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCtca                                                 <  3:530971/101‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                             GTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATG  >  NC_000913_3_pae_pgi/4399518‑4399803

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