Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,618,075	C→T	*234* (TAG→TAA)	hda ←	ATPase regulatory factor involved in DnaA inactivation

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,618,075	0	C	T	100.0%	253.1 / NA	77	*234* (TAG→TAA)	hda	ATPase regulatory factor involved in DnaA inactivation
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (41/36);  total (41/36)

ATCGCTAGTTGAAGCACATTGCCGGATGCGACGCTTGCCACGTCTTATCCGGCCTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGTCATAAATAGCGTGCGCATTTCTCTGTCGAGCCGCTTCAGCAAGAAACGCCCC  >  NC_000913/2617948‑2618210                                                                                                                                |                                                                                                                                        aTCGCTAGTTGAAGCACATTACCGGATGCGACGCTTGCCACGTCTTATCCGGCCTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAg                                                                                                                              >  1:946135/1‑139 (MQ=255)          tGAAGCACATTGCCGGATGCGACGCTTGCCACGTCTTATCCGGCCTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAAtttcttt                                                                                                                     <  1:1174539/139‑1 (MQ=255)          tGAAGCACATTGCCGGATGCGACGCTTGCCACGTCTTATCCGGCCTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAAtttcttt                                                                                                                     >  1:1277616/1‑139 (MQ=255)             aGCACATTGCCGGATGCGACGCTTGCCACGTCTTATCCGGCCTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTcac                                                                                                                  >  1:455509/1‑139 (MQ=255)                  aTTGCCGGATGCGACGCTTGCCACGTCTTATCCGGCCTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACg                                                                                                             >  1:207310/1‑139 (MQ=255)                   ttGCCGGATGCGACGCTTGCCACGTCTTATCCGGCCTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACgg                                                                                                            >  2:1352516/1‑139 (MQ=255)                       cGGATGCGACGCTTGCCACGTCTTATCCGGCCTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATg                                                                                                        >  1:386918/1‑139 (MQ=255)                        ggATGCGACGCTTGCCACGTCTTATCCGGCCTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATgg                                                                                                       <  2:455509/139‑1 (MQ=255)                         gATGCGACGCTTGCCACGTCTTATCCGGCCTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGt                                                                                                      >  2:1098527/1‑139 (MQ=255)                               aCGCTTGCCACGTCTTATCCGGCCTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCtt                                                                                                >  1:824630/1‑139 (MQ=255)                                   ttGCCACGTCTTATCCGGCCTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGt                                                                                            <  1:766494/139‑1 (MQ=255)                                               aTCCGGCCTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAAtttcttt                                                                                                                     <  1:1023736/102‑1 (MQ=255)                                               aTCCGGCCTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGACGTTCGCGCCGCATCGGACAATAAACACCTTATTTACAAGTTCAGAAtttcttt                                                                                                                     >  2:1023736/1‑102 (MQ=255)                                                     ccTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGa                                                                                                           <  1:905340/106‑1 (MQ=255)                                                     ccTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGa                                                                                                           >  2:905340/1‑106 (MQ=255)                                                     ccTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGc                                                                           <  1:842626/138‑1 (MQ=255)                                                     ccTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGc                                                                           >  2:842626/1‑138 (MQ=255)                                                     ccTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCa                                                                          <  1:1486400/139‑1 (MQ=255)                                                     ccTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCa                                                                          <  2:808897/139‑1 (MQ=255)                                                     ccTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCa                                                                          <  2:1483107/139‑1 (MQ=255)                                                        aCGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGa                                                                       >  1:1289222/1‑139 (MQ=255)                                                               cTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACaa                                                                                                                >  1:1515883/1‑91 (MQ=255)                                                               cTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACaa                                                                                                                <  2:1515883/91‑1 (MQ=255)                                                                     aTTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAATCACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAAt                                                                                 >  1:790863/1‑116 (MQ=255)                                                                     aTTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAATCACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAAt                                                                                 <  2:790863/116‑1 (MQ=255)                                                                     aTTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGt                                                                                    <  1:445974/113‑1 (MQ=255)                                                                     aTTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGt                                                                                    >  2:445974/1‑113 (MQ=255)                                                                     aTTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCa                                                          <  2:1289222/139‑1 (MQ=255)                                                                     aTTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCa                                                          <  2:386918/139‑1 (MQ=255)                                                                     aTTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCa                                                          >  1:1562128/1‑139 (MQ=255)                                                                       tAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTGCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAAc                                                        <  1:1098527/139‑1 (MQ=255)                                                                        aaCTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACg                                                       >  2:813215/1‑139 (MQ=255)                                                                          cTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGTc                                                     >  2:534756/1‑139 (MQ=255)                                                                          cTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGTc                                                     >  1:414314/1‑139 (MQ=255)                                                                          cTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTAATTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCa                                                          <  2:466103/134‑1 (MQ=255)                                                                          cTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTAATTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCa                                                          >  1:466103/1‑134 (MQ=255)                                                                               cgTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCa                                                          >  1:222056/1‑129 (MQ=255)                                                                               cgTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCa                                                          <  2:222056/129‑1 (MQ=255)                                                                                   gTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGTCATAAATAGc                                            <  2:1277616/139‑1 (MQ=255)                                                                                   gTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGTCATAAATAGc                                            <  1:1352516/139‑1 (MQ=255)                                                                                    ttCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACa                                                                                                                                     >  2:333256/1‑49 (MQ=255)                                                                                    ttCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACa                                                                                                                                     <  1:333256/49‑1 (MQ=255)                                                                                        gATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGTCATAAATAGCGTgcg                                       <  2:207310/139‑1 (MQ=255)                                                                                            aGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGTCATAAATAGCGTgcg                                       >  2:304561/1‑135 (MQ=255)                                                                                            aGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGTCATAAATAGCGTgcg                                       <  1:304561/135‑1 (MQ=255)                                                                                             ggCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATgg                                                                                                       >  2:35173/1‑70 (MQ=255)                                                                                             ggCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATgg                                                                                                       <  1:35173/70‑1 (MQ=255)                                                                                              gCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACaa                                                                                                                <  1:50972/60‑1 (MQ=255)                                                                                              gCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACaa                                                                                                                >  2:50972/1‑60 (MQ=255)                                                                                              gCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCa                                                                          >  2:638382/1‑98 (MQ=255)                                                                                              gCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCa                                                                          <  1:638382/98‑1 (MQ=255)                                                                                              gCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGTCATAAATAGCGTgcg                                       <  2:757945/133‑1 (MQ=255)                                                                                              gCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGTCATAAATAGCGTgcg                                       >  1:757945/1‑133 (MQ=255)                                                                                              gCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGTCATAAATAGCGTGCGCATTtc                                 >  2:1165609/1‑139 (MQ=255)                                                                                              gCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGTCATAAATAGCGTGCGCATTtc                                 >  1:66042/1‑139 (MQ=255)                                                                                              gCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGTCATAAATAGCGTGCGCATTtc                                 >  1:1114857/1‑139 (MQ=255)                                                                                                gTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGTCATAAATAGCGTGCGCATTtctc                               >  2:22385/1‑139 (MQ=255)                                                                                                   cgcgCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAAtttct                                                                                                                       >  1:807806/1‑48 (MQ=255)                                                                                                   cgcgCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAAtttct                                                                                                                       <  2:807806/48‑1 (MQ=255)                                                                                                   cgcgCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGAt                                                                      <  2:951982/97‑1 (MQ=255)                                                                                                   cgcgCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGAt                                                                      >  1:951982/1‑97 (MQ=255)                                                                                                   cgcgCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGTCATAAATAGCGTGCGCATTTCTCTGt                            >  2:90542/1‑139 (MQ=255)                                                                                                       ccgcATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGTCATAAATAGCGTGCGCATTTCTCTGTCGAg                        <  2:66042/139‑1 (MQ=255)                                                                                                            tCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGTCATAAATAGCGTGCGCATTTCTCTGTCGAGCCGCt                   <  2:824630/139‑1 (MQ=255)                                                                                                            tCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGTCATAAATAGCGTGCGCATTTCTCTGTCGAGCCGCt                   <  2:1114857/139‑1 (MQ=255)                                                                                                                    tAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGa                                                                                                           <  2:291666/43‑1 (MQ=255)                                                                                                                    tAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGa                                                                                                           >  1:291666/1‑43 (MQ=255)                                                                                                                      aaCACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGTCATAAATAGCGTGCGCATTTCTCTGTCGAGCCGCTTCAGCAAGaa         <  1:1165609/139‑1 (MQ=255)                                                                                                                       acacCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGTCATAAATAGCGTGCGCATTtct                                >  2:786120/1‑115 (MQ=255)                                                                                                                       acacCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGTCATAAATAGCGTGCGCATTtct                                <  1:786120/115‑1 (MQ=255)                                                                                                                         acCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGTCATAAATAGCGTGCGCATTTCTCTGTCGAGCCGCTTCAGCAAGAAACg      >  2:905770/1‑139 (MQ=255)                                                                                                                           cTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCa                                                                          >  1:1364617/1‑69 (MQ=255)                                                                                                                           cTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCa                                                                          <  2:1364617/69‑1 (MQ=255)                                                                                                                           cTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGAt                                                                      >  2:669732/1‑73 (MQ=255)                                                                                                                           cTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGAt                                                                      <  1:669732/73‑1 (MQ=255)                                                                                                                            ttATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGTCATAAATAGCGTGCGCATTTCTCTGTCGAGCCGCTTCAGCAAGAAACGccc   >  1:623624/1‑139 (MQ=255)                                                                                                                             tATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGTCATAAATAGCGTGCGCATTTCTCTGTCGAGCCGCTTCAGCAAGAAACGcccc  >  2:259654/1‑139 (MQ=255)                                                                                                                                |                                                                                                                                        ATCGCTAGTTGAAGCACATTGCCGGATGCGACGCTTGCCACGTCTTATCCGGCCTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGTCATAAATAGCGTGCGCATTTCTCTGTCGAGCCGCTTCAGCAAGAAACGCCCC  >  NC_000913/2617948‑2618210

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