Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,442,007	T→C	E57G (GAG→GGG)	msrA ←	methionine sulfoxide reductase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,442,007	0	T	C	100.0%	305.8 / NA	85	E57G (GAG→GGG)	msrA	methionine sulfoxide reductase A
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (41/44);  total (41/44)

TCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCTCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGC  >  NC_000913/4441872‑4442134                                                                                                                                        |                                                                                                                                tCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAACCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCca                                                                                                                              <  1:1664475/139‑1 (MQ=255)      aaCAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACAccc                                                                                                                         <  2:1658731/139‑1 (MQ=255)      aaCAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACAccc                                                                                                                         >  2:306455/1‑139 (MQ=255)          aTGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGa                                                                                                                     >  1:219711/1‑139 (MQ=255)          aTGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGa                                                                                                                     >  1:1623720/1‑139 (MQ=255)          aTGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGa                                                                                                                     >  1:1109441/1‑139 (MQ=255)           tGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCca                                                                                                                              >  2:387968/1‑129 (MQ=255)           tGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCca                                                                                                                              <  1:387968/129‑1 (MQ=255)           tGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGaa                                                                                                                    >  2:7864/1‑139 (MQ=255)           tGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGaa                                                                                                                    >  1:373754/1‑139 (MQ=255)           tGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGaa                                                                                                                    >  2:724855/1‑139 (MQ=255)            gCGTACCGCTTCGGCATGACCCGTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCca                                                                                                                              <  1:1585420/128‑1 (MQ=255)            gCGTACCGCTTCGGCATGACCCGTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCca                                                                                                                              >  2:1585420/1‑128 (MQ=255)             cGTACCGCTTCGGCATGACCCGTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAaac                                                                                                                  >  1:1359376/1‑139 (MQ=255)              gTACCGCTTCGGCATGACCCGTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCTAGAACAGACGCCCCACACCCCAGAaaca                                                                                                                 >  2:1001143/1‑139 (MQ=255)                      tCGGCATGACCCGTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCg                                                                                                         >  2:1039320/1‑139 (MQ=255)                      tCGGCATGACCCGTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCg                                                                                                         >  1:90555/1‑139 (MQ=255)                       cGGCATGACCCGTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGc                                                                                                        >  2:642016/1‑139 (MQ=255)                            tGACCCGTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCaaaaa                                                                                                   <  1:533422/139‑1 (MQ=255)                               cccGTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATgg                                                                                                >  2:475304/1‑139 (MQ=255)                                ccGTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCCTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGc                                                                                               <  1:1001143/139‑1 (MQ=255)                                   tATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGc                                                                                               >  2:77622/1‑136 (MQ=255)                                   tATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGc                                                                                               <  1:77622/136‑1 (MQ=255)                                   tATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAAt                                                                                            <  1:1039320/139‑1 (MQ=255)                                         cGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCAt                                                                                      >  2:468286/1‑139 (MQ=255)                                            aGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCca                                                                                                                              >  2:664839/1‑96 (MQ=255)                                            aGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCca                                                                                                                              <  1:664839/96‑1 (MQ=255)                                             gCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCg                                                                                  <  2:1109441/139‑1 (MQ=255)                                             gCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCg                                                                                  >  1:95427/1‑139 (MQ=255)                                             gCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCg                                                                                  <  1:475304/139‑1 (MQ=255)                                              cacaCTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGTTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGt                                                                                 >  1:101189/1‑139 (MQ=255)                                              cacaCTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGt                                                                                 <  2:1623720/139‑1 (MQ=255)                                                     cccGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTAc                                                                          >  2:456405/1‑139 (MQ=255)                                                       cGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACAt                                                                        <  1:306455/139‑1 (MQ=255)                                                       cGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACAt                                                                        <  1:468286/139‑1 (MQ=255)                                                           aaGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGt                                                                    <  1:1302662/139‑1 (MQ=255)                                                           aaGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGt                                                                    <  2:1465968/139‑1 (MQ=255)                                                             gTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGt                                                                            <  2:1080962/129‑1 (MQ=255)                                                             gTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGt                                                                            >  1:1080962/1‑129 (MQ=255)                                                             gTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCa                                                                  >  2:784737/1‑139 (MQ=255)                                                                ggATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTg                                                               <  2:95427/139‑1 (MQ=255)                                                                  aTTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGTTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAg                                                             <  2:101189/139‑1 (MQ=255)                                                                  aTTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAg                                                             <  2:1359376/139‑1 (MQ=255)                                                                   tttGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGt                                                            <  1:784737/139‑1 (MQ=255)                                                                       gcatataGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCAtt                                                                                     >  2:266904/4‑110 (MQ=255)                                                                       gcatataGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCAtt                                                                                     <  1:266904/107‑1 (MQ=255)                                                                             aGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGAc                                                  <  2:219711/139‑1 (MQ=255)                                                                                     gTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCa                                          <  1:724855/139‑1 (MQ=255)                                                                                       ataGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGa                                                          >  2:523964/1‑121 (MQ=255)                                                                                       ataGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGa                                                          <  1:523964/121‑1 (MQ=255)                                                                                       ataGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGc                                        <  1:456405/139‑1 (MQ=255)                                                                                           ccTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCAt                                             <  2:667365/130‑1 (MQ=255)                                                                                           ccTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCAt                                             >  1:667365/1‑130 (MQ=255)                                                                                             tGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCt                                  >  2:1056743/1‑139 (MQ=255)                                                                                              gcggcgGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCca                                                                                                                              >  1:581122/1‑46 (MQ=255)                                                                                              gcggcgGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCca                                                                                                                              <  2:581122/46‑1 (MQ=255)                                                                                               cggcggTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAg                                                             <  2:1319297/110‑1 (MQ=255)                                                                                               cggcggTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAg                                                             >  1:1319297/1‑110 (MQ=255)                                                                                                  cggTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACAcccc                                                                                                                        >  2:1540453/1‑48 (MQ=255)                                                                                                  cggTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACAcccc                                                                                                                        <  1:1540453/48‑1 (MQ=255)                                                                                                     tGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACAt                                                                        <  2:99750/93‑1 (MQ=255)                                                                                                     tGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACAt                                                                        >  1:99750/1‑93 (MQ=255)                                                                                                         gTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCggg                      <  1:642016/139‑1 (MQ=255)                                                                                                             aCGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGgtgt                  <  2:803951/139‑1 (MQ=255)                                                                                                             aCGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGgtgt                  <  1:7864/139‑1 (MQ=255)                                                                                                                 cGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAAcc                                                                                                              <  2:1576827/43‑1 (MQ=255)                                                                                                                 cGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAAcc                                                                                                              >  1:1576827/1‑43 (MQ=255)                                                                                                                 cGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGc                                                                                               >  1:470394/1‑58 (MQ=255)                                                                                                                 cGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGc                                                                                               <  2:470394/58‑1 (MQ=255)                                                                                                                  gggTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGt                                                            <  2:264911/92‑1 (MQ=255)                                                                                                                  gggTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGt                                                            >  1:264911/1‑92 (MQ=255)                                                                                                                   ggTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTc            >  1:288674/1‑139 (MQ=255)                                                                                                                    gTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACAt                                                                        <  2:1397087/78‑1 (MQ=255)                                                                                                                    gTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACAt                                                                        >  1:1397087/1‑78 (MQ=255)                                                                                                                        cTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCAtt                                                                                     <  1:801291/61‑1 (MQ=255)                                                                                                                        cTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCAtt                                                                                     >  2:1473305/1‑61 (MQ=255)                                                                                                                        cTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCAtt                                                                                     >  2:801291/1‑61 (MQ=255)                                                                                                                        cTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCAtt                                                                                     <  1:1473305/61‑1 (MQ=255)                                                                                                                             aGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGc  <  2:373754/139‑1 (MQ=255)                                                                                                                                   gACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGc                                                                                               >  2:153474/1‑40 (MQ=39)                                                                                                                                   gACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGc                                                                                               <  1:153474/40‑1 (MQ=255)                                                                                                                                   gACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTcc                                                                                   >  1:1657112/1‑52 (MQ=255)                                                                                                                                   gACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTcc                                                                                   <  2:1657112/52‑1 (MQ=255)                                                                                                                                   gACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACg                               <  2:577930/104‑1 (MQ=255)                                                                                                                                   gACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACg                               >  1:577930/1‑104 (MQ=255)                                                                                                                                        |                                                                                                                                TCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCTCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGC  >  NC_000913/4441872‑4442134

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