Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,624,357	(T)8→7	coding (141/393 nt)	ydeI ←	hydrogen peroxide resistance OB fold protein, putative periplasmic protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,624,350	0	T	.	100.0%	287.6 / NA	71	coding (148/393 nt)	ydeI	hydrogen peroxide resistance OB fold protein, putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (33/38);  total (33/38)

CCGCAGGAATAACGACATTAATTTCACCGCTCTTATCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCATTTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGCATAAGGCATGACAAGAAAACTTGCTAATACAATGG  >  NC_000913/1624215‑1624484                                                                                                                                        |                                                                                                                                       ccGCAGGAATAACGACATTAATTTCACCGCTCTTATCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCAtttt                                                                                                                                     >  1:1372267/1‑139 (MQ=255)   gCAGGAATAACGACATTAATTTCACCGCTCTTATCGCGAAAAACGTAACGGTCCTCTCC‑TTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCAttttttt                                                                                                                                  <  2:970076/139‑1 (MQ=255)    cAGGAATAACGACATTAATTTCACCGCTCTTATCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGC‑‑GGTTTTTGCGAAATCAACGGTCAtttttttgc                                                                                                                                <  2:1024803/139‑3 (MQ=255)      ggAATAACGACATTAATTTCACCGCTCTTATCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCAtttttttgc                                                                                                                                <  2:523165/139‑3 (MQ=255)         aTAACGACATTAATTTCACCGCTCTTATCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGcatc                                                                                                                            <  2:1724554/139‑1 (MQ=255)                 aTTAATTTCACCGCTCTTATCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTc                                                                                                                    >  1:390395/1‑139 (MQ=255)                    aaTTTCACCGCTCTTATCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCcgc                                                                                                                 >  1:170693/1‑139 (MQ=255)                         cACCGCTCTTATCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAg                                                                                                            <  1:620540/139‑1 (MQ=255)                          aCCGCTCTTATCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGc                                                                                                           >  2:1775525/1‑139 (MQ=255)                          aCCGCTCTTATCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGc                                                                                                           <  2:686649/139‑1 (MQ=255)                           ccGCTCTTATCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGcc                                                                                                          <  1:1506563/139‑1 (MQ=255)                           ccGCTCTTATCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGcc                                                                                                          <  1:154647/139‑1 (MQ=255)                           ccGCTCTTATCGCGAAAAACGTAAAGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGcc                                                                                                          >  1:1275444/1‑139 (MQ=255)                              ctctTATCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGcc                                                                                                          <  1:1883457/136‑1 (MQ=255)                              ctctTATCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGcc                                                                                                          >  2:1883457/1‑136 (MQ=255)                              ctctTATCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCAtc                                                                                                       <  2:902509/139‑1 (MQ=255)                              ctctTATCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCAtc                                                                                                       <  2:491481/139‑1 (MQ=255)                                   aTCGCGAAACACGTAACGGTCCAATCCTTTGTGAGAAATCAAATTACTGCTTAGTAAAACCGAGGCTGCCTCGGGCATGCTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGATGCCGCGATAGCCATCTTCTa                                                                                                  >  1:1222157/1‑139 (MQ=255)                                           aaaCGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGa                                                                                          >  2:825695/1‑139 (MQ=255)                                           aaaCGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGa                                                                                          >  2:1344255/1‑139 (MQ=255)                                                      tCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAg                                                                                                            >  1:1117908/1‑110 (MQ=255)                                                      tCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAg                                                                                                            <  2:1117908/110‑1 (MQ=255)                                                      tCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTgg                                                                               >  2:1432495/1‑139 (MQ=255)                                                                tgtgAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGcc                                                                                                          >  2:1400009/1‑102 (MQ=255)                                                                tgtgAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGcc                                                                                                          <  1:1400009/102‑1 (MQ=255)                                                                tgtgAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATcc                                                                     >  2:1883334/1‑139 (MQ=255)                                                                    agaAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGtt                                                                 >  1:1454518/1‑139 (MQ=255)                                                                       aaTCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAggcgg                                                                                     <  1:1569681/116‑1 (MQ=255)                                                                       aaTCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAggcgg                                                                                     >  2:1569681/1‑116 (MQ=255)                                                                       aaTCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTa                                                              >  2:1252152/1‑139 (MQ=255)                                                                         tCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTaaa                                                            <  1:1775525/139‑1 (MQ=255)                                                                         tCAAATTACCGCGTAGCGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTCGCGCTGCATCCTGTTTTaaa                                                            <  2:1275444/139‑1 (MQ=255)                                                                                           aaaCCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTaaa                                                            >  1:1322690/1‑121 (MQ=255)                                                                                           aaaCCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTaaa                                                            <  2:1322690/121‑1 (MQ=255)                                                                                           aaaCCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGcc                                          <  1:1574966/139‑1 (MQ=255)                                                                                           aaaCCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGcc                                          <  1:825695/139‑1 (MQ=255)                                                                                           aaaCCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGcc                                          <  1:921278/139‑1 (MQ=255)                                                                                            aaCCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCa                                         <  2:1372267/139‑1 (MQ=255)                                                                                                 aaGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCt                                                                           >  2:1716377/1‑100 (MQ=255)                                                                                                 aaGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCt                                                                           <  1:1716377/100‑1 (MQ=255)                                                                                                 aaGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGCa                                    <  1:106195/139‑1 (MQ=255)                                                                                                 aaGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGCa                                    <  1:1252152/139‑1 (MQ=255)                                                                                                  aGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGCa                                    >  2:1227632/1‑138 (MQ=255)                                                                                                  aGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGCa                                    <  1:1227632/138‑1 (MQ=255)                                                                                                   gcgcCATCGTGCATGGTTTTTGCGAAATCAACGGTCAttttttt                                                                                                                                  >  2:1002424/1‑44 (MQ=255)                                                                                                   gcgcCATCGTGCATGGTTTTTGCGAAATCAACGGTCAttttttt                                                                                                                                  <  1:1002424/44‑1 (MQ=255)                                                                                                   gcgcCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAg                                                                                                            <  2:304379/65‑1 (MQ=255)                                                                                                   gcgcCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAg                                                                                                            >  1:304379/1‑65 (MQ=255)                                                                                                   gcgcCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCAtc                                                                                                       >  1:922125/1‑70 (MQ=255)                                                                                                   gcgcCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCAtc                                                                                                       <  2:922125/70‑1 (MQ=255)                                                                                                   gcgcCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGCATa                                  >  2:1688950/1‑139 (MQ=255)                                                                                                      ccATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGCa                                    >  1:1002850/1‑134 (MQ=255)                                                                                                      ccATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGCa                                    <  2:1002850/134‑1 (MQ=255)                                                                                                       cATCGTGCATGGTTTTTGCGAAGTCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGCAt                                   >  2:817211/1‑134 (MQ=255)                                                                                                       cATCGTGCATGGTTTTTGCGAAATCAACGGTCAtttttttgcat                                                                                                                              >  1:1836606/1‑40 (MQ=255)                                                                                                       cATCGTGCATGGTTTTTGCGAAATCAACGGTCAtttttttgcat                                                                                                                              <  2:1836606/44‑5 (MQ=255)                                                                                                       cATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGCAt                                   <  1:817211/134‑1 (MQ=255)                                                                                                         tCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGCATAAGGCAt                            <  1:1344255/139‑1 (MQ=255)                                                                                                         tCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGCATAAGGCAt                            <  1:1883334/139‑1 (MQ=255)                                                                                                                   tttttGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGCATAAGGCATGACAAGaaaa                  <  2:390395/139‑1 (MQ=255)                                                                                                                    ttttGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGCATAAGGCATGACAAGAAAAc                 >  1:1715120/1‑139 (MQ=255)                                                                                                                      ttGCGAAATCAACGGTCA‑TTTTTTTGTATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCa                                         >  1:358857/1‑113 (MQ=255)                                                                                                                      ttGCGAAATCAACGGTCA‑TTTTTTTGTATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCa                                         <  2:358857/113‑1 (MQ=255)                                                                                                                      ttGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGCATAAGGCATGACAAGaaacttg               >  1:758193/1‑135 (MQ=255)                                                                                                                       tGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGCATAAGGCATGACAAGaaaa                  <  1:786089/135‑1 (MQ=255)                                                                                                                       tGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGCATAAGGCATGACAAGaaaa                  >  2:786089/1‑135 (MQ=255)                                                                                                                            aaTCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCAtctt                                                                                                     <  1:954019/47‑1 (MQ=38)                                                                                                                            aaTCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCAtctt                                                                                                     >  2:954019/1‑47 (MQ=38)                                                                                                                             aTCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGa                                                                                                               >  1:1866101/1‑36 (MQ=37)                                                                                                                             aTCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGa                                                                                                               <  2:1866101/36‑1 (MQ=37)                                                                                                                             aTCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCAtc                                                                                                       >  2:1473620/1‑44 (MQ=38)                                                                                                                             aTCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCAtc                                                                                                       <  1:1473620/44‑1 (MQ=38)                                                                                                                             aTCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGgcggcg                                                                                   >  2:488325/1‑64 (MQ=255)                                                                                                                             aTCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGgcggcg                                                                                   <  1:488325/64‑1 (MQ=255)                                                                                                                             aTCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGt                                                                                 >  2:553344/1‑66 (MQ=255)                                                                                                                             aTCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGt                                                                                 <  1:553344/66‑1 (MQ=255)                                                                                                                               cAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGCATAAGGCATGACAAGAAAACTTGCTAATACa      <  2:1454518/139‑1 (MQ=255)                                                                                                                                   ggTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGCATAAGGCATGACAAGAAAACTTGCTAATACAATgg  >  2:465350/1‑139 (MQ=255)                                                                                                                                   ggTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGCATAAGGCATGACAAGAAAACTTGCTAATACAATgg  <  2:863499/139‑1 (MQ=255)                                                                                                                                        |                                                                                                                                       CCGCAGGAATAACGACATTAATTTCACCGCTCTTATCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCATTTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGCATAAGGCATGACAAGAAAACTTGCTAATACAATGG  >  NC_000913/1624215‑1624484

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