Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,798,559	T→C	A128A (GCA→GCG)	pheS ←	phenylalanine tRNA synthetase, alpha subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,798,559	0	T	C	100.0%	230.4 / NA	66	A128A (GCA→GCG)	pheS	phenylalanine tRNA synthetase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (39/27);  total (39/27)

GTGCGGATCTGTACGCCAGAGGTCTGGGTACGCAGCAGGCGGGTAGTGTCAAACCAGAAAGTGTCGTGGTCAGCGCGCGCCGGGTGGTGACCAGGAATGTTCAGAGCATCGAAGTTATGATAATCGTCTTCGATTTCCGGCCCGGTTGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCGGATGCAGACCGCCGTTTTCAATGCGACGACCTGGCAGAGAGACATCAATCGTTTCCGCCGCCAGACGCGCATTCAGTGCAGCG  >  NC_000913/1798413‑1798703                                                                                                                                                   |                                                                                                                                                 gTGCGGATCTGTACGCCAGAGGTCTGGGTACGCAGCAGGCGGGTAGTGTCAAACCAGAAAGTGTCGTGGTCAGCGCGCGCCGGGGGGTGACCAGGAATGTTCAGAGCATCGAAGTTATGATAAACGTCTTCGATTTCCGGCCCGGTCGcc                                                                                                                                               >  1:1252637/1‑150 (MQ=255)     ggATCTGTACGCCAGAGGTCTGGGTACGCAGCAGGCGGGTAGTGTCAAACCAGAAAGTGTCGTGGTCAGCGCGCGCCGGGTGGTGACCAGGAATGTTCAGAGCATCGAAGTTATGATAATCGTCTTCGATTTCCGGCCCGGTCGCCACgg                                                                                                                                           >  2:1266593/1‑150 (MQ=255)       aTCTGTACGCCAGAGGTCTGGGTACGCAGCAGGCGGGTAGTGTCAAACCAGAAAGTGTCGTGGTCAGCGCGCGCCGGGTGGTGACCAGGAATGTTCAGAGCATCGAAGTTATGATAATCGTCTTCGATTTCCGGCCCGGTCGCCACGGTa                                                                                                                                         >  2:673038/1‑150 (MQ=255)       aTCTGTACGCCAGAGGTCTGGGTACGCAGCAGGCGGGTAGTGTCAAACCAGAAAGGGTCGTGGTCAGCGCGCGCCGGGTGGGGACCCGGAAAGTTCAGAGCATCGAAGTTATGATAATCGTCTTCGATTTTCGGCCCGGTCGCCACGGTa                                                                                                                                         >  1:50787/1‑150 (MQ=255)        tCTGTACGCCAGAGGTCTGGGTACGCAGCAGGCGGGTAGTGTCAAACCAGAAAGTGTCGTGGTCAGCGCGCGCCGGGGGGTGACCAGGAATGTTCAGAGCATCGAAGTTATGATAATCGTCTTCGATTTCCGGCCCGGTCGCCACGGTaa                                                                                                                                        >  1:391182/1‑150 (MQ=255)          tGTACGCCAGAGGTCTGGGTACGCAGCAGGCGGGTAGTGTCAAACCAGAAAGTGTCGTGGTCAGCGCGCGCCGGGTGGTGACCAGGAATGTTCAGAGCATCGAAGTTATGATAATCGTCTTCGATTTCCGGCCCGGTCGCCACGGTAAAg                                                                                                                                      >  2:586238/1‑150 (MQ=255)           gTACGCCAGAGGTCTGGGTACGCAGCAGGCGGGTAGTGTCAAACCAGAAAGTGTCGTGGTCAGCGCGCGCCGGGTGGTGACCAGGAATGTTCAGAGCATCGAAGTTATGATAATCGTCTTCGATTTCCGGCCCGGTCGCCACGGTAAAGcc                                                                                                                                    <  1:914487/151‑1 (MQ=255)           gTACGCCAGAGGTCTGGGTACGCAGCAGGCGGGTAGTGTCAAACCAGAAAGTGTCGTGGTCAGCGCGCGCCGGGTGGTGACCAGGAATGTTCAGAGCATCGAAGTTATGATAATCGTCTTCGATTTCCGGCCCGGTCGCCACGGTAAAGcc                                                                                                                                    <  1:746226/151‑1 (MQ=255)                      gTCTGGGTACGCAGCAGGCGGGTAGTGTCAAACCAGAAAGTGTCGTGGTCAGCGCGCGCCGGGTGGTGACCAGGAATGTTCAGAGCATCGAAGTTATGATAATCGTCTTCGATTTCCGGCCCGGTCGCCACGGTAAAGCCAAGCTCACCga                                                                                                                         <  2:479146/151‑1 (MQ=255)                       tCTGGGTACGCAGCAGGCGGGTAGTGTCAAACCAGAAAGTGTCGTGGTCAGCGCGCGCCGGGTGGTGACCAGGAATGTTCAGAGCATCGAAGTTATGATAATCGTCTTCGATTTCCGGCCCGGTCGCCACGGTAAAGCCAAGCTCACCga                                                                                                                         >  2:458420/1‑150 (MQ=255)                         tGGGTACGCAGCAGGCGGGTAGTGTCAAAACAGAAAGTGGCGTGGTCAGCGCGCGCCGGGCGGTGACCAGGAATGTTCAGAGCATCGAAGTTATGATAATCGTCTTCGATTTCCGGCCCGGTCGCCACGGTAAAGCCAAGCTCAACga                                                                                                                         >  1:1026539/1‑148 (MQ=255)                          gggTACGCAGCAGGCGGGTAGTGTCAAACCAGAAAGTGTCGTGGTCAGCGCGCGCCGGGTGGTGACCAGGAATGTTCAGAGCATCGAAGTTATGATAATCGTCTTCGATTTCCGGCCCGGTCGCCACGGTAAAGCCAAGCTCACCgaagaa                                                                                                                     <  1:975057/151‑1 (MQ=255)                            gTACGCAGCAGGCGGGTAGTGTCAAACCAGAAAGTGTCGTGGTCAGCGCGCGCCGGGTGGTGACCAGGAATGTTCAGAGCATCGAAGTTATGATAATCGTCTTCGATTTCCGGCCCGGTCGCCACGGTAAAGCCAAGCTCACCGAAGAAAc                                                                                                                   >  2:1142184/1‑151 (MQ=255)                               cGCAGCAGGCGGGTAGTGTCAAACCAGAAAGTGTCGTGGTCAGCGCGCGCCGGGTGGTGACCAGGAATGTTCAGAGCATCGAAGTTATGATAATCGTCTTCGATTTCCGGCCCGGTCGCCACGGTAAAGCCAAGCTCACCGAAGAAACttt                                                                                                                <  1:561858/151‑1 (MQ=255)                               cGCAGCAGGCGGGTAGTGTCAAACCAGAAAGTGTCGTGGTCAGCGCGCGCCGGGTGGTGACCAGGAATGTTCAGAGCATCGAAGTTATGATAATCGTCTTCGATTTCCGGCCCGGTCGCCACGGTAAAGCCAAGCTCACCGAAGAAACttt                                                                                                                <  2:525184/151‑1 (MQ=255)                                gcagcaGGCGGGTAGTGTCAAACCAGAAAGTGTCGTGGTCAGCGCGCGCCGGGTGGTGACCAGGAATGTTCAGAGCATCGAAGTTATGATAATCGTCTTCGATTTCCGGCCCGGTCGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTc                                                                                                               <  2:531249/151‑1 (MQ=255)                                         gggTAGTGTCAAACCAGAAAGTGTCGTGGTCAGCGCGCGCCGGGTGGTGACCAGGAATGTTCAGAGCATCGAAGTTATGATAATCGTCTTCGATTTCCGGCCCGGTCGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTc                                                                                                      <  1:113289/151‑1 (MQ=255)                                             aGTGTCAAACCAGAAAGTGTCGTGGTCAGCGCGCGCCGGGTGGTGACCAGGAATGTTCAGAGCATCGAAGTTATGATAATCGTCTTCGATTTCCGGCCCGGTCGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATg                                                                                                  >  1:1291355/1‑151 (MQ=255)                                             aGTGTCAAACCAGAAAGTGTCGTGGTCAGCGCGCGCCGGGTGGTGACCAGGAATGTTCAGAGCATCGAAGTTATGATAATCGTCTTCGATTTCCGGCCCGGTCGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATg                                                                                                  >  1:917235/1‑151 (MQ=255)                                                  cAAACCAGAAAGTGTCGTGGTCAGCGCGCGCCGGGTGGTGACCAGGAATGTTCAGAGCATCGAAGTTATGATAATCGTCTTCGATTTCCGGCCCGGTCGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACg                                                                                             >  2:120276/1‑151 (MQ=255)                                                      ccAGAAAGTGTCGTGGTCAGCGCGCGCCGGGTGGTGACCAGGAATGTTCAGAGCATCGAAGTTATGATAATCGTCTTCGATTTCCGGCCCGGTCGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTa                                                                                         <  1:603058/151‑1 (MQ=255)                                                       cAGAAAGTGTCGTGGTCAGCGCGCGCCGGGT‑GTGACCAGGAATGTTCAGAGCATCGAAGTTATGATAATCGTCTTCGATTTCCGGCCCGGTCGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAAc                                                                                       >  2:1120312/1‑151 (MQ=255)                                                          aaaGTGTCGTGGTCAGCGCGCGCCGGGTGGTGACCAGGAATGTTCAGAGCATCGAAGTTATGATAATCGTCTTCGATTTCCGGCCCGGTCGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAAcc                                                                                      >  1:22083/1‑150 (MQ=255)                                                          aaaGTGTCGTGGTCAGCGCGCGCCGGGTGGTGACCAGGAATGTTCAGAGCATCGAAGTTATGATAATCGTCTTCGATTTCCGGCCCGGTCGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCg                                                                                     <  1:835893/151‑1 (MQ=255)                                                          aaaGTGTCGTGGTCAGCGCGCGCCGGGTGGTGACCAGGAATGTTCAGAGCATCGAAGTTATGATAATCGTCTTCGATTTCCGGCCCGGTCGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCg                                                                                     >  2:1207092/1‑151 (MQ=255)                                                          aaaGTGTCGTGGTCAGCGCGCGCCGGGTGGTGACCAGGAATGTTCAGAGCATCGAAGTTATGATAATCGTCTTCGATTTCCGGCCCGGTCGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCg                                                                                     >  2:662342/1‑151 (MQ=255)                                                          aaaGTGTCGTGGTCAGCGCGCGCCGGGTGGTGACCAGGAATGTTCAGAGCATCGAAGTTATGATAATCGTCTTCGATTTCCGGCCCGGTCGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCg                                                                                     <  1:554695/151‑1 (MQ=255)                                                               gtCGTGGTCAGCGCGCGCCGGGTGGTGACCAGGAATGTTCAGAGCATCGAAGTTATGATAATCGTCTTCGATTTCCGGCCCGGTCGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCGGAt                                                                                  >  1:128409/1‑149 (MQ=255)                                                               gtCGTGGTCAGCGCGCGCCGGGTGGTGACCAGGAATGTTCAGAGCATCGAAGTTATGATAATCGTCTTCGATTTCCGGCCCGGTCGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCGGATGc                                                                                >  2:113632/1‑151 (MQ=255)                                                                 cGTGGTCAGCGCGCGCCGGGTGGTGACCAGGAATGTTCAGAGCATCGAAGTTATGATAATCGTCTTCGATTTCCGGCCCGGTCGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCGGATGCa                                                                               <  1:1234781/150‑1 (MQ=255)                                                                   tGGTCAGCGCGCGCCGGGTGGTGACCAGGAATGTTCAGAGCATCGAAGTTATGATAATCGTCTTCGATTTCCGGCCCGGTCGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCGGATGCAGAc                                                                            >  1:454040/1‑151 (MQ=255)                                                                    ggTCAGCGCGCGCCGGGTGGTGACCAGGAATGTTCAGAGCATCGAAGTTATGATAATCGTCTTCGATTTCCGGCCCGGTCGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCGGATGCAGAcc                                                                           <  1:292199/151‑1 (MQ=255)                                                                      tCAGCGCGCGCCGGGTGGTGACCAGGAATGTTCAGAGCATCGAAGTTATGATAATCGTCTTCGATTTCCGGCCCGGTCGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCGGATGCAGAccgc                                                                         >  2:460883/1‑151 (MQ=255)                                                                       cAGCGCGCGCCGGGTGGTGACCAGGAATGTTCAGAGCATCGAAGTTATGATAATCGTCTTCGATTTCCGGCCCGGTCGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCGGATGCAGAccgcc                                                                        >  1:197753/1‑151 (MQ=255)                                                                        aGCGCGCGCCGGGTGGTGACCAGGAATGTTCAGAGCATCGAAGTTATGATAATCGTCTTCGATTTCCGGCCCGGTCGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCGGATGCAGAccgccg                                                                       >  2:391090/1‑151 (MQ=255)                                                                         gcgcgcgcCGGGTGGTGACCAGGAATGTTCAGAGCATCGAAGTTATGATAATCGTCTTCGATTTCCGGCCCGGTCGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCGGATGCAGACCGCCGt                                                                      <  2:332612/151‑1 (MQ=255)                                                                            cgcgcCGGGTGGTGACCAGGAATGTTCAGAGCATCGAAGTTATGATAATCGTCTTCGATTTCCGGCCCGGTCGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCGGATGCAGACCGCCGtttt                                                                   <  2:407026/151‑1 (MQ=255)                                                                             gcgcCGGGTGGTGACCAGGAATGTTCAGAGCATCGAAGTTATGATAATCGTCTTCGATTTCCGGCCCGGTCGCCCCGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCGGATGCAGACCGCCGt                                                                      >  1:712859/1‑147 (MQ=255)                                                                             gcgcCGGGTGGTGACCAGGAATGTTCAGAGCATCGAAGTTATGATAATCGTCTTCGATTTCCGGCCCGGTCGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCGGATGCAGACCGCCGTTTTc                                                                  >  1:905190/1‑151 (MQ=255)                                                                                  gggTGGTGACCAGGAATGTTCAGAGCATCGAAGTTATGATAATCGTCTTCGATTTCCGGCCCGGTCGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCGGATGCAGACCGCCGTTTTCAATg                                                              >  1:1072764/1‑150 (MQ=255)                                                                                      ggtgACCAGGAATGTTCAGAGCATCGAAGTTATGATAATCGTCTTCGATTTCCGGCCGGGTCGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCGGATGCAGACCGCCGTTTTCAATGcgac                                                          <  1:81757/150‑1 (MQ=255)                                                                                           ccAGGAATGTTCAGAGCATCGAAGTTATGATAATCGTCTTCGATTTCCGGCCCGGTCGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCGGATGCAGACCGCCGTTTTCAATGCGACGACCt                                                     <  1:1082735/150‑1 (MQ=255)                                                                                           ccAGGAATGTTCAGAGCATCGAAGTTATGATAATCGTCTTCGATTTCCGGCCCGGTCGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCGGATGCAGACCGCCGTTTTCAATGCGACGACCTg                                                    >  2:1262355/1‑151 (MQ=255)                                                                                               gAATGTTCAGAGCATCGAAGTTATGATAATCGTCTTCGATTTCCGGCCCGGTCGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCGGATGCAGACCGCCGTTTTCAATGCGACGACCTgg                                                   >  1:747519/1‑148 (MQ=255)                                                                                                   gTTCAGAGCATCGAAGTTATGATAATCGTCTTCGATTTCCGGCCCGGTCGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCGGATGCAGACCGCCGTTTTCAATGCGACGACCTGGCagagag                                            <  2:277144/151‑1 (MQ=255)                                                                                                   gTTCAGAGCATCGAAGTTATGATAATCGTCTTCGATTTCCGGCCCGGTCGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCGGATGCAGACCGCCGTTTTCAATGCGACGACCTGGCagaga                                             <  1:792514/150‑1 (MQ=255)                                                                                                   gTTCAGAGCATCGAAGTTATGATAATCGTCTTCGATTTCCGGCCCGGTCGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCGGATGCAGACCGCCGTTTTCAATGCGACGACCTGGCagag                                              >  2:1143870/1‑149 (MQ=255)                                                                                                    ttCAGAGCATCGAAGTTATGATAATCGTCTTCGATTTCCGGCCCGGTCGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCGGATGCAGACCGCCGTTTTCAATGCGACGACCTGGCagag                                              >  2:903025/1‑148 (MQ=255)                                                                                                      cAGAGCATCGAAGTTATGATAATCGTCTTCGATTTCCGGCCCGGTCGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCGGATGCAGACCGCCGTTTTCAATGCGACGACCTGGCAGAGAGACa                                         <  1:167225/151‑1 (MQ=255)                                                                                                        gagCATCGAAGTTATGATAATCGTCTTCGATTTCCGGCCCGGTCGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCGGATGCAGACCGCCGTTTTCAATGCGACGACCTGGCAGAGAGACATc                                       <  2:1198125/151‑1 (MQ=255)                                                                                                         agCATCGAAGTTATGATAATCGTCTTCGATTTCCGGCCCGGTCGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCGGATGCAGACCGCCGTTTTCAATGCGACGACCTGGCAGAGAGACATc                                       >  2:501621/1‑150 (MQ=255)                                                                                                           cATCGAAGTTATGATAATCGTCTTCGATTTCCGGCCCGGTCGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCGGATGCAGACCGCCGTTTTCAATGCGACGACCTGGCAGAGAGACATCaa                                     <  1:1248909/150‑1 (MQ=255)                                                                                                           cATCGAAGTTATGATAATCGTCTTCGATTTCCGGCCCGGTCGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCGGATGCAGACCGCCGTTTTCAATGCGACGACCTGGCAGAGAGACATCAAt                                    <  1:830233/151‑1 (MQ=255)                                                                                                              cGAAGTTATGATAATCGTCTTCGATTTCCGGCCCGGTCGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCGGATGCAGACCGCCGTTTTCAATGCGACGACCTGGCAGAGAGACATCAATCGt                                 >  1:1061160/1‑151 (MQ=255)                                                                                                                aaGTTATGATAATCGTCTTCGATTTCCGGCCCGGTCGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACGGGATGCAGACCGCCGTTTTCAATGCGACGACCTGGCAGAGAGACATCAATCGttt                               <  1:211159/151‑1 (MQ=255)                                                                                                                aaGTTATGATAATCGTCTTCGATTTCCGGCCCGGTCGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCGGATGCAGACCGCCGTTTTCAATGCGACGACCTGGCAGAGAGACATCAATCGttt                               >  2:424855/1‑151 (MQ=255)                                                                                                                   ttATGATAATCGTCTTCGATTTCCGGCCCGGTCGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCGGATGCAGACCGCCGTTTTCAATGCGACGACCTGGCAGAGAGACATCAATCGTTTccg                            >  1:161423/1‑151 (MQ=255)                                                                                                                     aTGATAATCGTCTTCGATTTCCGGCCCGGTCGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCGGATGCAGACCGCCGTTTTCAATGCGACGACCTGGCAGAGAGACATCAATCGTTTccgcc                          >  1:34701/1‑151 (MQ=255)                                                                                                                      tGATAATCGTCTTCGATTTCCGGCCCGGTCGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCGGATGCAGACCGCCGTTTTCAATGCGACGACCTGGCAGAGAGACATCAATCGTTTccgccg                         >  2:842507/1‑151 (MQ=255)                                                                                                                           aTCGTCTTCGATTTCCGGCCCGGTCGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTGCGATACGGTCGATGGTACGGGTAACCGGATGCAGACCGCCGTTTTCAATGCGACGACCTGGCAGAGAGACATCAATCGTTTCcgccgcaga                     <  2:720930/150‑4 (MQ=255)                                                                                                                               tCTTCGATTTCCGGCCCGGTCGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCGGATGCAGACCGCCGTTTTCAATGCGACGACCTGGCAGAGAGACATCAATCGTTTCCGCCGCCAGAcgcg                >  1:41311/1‑151 (MQ=255)                                                                                                                                 ttCGATTTCCGGCCCGGTCGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCGGATGCAGACCGCCGTTTTCAATGCGACGACCTGGCAGAGAGACATCAATCGTTTCCGCCGCCAGACgcgc               >  1:1254992/1‑150 (MQ=255)                                                                                                                                  tCGATTTCCGGCCCGGTCGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCGGATGCAGACCGCCGTTTTCAATGCGACGACCTGGCAGAGAGACATCAATCGTTTCCGCCGCCAGACGCGCa              <  2:548651/150‑1 (MQ=255)                                                                                                                                  tCGATTTCCGGCCCGGTCGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCGGATGCAGACCGCCGTTTTCAATGCGACGACCTGGCAGAGAGACATCAATCGTTTCCGCCGCCAGACGCGCAt             >  1:196196/1‑151 (MQ=255)                                                                                                                                           ggCCCGGTCGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCGGATGCAGACCGCCGTTTTCAATGCGACGACCTGGCAGAGAGACATCAATCGTTTCCGCCGCCAGACGCGCATTCAGTGCAg    <  1:503135/151‑1 (MQ=255)                                                                                                                                             cccGGTCGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCGGATGCAGACCGCCGTTTTCAATGCGACGACCTGGCAGAGAGACATCAATCGTTTCCGCCGCCAGACGCGCATTCAGTGCAgcg  >  1:1103645/1‑151 (MQ=255)                                                                                                                                                   |                                                                                                                                                 GTGCGGATCTGTACGCCAGAGGTCTGGGTACGCAGCAGGCGGGTAGTGTCAAACCAGAAAGTGTCGTGGTCAGCGCGCGCCGGGTGGTGACCAGGAATGTTCAGAGCATCGAAGTTATGATAATCGTCTTCGATTTCCGGCCCGGTTGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCGGATGCAGACCGCCGTTTTCAATGCGACGACCTGGCAGAGAGACATCAATCGTTTCCGCCGCCAGACGCGCATTCAGTGCAGCG  >  NC_000913/1798413‑1798703

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