Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,654,069	G→C	T358S (ACT→AGT)	rspA ←	bifunctional D‑altronate/D‑mannonate dehydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,654,069	0	G	C	100.0%	263.3 / NA	65	T358S (ACT→AGT)	rspA	bifunctional D‑altronate/D‑mannonate dehydratase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (30/35);  total (30/35)

CTCCTTACCAGTTCCACAGCGTGCCATCTTCCAGACGTGCGACTGGTAGATAAGCAGGTTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAAGTCCAGTTGTGCGGGAAGACTTCGAGCATTTGTTCGGAATAACCCATGTATTCCTGGACACCGAAATTGGGGACCCACAGATCAAAGTGCAGCGCCGCAGCCATGCAGACTGGTGACAAATCGGAAGGACCGTGTGAGCCAGTACGTA  >  NC_000913/1653923‑1654216                                                                                                                                                   |                                                                                                                                                    cTCCTTACCAGTTCCACAGCGTGCCATCTTCCAGACGTGCGACTGGTAGATAAGCAGGTTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAACTcc                                                                                                                                                  <  1:570038/150‑1 (MQ=255)  tCCTTACCAGTTCCACAGCGTGCCATCTTCCAGACGTGCGACTGGTAGATAAGCAGGTTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAACTcc                                                                                                                                                  >  1:482930/1‑149 (MQ=255)    cTTACCAGTTCCACAGCGTGCCATCTTCCAGACGTGCGACTGGTAGATAAGCAGGTTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAACTCCAGtt                                                                                                                                              >  1:780262/1‑151 (MQ=255)     ttACCAGTTCCACAGCGTGCCATCTTCCAGACGTGCGACTGGTAGATAAGCAGGTTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAACTCCAGTtg                                                                                                                                             >  2:2014/1‑151 (MQ=255)       aCCAGTTCCACAGCGTGCCATCTTCCAGACGTGCGACTGGTAGATAAGCAGGTTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAACTCCAGTtgtg                                                                                                                                           >  1:727277/1‑151 (MQ=255)             tCCACAGCGTGCCATCTTCCAGACGTGCGACTGGTAGATAAGCAGGTTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAACTCCAGTTGTGCGGGa                                                                                                                                      <  2:282941/150‑1 (MQ=255)              ccACAGCGTGCCATCTTCCAGACGTGCGACTGGTAGATAAGCAGGTTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAACTCCAGTTGTGCGGGAAg                                                                                                                                    <  1:265037/151‑1 (MQ=255)                  aGCGTGCCATCTTCCAGACGTGCGACTGGTAGATAAGCAGGTTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAACTCCAGTTGTGCGGGAAGACtt                                                                                                                                <  1:21330/151‑1 (MQ=255)                    cGTGCCATCTTCCAGACGTGCGACTGGTAGATAAGCAGGTTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAACTCCAGTTGTGCGGGAAGACTTCg                                                                                                                              >  2:284877/1‑151 (MQ=255)                          aTCTTCCAGACGTGCGACTGGTAGATAAGCAGGTTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAACTCCAGTTGTGCGGGAAGACTTCGAGCAtt                                                                                                                        >  1:189701/1‑151 (MQ=255)                            cTTCCAGACGTGCGACTGGTAGATAAGCAGGTTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAACTCCAGTTGTGCGGGAAGACTTCGAGCATTTg                                                                                                                      >  1:723758/1‑151 (MQ=255)                              tCCAGACGTGCGACTGGTAGATAAGCAGGTTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAACTCCAGTTGTGCGGGAAGACTTCGAGCAttt                                                                                                                       <  2:1029775/148‑1 (MQ=255)                              tCCAGACGTGCGACTGGTAGATAAGCAGGTTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAACTCCAGTTGTGCGGGAAGACTTCGAGCAttt                                                                                                                       <  1:730093/148‑1 (MQ=255)                              tCCAGACGTGCGACTGGTAGATAAGCAGGTTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAACTCCAGTTGTGCGGGAAGACTTCGAGCAtt                                                                                                                        <  2:499563/147‑1 (MQ=255)                               ccAGACGTGCGACTGGTAGATAAGCAGGTTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAACTCCAGTTGTGCGGGAAGACTTCGAGCATTTGtt                                                                                                                    <  1:770704/150‑1 (MQ=255)                               ccAGACGTGCGACTGGTAGATAAGCAGGTTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAACTCCAGTTGTGCGGGAAGACTTCGAGCATTTGTTc                                                                                                                   >  2:909284/1‑151 (MQ=255)                                   aCGTGCGACTGGTAGATAAGCAGGTTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAACTCCAGTTGTGCGGGAAGACTTCGAGCATTTGTTCGGaa                                                                                                               <  1:411824/151‑1 (MQ=255)                                    cGTGCGACTGGTAGATAAGCAGGTTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAACTCCAGTTGTGCGGGAAGACTTCGAGCATTTGTTCGGAAt                                                                                                              <  2:490929/151‑1 (MQ=255)                                    cGTGCGACTGGTAGATAAGCAGGTTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAACTCCAGTTGTGCGGGAAGACTTCGAGCATTTGTTCGGAAt                                                                                                              <  2:935546/151‑1 (MQ=255)                                         gACTGGTAGATAAGCAGGTTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAACTCCAGTTGTGCGGGAAGACTTCGAGCATTTGTTCGGAATAAccc                                                                                                         >  1:766210/1‑151 (MQ=255)                                           cTGGTAGATAAGCAGGTTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAACTCCAGTTGTGCGGGAAGACTTCGAGCATTTGTTCGGAATAACCCAt                                                                                                       >  2:810616/1‑151 (MQ=255)                                           cTGGTAGATAAGCAGGTTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAACTCCAGTTGTGCGGGAAGACTTCGAGCATTTGTTCGGAATAACCCAt                                                                                                       >  1:973867/1‑151 (MQ=255)                                            tGGTAGATAAGCAGGTTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAACTCCAGTTGTGCGGGAAGACTTCGAGCATTTGTTCGGAATAACCCATg                                                                                                      >  1:1062489/1‑151 (MQ=255)                                              gTAGATAAGCAGGTTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAACTCCAGTTGTGCGGGAAGACTTCGAGCATTTGTTCGGAATAACCCATGTa                                                                                                    <  1:929054/151‑1 (MQ=255)                                               tAGATAAGCAGGTTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAACTCCAGTTGTGCGGGAAGACTTCGAGCATTTGTTCGGAATAACCCATGTAt                                                                                                   >  1:775368/1‑151 (MQ=255)                                                  aTAAGCAGGTTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAACTCCAGTTGTGCGGGAAGACTTCGAGCATTTGTTCGGAATAACCCATGTATTcc                                                                                                <  2:431530/151‑1 (MQ=255)                                                       cAGGTTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAACTCCAGTTGTGCGGGAAGACTTCGAGCATTTGTTCGGAATAACCCATGTATTCCTGGac                                                                                           <  2:647566/151‑1 (MQ=255)                                                                aGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAACTCCAGTTGTGCGGGAAGACTTCGAGCATTTGTTCGGAATAACCCATGTATTCCTGGACACCGAAAtt                                                                                  >  2:317521/1‑151 (MQ=255)                                                                 gggATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAACTCCAGTTGTGCGGGAAGACTTCGAGCATTTGTTCGGAATAACCCATGTATTCCTGGACACCGAAATTg                                                                                 <  1:175378/151‑1 (MQ=255)                                                                  ggATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAACTCCAGTTGTGCGGGAAGACTTCGAGCATTTGTTCGGAATAACCCATGTATTCCTGGACACCGAAATTgg                                                                                >  1:640347/1‑151 (MQ=255)                                                                   gATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAACTCCAGTTGTGCGGGAAGACTTCGAGCATTTGTTCGGAATAACCCATGTATTCCTGGACACCGAAATTggg                                                                               <  2:52312/151‑1 (MQ=255)                                                                     tatTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAACTCCAGTTGTGCGGGAAGACTTCGAGCATTTGTTCGGAATAACCCATGTATTCCTGGACACCGAAATTgggg                                                                              >  2:4289/1‑150 (MQ=255)                                                                       tttCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAACTCCAGTTGTGCGGGAAGACTTCGAGCATTTGTTCGGAATAACCCATGTATTCCTGGACACCGAAATTGGGGAcc                                                                           >  1:815822/1‑151 (MQ=255)                                                                          cgccgcCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAACTCCAGTTGTGCGGGAAGACTTCGAGCATTTGTTCGGAATAACCCATGTATTCCTGGACACCGAAATTGGGGACCca                                                                         <  1:930253/150‑1 (MQ=255)                                                                           gccgccAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAACTCCAGTTGTGCGGGAAGACTTCGAGCATTTGTTCGGAATAACCCATGTATTCCTGGACACCGAAATTGGGGACCcaca                                                                       <  2:1052856/151‑1 (MQ=255)                                                                              gccAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAACTCCAGTTGTGCGGGAAGACGTCGAGCATTTGTTCGGAATAACCCATGTATTCCTGGACACCGAAATTGGGGACCCACAGAt                                                                    <  2:46147/151‑1 (MQ=255)                                                                                     tttCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAACTCCAGTTGTGCGGGAAGACTTCGAGCATTTGTTCGGAATAACCCATGTATTCCTGGACACCGAAATTGGGGACCCACAGATCAAAg                                                               <  1:153874/149‑1 (MQ=255)                                                                                     tttCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAACTCCAGTTGTGCGGGAAGACTTCGAGCATTTGTTCGGAATAACCCATGTATTCCTGGACACCGAAATTGGGGACCCACAGATCAAAGt                                                              <  2:1020773/150‑1 (MQ=255)                                                                                     tttCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAACTCCAGTTGTGCGGGAAGACTTCGAGCATTTGTTCGGAATAACCCATGTATTCCTGGACACCGAAATTGGGGACCCACAGATCAAAGt                                                              <  2:294976/150‑1 (MQ=255)                                                                                         aTCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAACTCCAGTTGTGCGGGAAGACTTCGAGCATTTGTTCGGAATAACCCATGTATTCCTGGACACCGAAATTGGGGACCCACAGATCAAAGTGCAgc                                                         <  2:1052520/151‑1 (MQ=255)                                                                                             aaTTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAACTCCAGTTGTGCGGGAAGACTTCGAGCATTTGTTCGGAATAACCCATGTATTCCTGGACACCGAAATTGGGGACCCACAGATCAAAGTGCAGcgccg                                                     <  2:640487/151‑1 (MQ=255)                                                                                             aaTTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAACTCCAGTTGTGCGGGAAGACTTCGAGCATTTGTTCGGAATAACCCATGTATTCCTGGACACCGAAATTGGGGACCCACAGATCAAAGTGCAGcgccg                                                     >  1:5910/1‑151 (MQ=255)                                                                                             aaTTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAACTCCAGTTGTGCGGGAAGACTTCGAGCATTTGTTCGGAATAACCCATGTATTCCTGGACACCGAAATTGGGGACCCACAGATCAAAGTGCAGcgccg                                                     >  1:559146/1‑151 (MQ=255)                                                                                                 cGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAACTCCAGTTGTGCGGGAAGACTTCGAGCATTTGTTCGGAATAACCCATGTATTCCTGGACATCGAAATTGGGGACCCACAGATCAAAGTGCAGCGCCGCAGc                                                 >  2:712360/1‑151 (MQ=255)                                                                                                     aCCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAACTCCAGTTGTGCGGGAAGACTTCGAGCATTTGTTCGGAATAACCCATGTATTCCTGGACACCGAAATTGGGGACCCACAGATCAAAGTGCAGCGCCGCAGCCATg                                             >  1:453529/1‑151 (MQ=255)                                                                                                             ccGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAACTCCAGTTGTGCGGGAAGACTTCGAGCATTTGTTCGGAATAACCCATGTATTCCTGGACACCGAAATTGGGGACCCACAGATCAAAGTGCAGCGCCGCAGCCATGCAGACTgg                                     <  1:843576/151‑1 (MQ=255)                                                                                                              cGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAACTCCAGTTGTGCGGGAAGACTTCGAGCATTTGTTCGGAATAACCCATGTATTCCTGGACACCGAAATTGGGGACCCACAGATCAAAGTGCAGCGCCGCAGCCATGCAGACTGGt                                    <  2:83758/151‑1 (MQ=255)                                                                                                               ggTTTGTCTCCCGGATGCATATAGCCGTTATCGAAACTCCAGTTGTGCGGGAAGACTTCGAGCATTTGTTCGGAATAACCCATGTATTCCTGGACACCGAAATTGGGGACCCACAGATCAAAGTGCAGCGCCGCAGCCATGCAGACTGGTg                                   <  2:254068/151‑1 (MQ=255)                                                                                                                gTTTGTCTCCCGGATGCATATAGCCGTTATCGAAACTCCAGTTGTGCGGGAAGACTTCGAGCATTTGTTCGGAATAACCCATGTATTCCTGGACACCGAAATTGGGGACCCACAGATCAAAGTGCAGCGCCGCAGCCATGCAGACTGGTGa                                  <  1:226433/151‑1 (MQ=255)                                                                                                                gTTTGTCTCCCGGATGCATATAGCCGTTATCGAAACTCCAGTTGTGCGGGAAGACTTCGAGCATTTGTTCGGAATAACCCATGTATTCCTGGACACCGAAATTGGGGACCCACAGATCAAAGTGCAGCGCCGCAGCCATGCAGACTGGTGa                                  <  1:187675/151‑1 (MQ=255)                                                                                                                    gTCTCCCGGATGCATATAGCCGTTATCGAAACTCCAGTTGTGCGGGAAGACTTCGAGCATTTGTTCGGAATAACCCATGTATTCCTGGACACCGAAATTGGGGACCCACAGATCAAAGTGCAGCGCCGCAGCCATGCAGACTGGTGACaa                               <  1:353810/150‑1 (MQ=255)                                                                                                                     tctcCCGGATGCATATAGCCGTTATCGAAACTCCAGTTGTGCGGGAAGACTTCGAGCATTTGTTCGGAATAACCCATGTATTCCTGGACACCGAAATTGGGGACCCACAGATCAAAGTGCAGCGCCGCAGCCATGCAGACTGGTGACAAAt                             >  1:687605/1‑151 (MQ=255)                                                                                                                      ctcCCGGATGCATATAGCCGTTATCGAAACTCCAGTTGTGCGGGAAGACTTCGAGCATTTGTTCGGAATAACCCATGTATTCCTGGACACCGAAATTGGGGACCCACAGATCAAAGTGCAGCGCCGCAGCCATGCAGACTGGTGACAAATc                            >  1:1004819/1‑151 (MQ=255)                                                                                                                        cccGGATGCATATAGCCGTTATCGAAACTCCAGTTGTGCGGGAAGACTTCGAGCATTTGTTCGGAATAACCCATGTATTCCTGGACACCGAAATTGGGGACCCACAGATCAAAGTGCAGCGCCGCAGCCATGCAGACTGGTGACAAATCgg                          >  1:195779/1‑151 (MQ=255)                                                                                                                         ccGGATGCATATAGCCGTTATCGAAACTCCAGTTGTGCGGGAAGACTTCGAGCATTTGTTCGGAATAACCCATGTATTCCTGGACACCGAAATTGGGGACCCACAGATCAAAGTGCAGCGCCGCAGCCATGCAGACTGGTGACAAATCgg                          >  1:497402/1‑150 (MQ=255)                                                                                                                         ccGGATGCATATAGCCGTTATCGAAACTCCAGTTGTGCGGGAAGACTTCGAGCATTTGTTCGGAATAACCCATGTATTCCTGGACACCGAAATTGGGGACCCACAGATCAAAGTGCAGCGCCGCAGCCATGCAGACTGGTGACAAATCGGa                         <  2:779967/151‑1 (MQ=255)                                                                                                                           ggATGCATATAGCCGTTATCGAAACTCCAGTTGTGCGGGAAGACTTCGAGCATTTGTTCGGAATAACCCATGTATTCCTGGACACCGAAATTGGGGACCCACAGATCAAAGTGCAGCGCCGCAGCCATGCAGACTGGTGACAAATCGGAAg                       <  1:317720/151‑1 (MQ=255)                                                                                                                            gATGCATATAGCCGTTATCGAAACTCCAGTTGTGCGGGAAGACTTCGAGCATTTGTTCGGAATAACCCATGTATTCCTGGACACCGAAATTGGGGACCCACAGATCAAAGTGCAGCGCCGCAGCCATGCAGACTGGTGACAAATCGGAAgg                      >  1:826861/1‑151 (MQ=255)                                                                                                                              tGCATATAGCCGTTATCGAAACTCCAGTTGTGCGGGAAGACTTCGAGCATTTGTTCGGAATAACCCATGTATTCCTGGACACCGAAATTGGGGACCCACAGATCAAAGTGCAGCGCCGCAGCCATGCAGACTGGTGACAAATCGGAAGGAc                    >  2:522422/1‑151 (MQ=255)                                                                                                                               gCATATAGCCGTTATCGAAACTCCAGTTGTGCGGGAAGACTTCGAGCATTTGTTCGGAATAACCCATGTATTCCTGGACACCGAAATTGGGGACCCACAGATCAAAGTGCAGCGCCGCAGCCATGCAGACTGGTGACAAATCGGAAGGAcc                   >  2:235510/1‑151 (MQ=255)                                                                                                                                   ataGCCGTTATCGAAACTCCAGTTGTGCGGGAAGACTTCGAGCATTTGTTCGGAATAACCCATGTATTCCTGGACACCGAAATTGGGGACCCACAGATCAAAGTGCAGCGCCGCAGCCATGCAGACTGGTGACAAATCGGAAGGACCgtg                <  2:344409/150‑1 (MQ=255)                                                                                                                                     aGCCGTTATCGAAACTCCAGTTGTGCGGGAAGACTTCGAGCATTTGTTCGGAATAACCCATGTATTCCTGGACACCGAAATTGGGGACCCACAGATCAAAGTGCAGCGCCGCAGCCATGCAGACTGGTGACAAATCGGAAGGACCGtgtg              >  1:26432/1‑150 (MQ=255)                                                                                                                                         gTTATCGAAACTCCAGTTGTGCGGGAAGACTTCGAGCATTTGTTCGGAATAACCCATGTATTCCTGGACACCGAAATTGGGGACCCACAGATCAAAGTGCAGCGCCGCAGCCATGCAGACTGGTGACAAATCGGAAGGACCGTGTGAGcc          >  1:1023296/1‑150 (MQ=255)                                                                                                                                            aTCGAAACTCCAGTTGTGCGGGAAGACTTCGAGCATTTGTTCGGAATAACCCATGTATTCCTGGACACCGAAATTGGGGACCCACAGATCAAAGTGCAGCGCCGCAGCCATGCAGACTGGTGACAAATCGGAAGGACCGTGTGAGCCAgta      <  2:629661/151‑1 (MQ=255)                                                                                                                                                aaaCTCCAGTTGTGCGGGAAGACTTCGAGCATTTGTTCGGAATAACCCATGTATTCCTGGACACCGAAATTGGGGACCCACAGATCAAAGTGCAGCGCCGCAGCCATGCAGACTGGTGACAAATCGGAAGGACCGTGTGAGCCAgtacgta  <  1:520391/151‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                    CTCCTTACCAGTTCCACAGCGTGCCATCTTCCAGACGTGCGACTGGTAGATAAGCAGGTTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAAGTCCAGTTGTGCGGGAAGACTTCGAGCATTTGTTCGGAATAACCCATGTATTCCTGGACACCGAAATTGGGGACCCACAGATCAAAGTGCAGCGCCGCAGCCATGCAGACTGGTGACAAATCGGAAGGACCGTGTGAGCCAGTACGTA  >  NC_000913/1653923‑1654216

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