Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,337,147	T→C	S439S (TCT→TCC)	acnA →	aconitate hydratase 1; aconitase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,337,147	0	T	C	100.0%	271.6 / NA	72	S439S (TCT→TCC)	acnA	aconitate hydratase 1; aconitase A
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (40/32);  total (40/32)

AAAAGCATTTGCCGCCAGTAACGAACTGGAAGTGAATGCCACGCATAAAGATCGCCAGCCGGTCGATTATGTTATGAACGGACATCAGTATCAGTTACCTGATGGCGCTGTGGTCATTGCTGCGATAACCTCGTGCACCAACACCTCTAACCCAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGT  >  NC_000913/1337000‑1337292                                                                                                                                                    |                                                                                                                                                  aaaaGCATTTGCCGCCAGTAACGAACTGGAAGTGAATGCCACGCATAAAGATCGCCAGCCGGTCGATTATGTTATGAACGGACATCAGTATCAGTTACCTGATGGCGCTGTGGTCATTGCTGCGATAACCTCGTGCACCAACACCTCCAAc                                                                                                                                                <  1:435050/151‑1 (MQ=255)  aaaGCATTTGCCGCCAGTAACGAACTGGAAGTGAATGCCACGCATAAAGATCGCCAGCCGGTCGATTATGTTATGAACGGACATCAGTATCAGTTACCTGATGGCGCTGTGGTCATTGCTGCGATAACCTCGTGCACCAACACCTCCAAcc                                                                                                                                               <  2:462637/151‑1 (MQ=255)     gCATTTGCCGCCAGTAACGAACTGGAAGTGAATGCCACGCATAAAGATCGCCAGCCGGTCGATTATGTTATGAACGGACATCAGTATCAGTTACCTGATGGCGCTGTGGTCATTGCTGCGATAACCTCGTGCACCAACACCTCCAACCCaa                                                                                                                                            <  2:739137/151‑1 (MQ=255)     gCATTTGCCGCCAGTAACGAACTGGAAGTGAATGCCACGCATAAAGATCGCCAGCCGGTCGATTATGTTATGAACGGACATCAGTATCAGTTACCTGATGGCGCTGTGGTCATTGCTGCGATAACCTCGTGCACCAACACCTCCAACCCa                                                                                                                                             >  2:730765/1‑150 (MQ=255)     gCATTTGCCGCCAGTAACGAACTGGAAGTGAATGCCACGCATAAAGATCGCCAGCCGGTCGATTATGTTATGAACGGACATCAGTATCAGTTACCTGATGGCGCTGTGGTCATTGCTGCGATAACCTCGTGCACCAACACCTCCAACCCa                                                                                                                                             >  1:325676/1‑150 (MQ=255)       aTTTGCCGCCAGTAACGAACTGGAAGTGAATGCCACGCATAAAGATCGCCAGCCGGTCGATTATGTTATGAACGGACATCAGTATCAGTTACCTGATGGCGCTGTGGTCATTGCTGCGATAACCTCGTGCACCAACACCTACAACCCAAgt                                                                                                                                          <  2:467663/151‑1 (MQ=255)          tGCCGCCAGTAACGAACTGGAAGTGAATGCCACGCATAAAGATCGCCAGCCGGTCGATTATGTTATGAACGGACATCAGTATCAGTTACCTGATGGCGCTGTGGTCATTGCTGCGATAACCTCGTGCACCAACACCTCCAACCCAAgtgt                                                                                                                                        <  1:164211/150‑1 (MQ=255)           gccgccAGTAACGAACTGGAAGTGAATGCCACGCATAAAGATCGCCAGCCGGTCGATTATGTTATGAACGGACATCAGTATCAGTTCCCTGCTGGCGCTGTGATCATTGCTGCGATAACCTCGCGCACCAACACCTCCAGCCCAAGtgtg                                                                                                                                       <  1:107533/150‑1 (MQ=255)            ccgccAGTAACGAACTGGAAGTGAATGCCACGCATAAAGATCGCCAGCCGGTCGATTATGTTATGAACGGACATCAGTATCAGTTACCTGATGGCGCTGTGGTCATTGCTGCGATAACCTCGTGCACCAACACCTCCAACCCAAGTGTGCt                                                                                                                                     <  2:850450/151‑1 (MQ=255)              gccAGTAACGAACTGGAAGTGAATGCCACGCATAAAGATCGCCAGCCGGTCGATTATGTTATGAACGGACATCAGTATCAGTTACCTGATGGCGCTGTGGTCATTGCTGCGATAACCTCGTGCACCAACACCTCCAACCCAAGTGTGCTg                                                                                                                                    >  1:332084/1‑150 (MQ=255)                  gTAACGAACTGGAAGTGAATGCCACGCATAAAGATCGCCAGCCGGTCGATTATGTTATGAACGGACATCAGTATCAGTTACCTGATGGCGCTGTGGTCATTGCTGCGATAACCTCGTGCACCAACACCTCCAACCCAAGTGTGCTGATgg                                                                                                                                <  2:1193142/150‑1 (MQ=255)                     aCGAACTGGAAGTGAATGCCACGCATAAAGATCGCCAGCCGGTCGATTATGTTATGAACGGACATCAGTATCAGTTACCTGATGGCGCTGTGGTCATTGCTGCGATAACCTCGTGCACCAACACCTCCAACCCAAGTGTGCTGATGGCCGc                                                                                                                            >  2:944225/1‑151 (MQ=255)                      cGAACTGGAAGTGAATGCCACGCATAAAGATCGCCAGCCGGTCGATTATGTTATGAACGGACATCAGTATCAGTTACCTGATGGCGCTGTGGTCATTGCTGCGATAACCTCGTGCACCAACACCTCCAACCCAAGTGTGCTGATGGCCGCa                                                                                                                           <  2:515446/151‑1 (MQ=255)                        aaCTGGAAGTGAATGCCACGCATAAAGATCGCCAGCCGGTCGATTATGTTATGAACGGACATCAGTATCAGTTACCTGATGGCGCTGTGGTCATTGCTGCGATAACCTCGTGCACCAACACCTCCAACCCAAGTGTGCTGATGGCCGCAgg                                                                                                                         >  2:1010011/1‑151 (MQ=255)                             gAAGTGAATGCCACGCATAAAGATCGCCAGCCGGTCGATTATGTTATGAACGGACATCAGTATCAGTTACCTGATGGCGCTGTGGTCATTGCTGCGATAACCTCGTGCACCAACACCTCCAACCCAAGTGTGCTGATGGCCGCAGGCTTGc                                                                                                                    >  2:766344/1‑151 (MQ=255)                              aaGTGAATGCCACGCATAAAGATCGCCAGCCGGTCGATTATGTTATGAACGGACATCAGTATCAGTTACCTGATGGCGCTGTGGTCATTGCTGCGATAACCTCGTGCACCAACACCTCCAACCCAAGTGTGCTGATGGCCGCAGGCTTGc                                                                                                                    <  2:980086/150‑1 (MQ=255)                              aaGTGAATGCCACGCATAAAGATCGCCAGCCGGTCGATTATGTTATGAACGGACATCAGTATCAGTTACCTGATGGCGCTGTGGTCATTGCTGCGATAACCTCGTGCACCAACACCTCCAACCCAAGTGTGCTGATGGCCGCAGGCTTGCt                                                                                                                   >  2:687630/1‑151 (MQ=255)                              aaGTGAATGCCACGCATAAAGATCGCCAGCCGGTCGATTATGTTATGAACGGACATCAGTATCAGTTACCTGATGGCGCTGTGGTCATTGCTGCGATAACCTCGTGCACCAACACCTCCAACCCAAGTGTGCTGATGGCCGCAGGCTTGCt                                                                                                                   >  2:1184276/1‑151 (MQ=255)                                  gAATGCCACGCATAAAGATCGCCAGCCGGTCGATTATGTTATGAACGGACATCAGTATCAGTTACCTGATGGCGCTGTGGTCATTGCTGCGATAACCTCGTGCACCAACACCTCCAACCCAAGTGTGCTGATGGCCGCAGGCTTGCTGGCg                                                                                                               <  1:587874/151‑1 (MQ=255)                                   aaTGCCACGCATAAAGATCGCCAGCCGGTCGATTATGTTATGAACGGACATCAGTATCAGTTACCTGATGGCGCTGTGGTCATTGCTGCGATAACCTCGTGCACCAACACCTCCAACCCAAGTGTGCTGATGGCCGCAGGCTTGCTGGCg                                                                                                               >  2:824158/1‑150 (MQ=255)                                   aaTGCCACGCATAAAGATCGCCAGCCGGTCGATTATGTTATGAACGGACATCAGTATCAGTTACCTGATGGCGCTGTGGTCATTGCTGCGATAACCTCGTGCACCAACACCTCCAACCCAAGTGTGCTGATGGCCGCAGGCTTGCTGGCg                                                                                                               >  1:173856/1‑150 (MQ=255)                                   aaTGCCACGCATAAAGATCGCCAGCCGGTCGATTATGTTATGAACGGACATCAGTATCAGTTACCTGATGGCGCTGTGGTCATTGCTGCGATAACCTCGTGCACCAACACCTCCAACCCAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGa                                                                                                              <  2:1129462/151‑1 (MQ=255)                                   aaTGCCACGCATAAAGATCGCCAGCCGGTCGATTATGTTATGAACGGACATCAGTATCAGTTACCTGATGGCGCTGTGGTCATTGCTGCGATAACCTCGTGCACCAACACCTCCAACCCAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGa                                                                                                              <  2:111543/151‑1 (MQ=255)                                        cACGCATAAAGATCGCCAGCCGGTCGATTATGTTATGAACGGACATCAGTATCAGTTACCTGATGGCGCTGTGGTCATTGCTGCGATAACCTCGTGCACCAACACCTCCAACCCAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGaaaaaa                                                                                                         <  2:807508/151‑1 (MQ=255)                                         aCGCATAAAGATCGCCAGCCGGTCGATTATGTTATGAACGGACATCAGTATCAGTTACCTGATGGCGCTGTGGTCATTGCTGCGATAACCTCGTGCACCAACACCTCCAACCCAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAg                                                                                                        <  2:523149/151‑1 (MQ=255)                                                  gATCGCCAGCCGGTCGATTATGTTATGAACGGACATCAGTATCAGTTACCTGATGGCGCTGTGGTCATTGCTGCGATAACCTCGTGCACCAACACCTCCAACCCAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAActc                                                                                               >  1:134587/1‑151 (MQ=255)                                                  gATCGCCAGCCGGTCGATTATGTTATGAACGGACATCAGTATCAGTTACCTGATGGCGCTGTGGTCATTGCTGCGATAACCTCGTGCACCAACACCTCCAACCCAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAActc                                                                                               >  2:154312/1‑151 (MQ=255)                                                  gATCGCCAGCCGGTCGATTATGTTATGAACGGACATCAGTATCAGTTACCTGATGGCGCTGTAGTCATTGCTGCGATAACCTCGTGCACCAACACCTCCAACCCAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAActc                                                                                               >  1:1086748/1‑151 (MQ=255)                                                   aTCGCCAGCCGGTCGATTATGTTATGAACGGACATCAGTATCAGTTACCTGATGGCGCTGTGGTCATTGCTGCGATAACCTCGTGCACCAACACCTCCAACCCAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAActct                                                                                              >  2:472221/1‑151 (MQ=255)                                                   aTCGCCAGCCGGTCGATTATGTTATGAACGGACATCAGTATCAGTTACCTGATGGCGCTGTGGTCATTGCTGCGATAACCTCGTGCACCAACACCTCCAACCCAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAActct                                                                                              <  2:272002/151‑1 (MQ=255)                                                     cGCCAGCCGGTCGATTATGTTATGAACGGACATCAGTATCAGTTACCTGATGGCGCTGTGGTCATTGCTGCGATAACCTCGTGCACCAACACCTCCAACCCAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTgg                                                                                            >  2:39260/1‑151 (MQ=255)                                                     cGCCAGCCGGTCGATTATGTTATGAACGGACATCAGTATCAGTTACCTGATGGCGCTGTGGTCATTGCTGCGATAACCTCGTGCACCAACACCTCCAACCCAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTgg                                                                                            <  1:944664/151‑1 (MQ=255)                                                      gCCAGCCGGTCGATTATGTTATGAACGGACATCAGTATCAGTTACCTGATGGCGCTGTGGTCATTGCTGCGATAACCTCGTGCACCAACACCTCCAACCCAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTggg                                                                                           >  2:1115306/1‑151 (MQ=255)                                                      gCCAGCCGGTCGATTATGTTATGAACGGACATCAGTATCAGTTACCTGATGGCGCTGTGGTCATTGCTGCGATAACCTCGTGCACCAACACCTCCAACCCAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTggg                                                                                           <  1:489277/151‑1 (MQ=255)                                                         aGCCGGTCGATTATGTTATGAACGGACATCAGTATCAGTTACCTGATGGCGCTGTGGTCATTGCTGCGATAACCTCGTGCACCAACACCTCCAACCCAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCt                                                                                        <  1:1114002/151‑1 (MQ=255)                                                         aGCCGGTCGATTATGTTATGAACGGACATCAGTATCAGTTACCTGATGGCGCTGTGGTCATTGCTGCGATAACCTCGTGCACCAACACCTCCAACCCAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCt                                                                                        <  2:132828/151‑1 (MQ=255)                                                             ggTCGATTATGTTATGAACGGACATCAGTATCAGTTACCTGATGGCGCTGTGGTCATTGCTGCGATAACCTCGTGCACCAACACCTCCAACCCAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCa                                                                                      >  2:835605/1‑149 (MQ=255)                                                                  aTTATGTTATGAACGGACATCAGTATCAGTTACCTGATGGCGCTGTGGTCATTGCTGCGATAACCTCGTGCACCAACACCTCCAACCCAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGc                                                                                >  1:87957/1‑150 (MQ=255)                                                                   ttatgttatgAACGGACATCAGTATCAGTTACCTGATGGCGCTGTGGTCATTGCTGCGATAACCTCGTGCACCAACACCTCCAACCCAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCa                                                                               >  2:739568/1‑150 (MQ=255)                                                                    tatgttatgAACGGACATCAGTATCAGTTACCTGATGGCGCTGTGGTCATTGCTGCGATAACCTCGTGCACCAACACCTCCAACCCAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAA‑‑CTGGGCCTCAAGCGGCAAcc                                                                            >  1:320876/1‑150 (MQ=255)                                                                       gttatgAACGGACATCAGTATCAGTTACCTGATGGCGCTGTGGTCATTGCTGCGATAACCTCGTGCACCAACACCTCCAACCCAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCAt                                                                          <  1:906919/151‑1 (MQ=255)                                                                            gAACGGACATCAGTATCAGTTACCTGATGGCGCTGTGGTCATTGCTGCGATAACCTCGTGCACCAACACCTCCAACCCAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGt                                                                      <  2:849198/150‑1 (MQ=255)                                                                              aCGGACATCAGTATCAGTTACCTGATGGCGCTGTGGTCATTGCTGCGATAACCTCGTGCACCAACACCTCCAACCCAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCa                                                                    >  2:645279/1‑150 (MQ=255)                                                                              aCGGACATCAGTATCAGTTACCTGATGGCGCTGTGGTCATTGCTGCGATAACCTCGTGCACCAACACCTCCAACCCAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCa                                                                    >  2:634695/1‑150 (MQ=255)                                                                              aCGGACATCAGTATCAGTTACCTGATGGCGCTGTGGTCATTGCTGCGATAACCTCGTGCACCAACACCTCCAACCCAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCa                                                                    >  1:18596/1‑150 (MQ=255)                                                                               cGGACATCAGTATCAGTTACCTGATGGCGCTGTGGTCATTGCTGCGATAACCTCGTGCACCAACACCTCCAACCCAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCaa                                                                   >  2:598998/1‑150 (MQ=255)                                                                                 gACATCAGTATCAGTTACCTGATGGCGCTGTGGTCATTGCTGCGATAACCTCGTGCACCAACACCTCCAACCCAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAg                                                                 >  1:309050/1‑150 (MQ=255)                                                                                  aCATCAGTATCAGTTACCTGATGGCGCTGTGGTCATTGCTGCGATAACCTCGTGCACCAACACCTCCAACCCAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCg                                                               >  1:923922/1‑151 (MQ=255)                                                                                      cagtatcagtTACCTGATGGCGCTGTGGTCATTGCTGCGATAACCTCGTGCACCAACACCTCCAACCCAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGc                                                           >  2:439183/1‑151 (MQ=255)                                                                                       agtatcagtTACCTGATGGCGCTGTGGTCATTGCTGCGATAACCTCGTGCACCAACACCTCCAACCCAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCt                                                          >  2:289401/1‑151 (MQ=255)                                                                                               ttACCTGATGGCGCTGTGGTCATTGCTGCGATAACCTCGTGCACCAACACCTCCAACCCAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCgg                                                  >  2:933413/1‑151 (MQ=255)                                                                                                tACCTGATGGCGCTGTGGTCATTGCTGCGATAACCTCGTGCACCAACACCTCCAACCCAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCggg                                                 >  2:1077092/1‑151 (MQ=255)                                                                                                 aCCTGATGGCGCTGTGGTCATTGCTGCGATAACCTCGTGCACCAACACCTCCAACCCAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGt                                                >  1:750869/1‑151 (MQ=255)                                                                                                  ccTGATGGCGCTGTGGTCATTGCTGCGATAACCTCGTGCACCAACACCTCCAACCCAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGtt                                               >  1:1124265/1‑151 (MQ=255)                                                                                                  ccTGATGGCGCTGTGGTCATTGCTGCGATAACCTCGTGCACCAACACCTCCAACCCAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGtt                                               >  2:938388/1‑151 (MQ=255)                                                                                                     gATGGCGCTGTGGTCATTGCTGCGATAACCTCGTGCACCAACACCTCCAACCCAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTc                                              <  2:420226/149‑1 (MQ=255)                                                                                                         gcgcTGTGGTCATTGCTGCGATAACCTCGTGCACCAACACCTCCAACCCAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGt                                        >  2:949826/1‑151 (MQ=255)                                                                                                          cgcTGTGGTCATTGCTGCGATAACCTCGTGCACCAACACCTCCAACCCAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTc                                       <  1:739798/151‑1 (MQ=255)                                                                                                             tgtgGTCATTGCTGCGATAACCTCGTGCACCAACACCTCCAACCCAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGtt                                    >  2:163029/1‑151 (MQ=255)                                                                                                              gtgGTCATTGCTGCGATAACCTCGTGCACCAACACCTCCAACCCAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGtt                                    <  1:1179303/150‑1 (MQ=255)                                                                                                              gtgGTCATTGCTGCGATAACCTCGTGCACCAACACCTCCAACCCAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGtt                                    <  2:87981/150‑1 (MQ=255)                                                                                                                 gTCATCGCTGTGATCCCCCAGTGCACTAACACCTCCAACCAAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTCACTCTGGGCCTCAAGCGGCAACCATGGGTCACAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTg                                <  1:852185/151‑1 (MQ=255)                                                                                                                  tCATTGCTGCGATAACCTCGTGCACCAACACCTCCAACCCAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTg                                >  2:1191889/1‑150 (MQ=255)                                                                                                                      tgctgcGATAACCTCGTGCACCAACACCTCCAACCCAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTa                            <  2:655502/150‑1 (MQ=255)                                                                                                                          gcGATAACCTCGTGCACCAACACCTCCAACCCAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCTTTTCTGATTGTCt                         <  1:56466/149‑1 (MQ=255)                                                                                                                           cGATAACCTCGTGCACCAACACCTCCAACCCAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGc                      <  1:730597/151‑1 (MQ=255)                                                                                                                               aaCCTCGTGCACCAACACCTCCAACCCAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCaaaa                  <  2:893799/151‑1 (MQ=255)                                                                                                                               aaCCTCGTGCACCAACACCTCCAACCCAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCaaaa                  <  1:1040873/151‑1 (MQ=255)                                                                                                                                  cTCGTGCACCAACACCTCCAACCCAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCg               >  1:247159/1‑151 (MQ=255)                                                                                                                                   tCGTGCACCAACACCTCCAACCCAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCg               >  1:479796/1‑150 (MQ=255)                                                                                                                                          ccAACACCTCCAACCCAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGac       >  1:203317/1‑151 (MQ=255)                                                                                                                                               acCTCCAACCCAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGt  >  1:196529/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                                  AAAAGCATTTGCCGCCAGTAACGAACTGGAAGTGAATGCCACGCATAAAGATCGCCAGCCGGTCGATTATGTTATGAACGGACATCAGTATCAGTTACCTGATGGCGCTGTGGTCATTGCTGCGATAACCTCGTGCACCAACACCTCTAACCCAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGT  >  NC_000913/1337000‑1337292

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