Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,796,034	A→G	R637R (CGT→CGC)	pheT ←	phenylalanine tRNA synthetase, beta subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,796,034	0	A	G	100.0%	217.7 / NA	64	R637R (CGT→CGC)	pheT	phenylalanine tRNA synthetase, beta subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (35/29);  total (35/29)

TCCAGTTCGAACACCAGAGTGCGACCGTTAAGATCCAGTTTACGTTCCAGTTCAGGATGAACAACCCCAACAAAACCAATACGTTCACCTTTCAGATAAATCGCTGCGGATTGCCCCGGATGCAGTGCCGGATTCGCTTCTGCACGGAACTCAACCTCATTCAGTTTACCGGTCAGGTCGAGAACGGATTCAAGATCGCCTTTCAAATCATAGAAATCAACGGTCTCTTTTGCCAGGTTCCAGTGCTCTTCGTAACGGTTACCGCAAATCACACCGGCTAACATC  >  NC_000913/1795891‑1796175                                                                                                                                                |                                                                                                                                              tCCAGTTCGAACACCAGAGTGCGACCGTTAAGATCCAGTTTACGTTCCAGTTCAGGATGAACAACCCCAACAAAACCAATACGTTCACCTTTCAGATAAATCGCTGCGGATTGCCCCGGATGCAGTGCCGGATTCGCTTCTGCGCGGAACt                                                                                                                                        >  2:850895/1‑151 (MQ=255)      ttCGAACACCAGAGTGCGACCGTTAAGATCCAGTTTACGTTCCAGTTCAGGATGAACAACCCCAACAAAACCAATACGTTCACCTTTCAGATAAATCGCTGCGGATTGCCCCGGATGCAGTGCCGGATTCGCTTCTGCGCGGAACTCAAcc                                                                                                                                   >  1:444386/1‑151 (MQ=255)           acacCAGAGTGCGACCGTTAAGATCCAGTTTACGTTCCAGTTCAGGATGAACAACCCCAACAAAACCAATACGTTCACCTTTCAGATAAATCGCTGCGGATTGCCCCGGATGCAGTGCCGGATTCGCTTCTGCGCGGAACTCAACCTCAtt                                                                                                                              >  1:847610/1‑151 (MQ=255)           acacCAGAGTGCGACCGTTAAGATCCAGTTTACGTTCCAGTTCAGGATGAACAACCCCAACAAAACCAATACGTTCACCTTTCAGATAAATCGCTGCGGATTGCCCCGGATGCAGTGCCGGATTCGCTTCTGCGCGGAACTCAACCTCAtt                                                                                                                              <  1:611396/151‑1 (MQ=255)           acacCAGAGTGCGACCGTTAAGATCCAGTTTACGTTCCAGTTCAGGATGAACAACCCCAACAAAACCAATACGTTCACCTTTCAGATAAATCGCTGCGGATTGCCCCGGATGCAGTGCCGGATTCGCTTCTGCGCGGAACTCAACCTCAtt                                                                                                                              >  2:775275/1‑151 (MQ=255)            cacCAGAGTGCGACCGTTAAGATCCAGTTTACGTTCCAGTTCAGGATGAACAACCCCAACAAAACCAATACGTTCACCTTTCAGATAAATCGCTGCGGATTGCCCCGGATGCAGTGCCGGATTCGCTTCTGCGCGGAACTCAACCTCATTc                                                                                                                             >  2:537445/1‑151 (MQ=255)              ccAGAGTGCGACCGTTAAGATCCAGTTTACGTTCCAGTTCAGGATGAACAACCCCAACAAAACCAATACGTTCACCTTTCAGATAAATCGCTGCGGATTGCCCCGGATGCAGTGCCGGATTCGCTTCTGCGCGGAACTCAACCTCATTc                                                                                                                             >  1:636464/1‑149 (MQ=255)              ccAGAGTGCGACCGTTAAGATCCAGTTTACGTTCCAGTTCAGGATGAACAACCCCAACAAAACCAATACGTTCACCTTTCAGATAAATCGCTGCGGATTGCCCCGGATGCAGTGCCGGATTCGCTTCTGCGCGGAACTCAACCTCATTCAg                                                                                                                           <  2:678974/151‑1 (MQ=255)                 gagTGCGACCGTTAAGATCCAGTTTACGTTCCAGTTCAGGATGAACAACCCCAACAAAACCAATACGTTCACCTTTCAGATAAATCGCTGCGGATTGCCCCGGATGCAGTGCCGGATTCGCTTCTGCGCGGAACTCAACCTCATTCAGttt                                                                                                                        <  1:816442/151‑1 (MQ=255)                 gagTGCGACCGTTAAGATCCAGTTTACGTTCCAGTTCAGGATGAACAACCCCAACAAAACCAATACGTTCACCTTTCAGATAAATCGCTGCGGATTGCCCCGGATGCAGTGCCGGATTCGCTTCTGCGCGGAACTCAACCTCATTCAGttt                                                                                                                        <  2:480788/151‑1 (MQ=255)                        aCCGTTAAGATCCAGTTTACGTTCCAGTTCAGGATGAACAACCCCAACAAAACCAATACGTTCACCTTTCAGATAAATCGCTGCGGATTGCCCCGGATGCAGTGCCGGATTCGCTTCTGCGCGGAACTCAACCTCATTCAGTTTACCGGTc                                                                                                                 >  2:658661/1‑151 (MQ=255)                          cGTTAAGATCCAGTTTACGTTCCAGTTCAGGATGAACAACCCCAACAAAACCAATACGTTCACCTTTCAGATAAATCGCTGCGGATTGCCCCGGATGCAGTGCCGGATTCGCTTCTGCGCGGAACTCAACCTCATTCAGTTTACCGGTCAg                                                                                                               <  2:199693/151‑1 (MQ=255)                           gTTAAGATCCAGTTTACGTTCCAGTTCAGGATGAACAACCCCAACAAAACCAATACGTTCACCTTTCAGATAAATCGCTGCGGATTGCCCCGGATGCATTGCCGGATTCACTTCTGCGCGGAACTCAACCTCATTCAGTTTACCGGTCAgg                                                                                                              <  2:186523/151‑1 (MQ=255)                              aaGATCCAGTTTACGTTCCAGTTCAGGATGAACAACCCCAACAAAACCAATACGTTCACCTTTCAGATAAATCGCTGCGGATTGCCCCGGATGCAGTGCCGGATTCGCTTCTGCGCGGAACTCAACCTCATTCAGTTTACCGGTCAGGTc                                                                                                            <  1:401599/150‑1 (MQ=255)                                 aTCCAGTTTACGTTCCAGTTCAGGATGAACAACCCCAACAAAACCAATACGTTCACCTTTCAGATAAATCGCTGCGGATTGCCCCGGATGCAGTGCCGGATTCGCTTCTGCGCGGAACTCAACCTCATTCAGTTTACCGGTCAGGTCg                                                                                                           >  2:802206/1‑148 (MQ=255)                                  tCCAGTTTACGTTCCAGTTCAGGATGAACAACCCCAACAAAACCAATACGTTCACCTTTCAGATAAATCGCTGCGGATTGCCCCGGATGCAGTGCCGGATTCGCTTCTGCGCGGAACTCAACCTCATTCAGTTTACCGGTCAGGTCgaga                                                                                                        >  2:420958/1‑150 (MQ=255)                                         tACGTTCCAGTTCAGGATGAACAACCCCAACAAAACCAATACGTTCACCTTTCAGATAAATCGCTGCGGATTGCCCCGGATGCAGTGCCGGATTCGCTTCTGCGCGGAACTCAACCTCATTCAGTTTACCGGTCAGGTCGAGAACGGATTc                                                                                                >  1:121719/1‑151 (MQ=255)                                         tACGTTCCAGTTCAGGATGAACAACCCCAACAAAACCAATACGTTCACCTTTCAGATAAATCGCTGCGGATTGCCCCGGATGCAGTGCCGGATTCGCTTCTGCGCGGAACTCAACCTCATTCAGTTTACCGGTCAGGTCGAGAACGGATTc                                                                                                >  2:419500/1‑151 (MQ=255)                                          aCGTTCCAGTTCAGGATGAACAACCCCAACAAAACCAATACGTTCACCTTTCAGATAAATCGCTGCGGATTGCCCCGGATGCAGTGCCGGATTCGCTTCTGCGCGGAACTCAACCTCATTCAGTTTACCGGTCAGGTCGAGAACGGATTc                                                                                                <  2:191180/150‑1 (MQ=255)                                          aCGTTCCAGTTCAGGATGAACAACCCCAACAAAACCAATACGTTCACCTTTCAGATAAATCGCTGCGGATTGCCCCGGATGCAGTGCCGGATTCGCTTCTGCGCGGAACTCAACCTCATTCAGTTTACCGGTCAGGTCGAGAACGGATTc                                                                                                <  2:856017/150‑1 (MQ=255)                                           cGTTCCAGTTCAGGATGAACAACCCCAACAAAACCAATACGTTCACCTTTCAGATAAAGCGCTGCGGATTGCCCCGGATGCAGTGCCGGATTCGCTTCTGCGCGGAACTCAACCTCATTCAGTTTACCGGTCAGGTCGAGAACGGATTCaa                                                                                              <  1:676181/151‑1 (MQ=255)                                            gTTCCAGTTCAGGATGAACAACCCCAACAAAACCAATACGTTCACCTTTCAGATAAATCGCTGCGGATTGCCCCGGATGCAGTGCCGGATTCGCTTCTGCGCGGAACTCAACCTCATTCAGTTTACCGGTCAGGTCGAGAACGGATTCa                                                                                               >  1:652745/1‑149 (MQ=255)                                                     cAGGATGAACAACCCCAACAAAACCAATACGTTCACCTTTCAGATAAATCGCTGCGGATTGCCCCGGATGCAGTGCCGGATTCGCTTCTGCGCGGAACTCAACCTCATTCAGTTTACCGGTCAGGTCGAGAACGGATTCAAGATCGCCttt                                                                                    <  1:529337/151‑1 (MQ=255)                                                        gATGAACAACCCCAACAAAACCAATACGTTCACCTTTCAGATAAATCGCTGCGGATTGCCCCGGATGCAGTGCCGGATTCGCTTCTGCGCGGAACTCAACCTCATTCAGTTTACCGGTCAGGTCGAGAACGGATTCAAGATCGCCTTTCaa                                                                                 <  1:17004/151‑1 (MQ=255)                                                         aTGAACAACCCCAACAAAACCAATACGTTCACCTTTCAGATAAATCGCTGCGGATTGCCCCGGATGCAGTGCCGGATTCGCTTCTGCGCGGAACTCAACCTCATTCAGTTTACCGGTCAGGTCGAGAACGGATTCAAGATCGCCTTTCaaa                                                                                <  2:700483/151‑1 (MQ=255)                                                         aTGAACAACCCCAACAAAACCAATACGTTCACCTTTCAGATAAATCGCTGCGGATTGCCCCGGATGCAGTGCCGGATTCGCTTCTGCGCGGAACTCAACCTCATTCAGTTTACCGGTCAGGTCGAGAACGGATTCAAGATCGCCTTTCaa                                                                                 >  1:826961/1‑150 (MQ=255)                                                           gAACAACCCCAACAAAACCAATACGTTCACCTTTCAGATAAATCGCTGCGGATTGCCCCGGATGCAGTGCCGGATTCGCTTCTGCGCGGAACTCAACCTCATTCAGTTTACCGGTCAGGTCGAGAACGGATTCAAGATCGCCTTTCAAAt                                                                               >  1:93248/1‑150 (MQ=255)                                                           gAACAACCCCAACAAAACCAATACGTTCACCTTTCAGATAAATCGCTGCGGATTGCCCCGGATGCAGTGCCGGATTCGCTTCTGCGCGGAACTCAACCTCATTCAGTTTACCGGTCAGGTCGAGAACGGATTCAAGATCGCCTTTCAAAt                                                                               <  1:151247/150‑1 (MQ=255)                                                              caacCCCAACGAAACTAATACGTTCACCTTTCAGATAAATCGCGGCGGATTGCCCCGGATGCAGTGCCGGATTCGCCTCTACGCGGAACTCGACCTCCTTCAGTTTACCGGTCCGGTCGAGAACGGATTCAAGATCGCCTTTCAAATc                                                                              >  2:94073/1‑148 (MQ=255)                                                                acCCCAACAAAACCAATACGTTCACCTTTCAGATAAATCGCTGCGGATTGCCCCGGATGCAGTGCCGGATTCGCTTCTGCGCGGAACTCAACCTCATTCAGTTTACCGGTCAGGTCGAGAACGGATTCAAGATCGCCTTTCAAATCAt                                                                            >  1:303279/1‑148 (MQ=255)                                                                 ccccAACAAAACCAATACGTTCACCTTTCAGATAAATCGCTGCGGATTGCCCCGGATGCAGTGCCGGATTCGCTTCTGCGCGGAACTCAACCTCATTCAGTTTACCGGTCAGGTCGAGAACGGATTCAAGATCGCCTTTCAAATCATAGa                                                                         >  2:236906/1‑150 (MQ=255)                                                                   ccAACAAAACCAATACGTTCACCTTTCAGATAAATCGCTGCGGATTGCCCCGGATGCAGTGCCGGATTCGCTTCTGCGCGGAACTCAACCTCATTCAGTTTACCGGTCAGGTCGAGAACGGATTCAAGATCGCCTTTCAAATCATAGAAAt                                                                      >  1:312143/1‑151 (MQ=255)                                                                          aaCCAATACGTTCACCTTTCAGATAAATCGCTGCGGATTGCCCCGGATGCAGTGCCGGATTCGCTTCTGCGCGGAACTCAACCTCATTCAGTTTACCGGTCAGGTCGAGAACGGATTCAAGATCGCCTTTCAAATCATAGAAATCAACGGt                                                               <  1:48708/151‑1 (MQ=255)                                                                          aaCCAATACGTTCACCTTTCAGATAAATCGCTGCGGATTGCCCCGGATGCAGTGCCGGATTCGCTTCTGCGCGGAACTCAACCTCATTCAGTTTACCGGTCAGGTCGAGAACGGATTCAAGATCGCCTTTCAAATCATAGAAATCAACGGt                                                               <  2:324978/151‑1 (MQ=255)                                                                          aaCCAATACGTTCACCTTTCAGATAAATCGCTGCGGATTGCCCCGGATGCAGTGCCGGATTCGCTTCTGCGCGGAACTCAACCTCATTCAGTTTACCGGTCAGGTCGAGAACGGATTCAAGATCGCCTTTCAAATCATAGAAATCAACGGt                                                               <  2:313487/151‑1 (MQ=255)                                                                             cAATACGTTCACCTTTCAGATAAATCGCTGCGGATTGCCCCGGATGCAGTGCCGGATTCGCTTCTGCGCGGAACTCAACCTCATTCAGTTTACCGGTCAGGTCGAGAACGGATTCAAGATCGCCTTTCAAATCATAGAAATCAACGGtctc                                                            >  1:751153/1‑151 (MQ=255)                                                                             cAATACGTTCACCTTTCAGATAAATCGCTGCGGATTGCCCCGGATGCAGTGCCGGATTCGCTTCTGCGCGGAACTCAACCTCATTCAGTTTACCGGTCAGGTCGAGAACGGATTCAAGATCGCCTTTCAAATCATAGAAATCAACGGtctc                                                            >  2:713558/1‑151 (MQ=255)                                                                              aaTACGTTCACCTTTCAGATAAATCGCTGCGGATTGCCCCGGATGCAGTGCCGGATTCGCTTCTGCGCTGAACTCAACCTCATTCAGTTTACCTGTCAGGTCGAGAACGGCTTCAAGATCGCCTTTCCAATCATAGAAATCAACGGTctct                                                           >  2:573082/1‑151 (MQ=255)                                                                                tACGTTCACCTTTCAGATAAATCGCTGCGGATTGCCCCGGATGCAGTGCCGGATTCGCTTCTGCGCGGAACTCAACCTCATTCAGTTTACCGGTCAGGTCGAGAACGGATTCAAGATCGCCTTTCAAATCATAGAAATCAACGGTCTCttt                                                         >  1:668176/1‑151 (MQ=255)                                                                                 aCGTTCACCTTTCAGATAAATCGCTGCGGATTGCCCCGGATGCAGTGCCGGATTCGCTTCTGCGCGGAACTCAACCTCATTCAGTTTACCGGTCAGGTCGAGAACGGATTCAAGATCGCCTTTCAAATCATAGAAATCAACGGTCTCtttt                                                        <  1:741735/151‑1 (MQ=255)                                                                                     tCACCTTTCAGATAAATCGCTGCGGATTGCCCCGGATGCAGTGCCGGATTCGCTTCTGCGCGGAACTCAACCTCATTCAGTTTACCGGTCAGGTCGAGAACGGATTCAAGATCGCCTTTCAAATCATAGAAATCAACGGTCTCTTTTGcc                                                     >  1:312731/1‑150 (MQ=255)                                                                                        ccTTTCAGATAAATCGCTGCGGATTGCCCCGGATGCAGTGCCGGATTCGCTTCTGCGCGGAACTCAACCTCATTCAGTTTACCGGTCAGGTCGAGAACGGATTCAAGATCGCCTTTCAAATCATAGAAATCAACGGTCTCTTTTGCCAGGt                                                 <  1:866110/151‑1 (MQ=255)                                                                                          tttCAGATAAATCGCTGCGGATTGCCCCGGATGCAGTGCCGGATTCGCTTCTGCGCGGAACTCAACCTCATTCAGTTTACCGGTCAGGTCGAGAACGGATTCAAGATCGCCTTTCAAATCATAGAAATCAACGGTCTCTTTTGCCAGGTTc                                               >  1:37412/1‑151 (MQ=255)                                                                                           ttCAGATAAATCGCTGCGGATTGCCCCGGATGCAGTGCCGGATTCGCTTCTGCGCGGAACTCAACCTCATTCAGTTTACCGGTCAGGTCGAGAACGGATTCAAGATCGCCTTTCAAATCATAGAAATCAACGGTCTCTTTTGCCAGGtt                                                <  2:45700/149‑1 (MQ=255)                                                                                           ttCAGATAAATCGCTGCGGATTGCCCCGGATGCAGTGCCGGATTCGCTTCTGCGCGGAACTCAACCTCATTCAGTTTACCGGTCAGGTCGAGAACGGATTCAAGATCGCCTTTCAAATCATAGAAATCAACGGTCTCTTTTGCCAGGtt                                                <  1:419545/149‑1 (MQ=255)                                                                                           ttCAGATAAATCGCTGCGGATTGCCCCGGATGCAGTGCCGGATTCGCTTCTGCGCGGAACTCAACCTCATTCAGTTTACCGGTCAGGTCGAGAACGGATTCAAGATCGCCTTTCAAATCATAGAAATCAACGGTCTCTTTTGCCAGGTTc                                               <  2:665854/150‑1 (MQ=255)                                                                                             cAGATAAATCGCTGCGGATTGCCCCGGATGCAGTGCCGGATTCGCTTCTGCGCGGAACTCAACCTCATTCAGTTTACCGGTCAGGTCGAGAACGGATTCAAGATCGCCTTTCAAATCATAGAAATCAACGGTCTCTTTTGCCAGGTTCCAg                                            <  1:332716/151‑1 (MQ=255)                                                                                                  aaaTCGCTGCGGATTGCCCCGGATGCAGTGCCGGATTCGCTTCTGCGCGGAACTCAACCTCATTCAGTTTACCGGTCAGGTCGAGAACGGATTCAAGATCGCCTTTCAAATCATAGAAATCAACGGTCTCTTTTGCCAGGTTCCAGTGctc                                       >  1:428183/1‑151 (MQ=255)                                                                                                  aaaTCGCTGCGGATTGCCCCGGATGCAGTGCCGGATTCGCTTCTGCGCGGAACTCAACCTCATTCAGTTTACCGGTCAGGTCGAGAACGGATTCAAGATCGCCTTTCAAATCATAGAAATCAACGGTCTCTTTTGCCAGGTTCCAGTGct                                        <  2:317951/150‑1 (MQ=255)                                                                                                       gCTGCGGATTGCCCCGGATGCAGTGCCGGATTCGCTTCTGCGCGGAACTCAACCTCATTCAGTTTACCGGTCAGGTCGAGAACGGATTCAAGATCGCCTTTCAAATCATAGAAATCAACGGTCTCTTTTGCCAGGTTCCAGTGCTCTTCGt                                  >  1:590425/1‑151 (MQ=255)                                                                                                       gCTGCGGATTGCCCCGGATGCAGTGCCGGATTCGCTTCTGCGCGGAACTCAACCTCATTCAGTTTACCGGTCAGGTCGAGAACGGATTCAAGATCGCCTTTCAAATCATAGAAATCAACGGTCTCTTTTGCCAGGTTCCAGTGCTCTTCGt                                  <  1:885084/151‑1 (MQ=255)                                                                                                            ggATTGCCCCGGATGCAGTGCCGGATTCGCTTCTGCGCGGAACTCAACCTCATTCAGTTTACCGGTCAGGTCGAGAACGGATTCAAGATCGCCTTTCAAATCATAGAAATCAACGGTCTCTTTTGCCAGGTTCCAGTGCTCTTCGTAACgg                             <  2:368040/151‑1 (MQ=255)                                                                                                            ggATTGCCCCGGATGCAGTGCCGGATTCGCTTCTGCGCGGAACTCAACCTCATTCAGTTTACCGGTCAGGTCGAGAACGGATTCAAGATCGCCTTTCAAATCATAGAAATCAACGGTCTCTTTTGCCAGGTTCCAGTGCTCTTCGTAACgg                             <  2:272590/151‑1 (MQ=255)                                                                                                                 gCCCCGGATGCAGTGCCGGGTTCGCTTCTGCGCGGAACTCAACCTCATTCAGTTTACCGGTCAGGTCGAGAACGGATTCAAGATCGCCTTTCAAATCATAGAAATCAACGGTCTCTTTTGCCAGGTTCCAGTGCTCTTCGTAACGGTTAcc                        >  2:244555/1‑151 (MQ=255)                                                                                                                 gCCCCGGATGCAGTGCCGGATTCGCTTCTGCGCGGAACTCAACCTCATTCAGTTTACCGGTCAGGTCGAGAACGGATTCAAGATCGCCTTTCAAATCATAGAAATCAACGGTCTCTTTTGCCAGGTTCCAGTGCTCTTCGTAACGGTTAcc                        >  1:788431/1‑151 (MQ=255)                                                                                                                  ccccGGATGCAGTGCCGGATTCGCTTCTGCGCGGAACTCAACCTCATTCAGTTTACCGGTCAGGTCGAGAACGGATTCAAGATCGCCTTTCAAATCATAGAAATCAACGGTCTCTTTTGCCAGGTTCCAGTGCTCTTCGTAACGGTTACCg                       >  2:614627/1‑151 (MQ=255)                                                                                                                  ccccGGATGCAGTGCCGGATTCGCTTCTGCGCGGAACTCAACCTCATTCAGTTTACCGGTCAGGTCGAGAACGGATTCAAGATCGCCTTTCAAATCATAGAAATCAACGGTCTCTTTTGCCAGGTTCCAGTGCTCTTCGTAACGGTTACCg                       >  2:119927/1‑151 (MQ=255)                                                                                                                        aTGCAGTGCCGGATTCGCTTCTGCGCGGAACTCAACCTCATTCAGTTTACCGGTCAGGTCGAGAACGGATTCAAGATCGCCTTTCAAATCATAGAAATCAACGGTCTCTTTTGCCAGGTTCCAGTGCTCTTCGTAACGGTTACCGCAAATc                 >  2:850222/1‑151 (MQ=255)                                                                                                                               gCCGGATTCGCTTCTGCGCGGAACTCAACCTCATTCAGTTTACCGGTCAGGTCGAGAACGGATTCAAGATCGCCTTTCAAATCATAGAAATCAACGGTCTCTTTTGCCAGGTTCCAGTGCTCTTCGTAACGGTTACCGCAAATCACACCgg          <  1:433945/151‑1 (MQ=255)                                                                                                                                ccGGATTCGCTTCTGCGCGGAACTCAACCTCATTCAGTTTACCGGTCAGGTCGAGAACGGATTCAAGATCGCCTTTCAAATCATAGAAATCAACGGTCTCTTTTGCCAGGTTCCAGTGCTCTTCGTAACGGTTACCGCAAATCACACCGGc         >  2:605966/1‑151 (MQ=255)                                                                                                                                 cGGATTCGCTTCTGCGCGGAACTCAACCTCATTCAGTTTACCGGTCAGGTCGAGAACGGATTCAAGATCGCCTTTCAAATCATAGAAATCAACGGTCTCTTTTGCCAGGTTCCAGTGTTCTTCGTAACGGTTACGGCAAATCACACCg           >  1:20804/1‑148 (MQ=255)                                                                                                                                 cGGATTCGCTTCTGCGCGGAACTCAACCTCATTCAGTTTACCGGTCAGGTCGAGAACGGATTCAAGATCGCCTTTCAAATCATAGAAATCAACGGTCTCTTTTGCCAGGTTCCAGTGCTCTTCGTAACGGTTACCGCAAATCACACCGGCt        <  1:506916/151‑1 (MQ=255)                                                                                                                                      tCGCTTCTGCGCGGAACTCAACCTCATTCAGTTTACCGGTCAGGTCGAGAACGGATTCAAGATCGCCTTTCAAATCATAGAAATCAACGGTCTCTTTTGCCAGGTTCCAGTGCTCTTCGTAACGGTTACCGCAAATCACACCGGCTAACAt   >  1:490134/1‑151 (MQ=255)                                                                                                                                        gCTTCTGCGCGGAACTCAACCTCATTCAGTTTACCGGTCAGGTCGAGAACGGATTCAAGATCGCCTTTCAAATCATAGAAATCAACGGTCTCTTTTGCCAGGTTCCAGTGCTCTTCGTAACGGTTACCGCAAATCACACCGGCTAACATc  >  1:275580/1‑150 (MQ=255)                                                                                                                                                |                                                                                                                                              TCCAGTTCGAACACCAGAGTGCGACCGTTAAGATCCAGTTTACGTTCCAGTTCAGGATGAACAACCCCAACAAAACCAATACGTTCACCTTTCAGATAAATCGCTGCGGATTGCCCCGGATGCAGTGCCGGATTCGCTTCTGCACGGAACTCAACCTCATTCAGTTTACCGGTCAGGTCGAGAACGGATTCAAGATCGCCTTTCAAATCATAGAAATCAACGGTCTCTTTTGCCAGGTTCCAGTGCTCTTCGTAACGGTTACCGCAAATCACACCGGCTAACATC  >  NC_000913/1795891‑1796175

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