Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,225,014	T→A	7.0%	intergenic (‑144/+30)	pbpG ← / ← yohC	peptidoglycan DD‑endopeptidase PbpG/putative inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,225,014	0	T	A	7.0%	260.3 / 7.7	85	intergenic (‑144/+30)	pbpG/yohC	peptidoglycan DD‑endopeptidase PbpG/putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (35/44);  new base A (4/2);  total (39/46)
Fisher's exact test for biased strand distribution p-value = 4.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-01

GGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACTAGAGAAGCTAAGGCCTTCCTCTTTGTTGAT  >  NC_000913/2224866‑2225163                                                                                                                                                     |                                                                                                                                                      ggCATGATGAGCATTCAGATAGTGGTTCAAGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTa                                                                                                                                                      <  1:928450/150‑1 (MQ=255)     tgatgaGCATTCAGATAGTGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGTCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGcc                                                                                                                                                  <  2:657035/150‑1 (MQ=255)        tgaGCATTCAGATAGTGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGc                                                                                                                                               <  2:920150/150‑1 (MQ=255)        tgaGCATTCAGATAGTGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGc                                                                                                                                               <  2:1050733/150‑1 (MQ=255)          aGCATTCAGATAGTGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCaa                                                                                                                                             <  2:937852/150‑1 (MQ=255)             aTTCAGATAGTGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGGCGTAAATGTCAGCAATGc                                                                                                                                          <  2:414195/150‑1 (MQ=255)               tCAGATAGTGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGGCGTAAATGTCAGCAATGCag                                                                                                                                        <  2:1192190/150‑2 (MQ=255)                 aGATAGTGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACg                                                                                                                                      >  1:386146/1‑150 (MQ=255)                  gATAGTGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTATGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGt                                                                                                                                     <  1:481664/150‑1 (MQ=255)                   aTAGTGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTc                                                                                                                                    <  2:73314/150‑1 (MQ=255)                   aTAGTGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTc                                                                                                                                    <  2:65095/150‑1 (MQ=255)                      gTGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCaaa                                                                                                                                 <  2:1217503/150‑1 (MQ=255)                       tGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAg                                                                                                                                >  1:223620/1‑150 (MQ=255)                        ggTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCaacg                                                                                                                                           >  1:856342/1‑136 (MQ=255)                        ggTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCaacg                                                                                                                                           <  2:856342/138‑3 (MQ=255)                         gTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGaa                                                                                                                              >  2:417411/1‑150 (MQ=255)                         gTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGaa                                                                                                                              >  1:297316/1‑150 (MQ=255)                          ttCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAAc                                                                                                                             <  2:506908/150‑1 (MQ=255)                              cgcgCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTc                                                                                                                         >  1:1101713/1‑150 (MQ=255)                              cgcgCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTc                                                                                                                         >  1:404148/1‑150 (MQ=255)                                cgcACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTa                                                                                                                       <  2:887854/150‑1 (MQ=255)                                  cACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGa                                                                                                                     <  1:348411/150‑1 (MQ=255)                                  cACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGa                                                                                                                     >  2:983617/1‑150 (MQ=255)                                        ttGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAAc                                                                                                               <  2:1195070/150‑1 (MQ=255)                                        ttGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAAc                                                                                                               <  1:748951/150‑1 (MQ=255)                                          gcgcACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACaaa                                                                                                                >  1:257704/1‑147 (MQ=255)                                          gcgcACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACaaa                                                                                                                <  2:257704/147‑1 (MQ=255)                                            gcACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATa                                                                                                           <  1:898586/150‑1 (MQ=255)                                            gcACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATa                                                                                                           >  1:1252943/1‑150 (MQ=255)                                            gcACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATCCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACATTCTAGAACAAACGATa                                                                                                           >  2:1118694/1‑150 (MQ=255)                                                 gccgGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGt                                                                                                      <  2:166693/150‑1 (MQ=255)                                                       gTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGcccc                                                                                                >  2:1427338/1‑150 (MQ=255)                                                       gTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGcccc                                                                                                >  2:1422627/1‑150 (MQ=255)                                                           ggATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAg                                                                                            <  1:603982/150‑1 (MQ=255)                                                             aTTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGta                                                                                          <  1:1307140/150‑1 (MQ=255)                                                              tttACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGtat                                                                                         <  1:1328797/150‑1 (MQ=255)                                                              tttACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGtat                                                                                         >  1:431508/1‑150 (MQ=255)                                                               ttACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGtata                                                                                        <  2:1157796/150‑1 (MQ=255)                                                                tACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATaa                                                                                       >  1:1270531/1‑150 (MQ=255)                                                                 aCTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCTCATAGTATAAc                                                                                      <  1:440280/150‑1 (MQ=255)                                                                    gAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTCCATTCCAGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGt                                                                                   >  2:102609/1‑150 (MQ=255)                                                                          aGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGc                                                                             >  1:323383/1‑150 (MQ=255)                                                                             gccGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGag                                                                          >  1:239695/3‑150 (MQ=255)                                                                             gACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGag                                                                          >  1:804861/1‑150 (MQ=255)                                                                             gACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGag                                                                          >  1:239694/1‑150 (MQ=255)                                                                              aCGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGagc                                                                         <  2:116653/150‑1 (MQ=255)                                                                               cGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGagca                                                                        <  1:164048/150‑1 (MQ=255)                                                                                gCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGagcag                                                                       >  1:643875/1‑150 (MQ=255)                                                                                 ccATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGagcagc                                                                      <  1:983617/150‑1 (MQ=255)                                                                                   atcatAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCAc                                                                    >  2:1092412/1‑150 (MQ=255)                                                                                    tcatAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACt                                                                   >  2:1115544/1‑150 (MQ=255)                                                                                     catAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACtt                                                                  <  2:318871/150‑1 (MQ=255)                                                                                      atAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTc                                                                 >  2:57068/1‑150 (MQ=255)                                                                                       tAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTcc                                                                <  2:404148/150‑1 (MQ=255)                                                                                       tAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTcc                                                                >  1:999629/1‑150 (MQ=255)                                                                                       tAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTcc                                                                >  2:1261199/1‑150 (MQ=255)                                                                                        aaCGAGCAAAAAGTACGATAGTCAAAGTAGAATCTTGAGCAAATGCTGCATTGCTTACATTTAACACAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCa                                                               >  2:498624/1‑150 (MQ=255)                                                                                         aCGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAg                                                              <  2:223620/150‑1 (MQ=255)                                                                                         aCGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAg                                                              >  1:124808/1‑150 (MQ=255)                                                                                             gCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTAcc                                                          <  2:999629/150‑1 (MQ=255)                                                                                             gCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTAcc                                                          <  2:1101713/150‑1 (MQ=255)                                                                                                aaaaGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGc                                                       >  1:54979/1‑150 (MQ=255)                                                                                                  aaGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCAc                                                     >  2:1287182/1‑150 (MQ=255)                                                                                                     tGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGcc                                                  <  1:1427338/150‑1 (MQ=255)                                                                                                     tGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGcc                                                  <  1:1422627/150‑1 (MQ=255)                                                                                                     tGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGcc                                                  <  1:1016385/150‑1 (MQ=255)                                                                                                      gCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCa                                                 >  1:197996/1‑150 (MQ=255)                                                                                                      gCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCa                                                 >  1:1271626/1‑150 (MQ=255)                                                                                                             tCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCa                                          <  2:804861/150‑1 (MQ=255)                                                                                                                  ggAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTc                                     <  2:386146/150‑1 (MQ=255)                                                                                                                         ggTGATGAAATGCTTCATTGCTGAAATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCTTAGCCCCATAGTATAACGGTTAGGGCGAGCATCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACTag                              <  2:1252943/149‑1 (MQ=255)                                                                                                                          gTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACTaga                             <  1:986316/150‑1 (MQ=255)                                                                                                                              ggAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACTAGAGAAg                         >  2:1211236/1‑150 (MQ=255)                                                                                                                               ttgacgttgcATTGCTGGCGTAAATGTCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACTAGAGAAGc                        >  1:1409851/8‑150 (MQ=38)                                                                                                                                tgacgttgcATTGCTGGCGTAAATGTCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACTAGAGAAGCt                       >  1:521121/7‑150 (MQ=38)                                                                                                                                 gacgttgcATTGCTGGCGTAAATGTCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACTAGAGAAGCTa                      >  2:916942/6‑150 (MQ=255)                                                                                                                                   cgttgcATTGCTGGCGTAAATGTCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACTAGAGAAGCTAAg                    >  2:937701/4‑150 (MQ=255)                                                                                                                                    gctgcATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACTAGAGAAGCTAAgg                   <  2:731336/150‑1 (MQ=255)                                                                                                                                         aTTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACTAGAGAAGCTAAGGCCTTc              <  1:417411/150‑1 (MQ=255)                                                                                                                                         aTTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACTAGAGAAGCTAAGGCCTTc              >  1:217072/1‑150 (MQ=255)                                                                                                                                          ttGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACTAGAGAAGCTAAGGCCTTcc             <  1:574266/150‑1 (MQ=255)                                                                                                                                           tGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACTAGAGAAGCTAAGGCCTTCct            <  2:217072/150‑1 (MQ=255)                                                                                                                                             cTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACTAGAGAAGCTAAGGCCTTCctct          >  1:1261946/1‑150 (MQ=255)                                                                                                                                                aCATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACTAGAGAAGCTAAGGCCTTCCTCTttg       <  2:1270531/150‑1 (MQ=255)                                                                                                                                                aCATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACTAGAGAAGCTAAGGCCTTCCTCTttg       <  2:197996/150‑1 (MQ=255)                                                                                                                                                     tACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACTAGAGAAGCTAAGGCCTTCCTCTTTGTTGag  >  1:697967/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                      GGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACTAGAGAAGCTAAGGCCTTCCTCTTTGTTGAT  >  NC_000913/2224866‑2225163

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