Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,930,868	A→G	F76L (TTC→CTC)	ravA ←	hexameric AAA+ MoxR family ATPase, putative molecular chaperone

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,930,868	0	A	G	100.0%	229.3 / NA	64	F76L (TTC→CTC)	ravA	hexameric AAA+ MoxR family ATPase, putative molecular chaperone
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (35/29);  total (35/29)

GAATTGCCGGGCCCGCTTTCCAGATCTCATCCAGAAAGACAATTTCTGCTTCCGGCAGGTAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGACCAAAAACTTCTTCCGGCGTGGAGAAGCGGGTCATCAGATATTCAAACGCGCGGGCATTCTGAAAGGCGAATTTTAAGCGCCGGGCGATCAAACTTTTGGCAATACCTGGC‑GGGCCAAGGAGGAACACA‑CTTTCACCACTTAATGCCGCTAATAAACACAAGCGGATGGCGTG  >  NC_000913/3930724‑3931015                                                                                                                                                 |                                                                                                                                                      gAATTGCCGGGCCCGCTTTCCAGATCTCATCCAGAAAGACAATTTCTGCTTCCGGCAGGTAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGACCAAAAACTTCTTCCGGCGTGGAGAGGCGGGt                                                                                                                                                 <  1:197443/151‑1 (MQ=255)  aaTTGCCGGGCCCGCTTTCCAGATCTCATCCAGAAAGACAATTTCTGCTTCCGGCAGGTAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGACCAAAAACTTCTTCCGGCGTGGAGAGGCgg                                                                                                                                                   >  1:853812/1‑148 (MQ=255)     tGCCGGGCCCGCTTTCCAGATCTCATCCAGAAAGACAATTTCTGCTTCCGGCAGGTAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGACCAAAAACTTCTTCCGGCGTGGAGAGGCGGGtta                                                                                                                                               >  2:290254/1‑147 (MQ=255)       ccGGGCCCGCTTTCCAGATCTCATCCAGAAAGACAATTTCTGCTTCCGGCAGGTAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGACCAAAAACTTCTTCCGGCGTGGAGAGGCGGGgcatca                                                                                                                                            >  1:439625/1‑150 (MQ=255)         gggCCCGCTTTCCAGATCTCATCCAGAAAGACAATTTCTGCTTCCGGCAGGTAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGACCAAAAACTTCTTCCGGCGTGGAGAGGCGGGTCATCAGat                                                                                                                                         <  1:704828/151‑1 (MQ=255)          ggCCCGCTTTCCAGATCTCATCCAGAAAGACAATTTCTGCTTCCGGCAGGTAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGACCAAAAACTTCTTCCGGCGTGGAGAGGCGGGTCATCAGata                                                                                                                                        <  2:849472/151‑1 (MQ=255)            cccGCTTTCCAGATCTCATCCAGAAAGACAATTTCTGCTTCCGGCAGGTAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGACCAAAAACTTCTTCCGGCGTGGAGAGGCGGGTCATCAGATAtt                                                                                                                                      <  2:492920/151‑1 (MQ=255)             ccGCTTTCCAGATCTCATCCAGAAAGACAATTTCTGCTTCCGGCAGGTAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGACCAAAAACTTCTTCCGGCGTGGAGTGGCGGATTAACAGATAtt                                                                                                                                      >  1:930523/1‑150 (MQ=255)               gCTTTCCAGATCTCATCCAGAAAGACAATTTCTGCTTCCGGCAGGTAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGACCAAAAACTTCTTCCGGCGTGGAGAGGCGGGTGATCAGAAAtt                                                                                                                                      >  1:34467/1‑148 (MQ=255)                cTTTCCAGATCTCATCCAGAAAGACAATTTCTGCTTCCGGCAGGTAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGACCAAAAACTTCTTCCGGCGTGGAGAGGCGGGTCATCAGATAtt                                                                                                                                      >  2:909882/1‑147 (MQ=255)                  ttCCAGATCTCATCCAGAAAGACAATTTCTGCTTCCGGCAGGTAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGACCAAAAACTTCTTCCGGCGTGGAGAGGCGGGGCCATAGAAATTCaaa                                                                                                                                  >  1:248203/1‑149 (MQ=255)                    ccAGATCTCATCCAGAAAGACAATTTCTGCTTCCGGCAGGTAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGACCAAAAACTTCTTCCGGCGTGGAGAGGCGGGTCATCAGATATTCAAAcgcg                                                                                                                              <  2:878212/151‑1 (MQ=255)                     cAGATCTCAGCCAGAAAGACAATGGCTGCTTCCGGCAGGTAAGCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGACCAAAAACTTCTTCCGGCGTGGAGAGGCGGGTCATCAGATATTCAAAcgc                                                                                                                               <  1:694169/149‑1 (MQ=255)                          ctcATCCAGAAAGACAATTTCTGCTTCCGGCAGGTAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGACCAAAAACTTCTTCCGGCGTGGAGAGGCGGGGCATCAGAAAATCAAACGCGCggg                                                                                                                          >  1:933308/1‑149 (MQ=255)                            cATCCAGAAAGACAATTTCTGCTTCCGGCAGGTAACCGCTGGTTAAACGTTCATAGCTCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGACGAAAAGCTTCTGCCGGCGTGGAGAGGCGGTTCATCAGATATTCAAACGCGCGGGCag                                                                                                                       <  1:991065/150‑2 (MQ=255)                                cAGAAAGACAATTTCTGCTTCCGGCAGGTAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGACCAAAAACTTCTTCCGGCGTGGAGAGGCGGGTCATCAGATATTCAAACGCGCGGGCATTCTGa                                                                                                                  <  2:699403/151‑1 (MQ=255)                                cAGAAAGACAATTTCTGCTTCCGGCAGGTAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGACCAAAAACTTCTTCCGGCGTGGAGAGGCGGGTCATCAGAAATTTAAAAGCGCGGGCAATCt                                                                                                                    >  1:1106131/1‑149 (MQ=255)                                 aGAAAGACAATTTCTGCTTCCGGCAGGTAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGACCAAAAACTTCTTCCGGCGTGGAGAGGCGGGGCCTCAGATATTCAAAAACGCGGGGACTCTg                                                                                                                   >  1:788824/1‑149 (MQ=255)                                     aGACAATTTCTGCTTCCGGCAGGTAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGACCAAAAACTTCTTCCGGCGTGGAGAGGCGGGTCATCAGATATTCAAACGCGCGGGCATTCTGAAAGGc                                                                                                             <  1:969904/151‑1 (MQ=255)                                       aCAATTTCTGCTTCCGGCAGGTAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGACCAAAAACTTCTTCCGGCGTGGAGAGGCGGGTCATCAGATATTCAAACGCGCGGGCATTCTGAAAGGCGa                                                                                                           <  1:1210538/151‑1 (MQ=255)                                       aCAATTTCTGCTTCCGGCAGGTAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGACCAAAAACTTCTTCCGGCGTGGAGAGGCGGGTCATCAGATATTCAAACGCGCGGGCATTCTGAAAGGCGa                                                                                                           <  2:779076/151‑1 (MQ=255)                                         aaTTTCTGCTTCCGGCAGGTAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGACCAAAAACTTCTTCCGGCGTGGAGAGGCGGGTCATCAGATATTCAAACGCGCGGGCATTCTGAAAGGCGAAt                                                                                                         <  1:64777/151‑1 (MQ=255)                                         aaTTTCTGCTTCCGGCAGGTAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGACCAAAAACTTCTTCCGGCGTGGAGAGGCGGGTCATCAGATATTCAAACGCGCGGGCATTCTGAAAGGCGAAt                                                                                                         <  2:530439/151‑1 (MQ=255)                                         aaTTTCTGCTTCCGGCAGGTAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGACCAAAAACTTCTTCCGGCGTGGAGAGGCGGGTCATCAGATAATCAAACACGCCGGCATTCTGAAAGGGGa                                                                                                           >  1:266553/1‑149 (MQ=255)                                            ttCTGCTTCCGGCAGGTAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGACCAAAAACTTCTTCCGGCGTGGAGAGGCGGGTCATCAGATATTCAAACGCGCGGGCATTCTGAAAGGCGAAttt                                                                                                       <  1:891605/150‑1 (MQ=255)                                             tCTGCTTCCGGCAGGTAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGACCAAAAACTTCTTCCGGCGTGGAGAGGCGGGTCATCAGATATTCAAACGCGCGGGCATTCTGAAAGGCGAAtttt                                                                                                      >  2:1080873/1‑150 (MQ=255)                                              cTGCTTCCGGCAGGTAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGACCAAAAACTTCTTCCGGCGTGGAGAGGCGGGGCCTCCGATATTTCAAAACGCGGGGGTTTTTAAAGGGGGAtttt                                                                                                      >  1:1001318/1‑149 (MQ=255)                                              cTGCTTCCGGCAGGTAACCGCTGGTTAAACGTTCATAGCG‑CCTTCATCTTTTAGCGCCTGAATAGAAAGGGGACCAAAAACTTCTTCCGGCGTGGAGAGGCGGGTCATCAGATATTCAAACGCGCGGGCATTCTGAAAGGCGAATTTTa                                                                                                     <  2:109212/149‑1 (MQ=255)                                                gCTTCCGGCAGGTAACCGCTGGTTAAATGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGACCAAAAACTTCTTCCGGCGTGGAGAGGCGGGTCATCAGAAATTCAAACACGCGGGCATTCTGAAAGGGGAATTTTTAgc                                                                                                  >  1:1118310/1‑151 (MQ=255)                                                  ttCCGGCAGGTAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGACCAAAAACTTCTTCCGGCGTGGAGAGGCGGGTCATCAGATATTTAAAAGCGCGGGCATTCTGAAAGGGGGAtttttaaag                                                                                                 >  2:24211/1‑145 (MQ=255)                                                    ccGGCAGGTAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGACCAAAAACTTCTTCCGGCGTGGAGAGGCGGGGTCTCAGATATTTAAAAGCGCGGGCATTCTGAAAAGCGAATTTTaatggcc                                                                                               >  2:402871/1‑145 (MQ=255)                                                       gCAGGTAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGACCAAAAACTTCTTCCGGCGTGGAGAGGCGGGTCATCAGATATTCAAACGCGCGGGCATTCTGAAAGGCGAATTTTAAGCGCCGGGc                                                                                           <  2:155308/151‑1 (MQ=255)                                                        cAGGTAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGACCAAAAACTTCTTCCGGCGTGGAGAGGCGGGTCATCAGATATTCAAACGCGCGGGCATTTTGAAAGGGGAATTTTAAGCGCCGGGCg                                                                                          >  1:165777/1‑151 (MQ=255)                                                           gTAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGACCAAAAACTTCTTCCGGCGTGGAGAGGCGGGTCATCAGATATTCAAACGCGCGGGGATTTTGAAAGGCGAATTTTAAGCGCCCGGCGATc                                                                                       >  2:434742/1‑151 (MQ=255)                                                            tAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGACCAAAAACTTCTTCCGGCGTGGAGAGGCGGGTTAATAAATATTTTAACGCGCCGGGCTTTTCTGAGGCGAATTTTTAGCGCCGGGCg                                                                                          >  1:641180/1‑147 (MQ=255)                                                            tAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGACCAAAAACTTCTTCCGGCGTGGAGAGGCGGGTCATCAGATATTCAAACGCGCGGGCATTCTGGAAGGCGAATTTTAAGCGCCGGGCGATc                                                                                       >  1:875608/1‑150 (MQ=255)                                                             aaCCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGACCAAAAACTTCTTCCGGCGTGGAGAGGCGGGGCATCCGATATTCAAACGCGCGGGCATTCTTTAAGGCGAATTTTAAGCGCCGGGCg                                                                                          >  2:899742/1‑146 (MQ=255)                                                              aCCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGACCAAAAACTTCTTCCGGCGTGGAGAGGCGGGGCATCAGATATTCAAACGCGCGGGCATTCTGAAAGGCGAATTTTAAGCGCCGGGCGATCaa                                                                                     >  2:257035/1‑150 (MQ=255)                                                                   tGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGACCAAAAACTTCTTCCGGCGTGGAGAGGCGGGTCATCAGATATTTAAACGCGCGGGCATTCTGAAAGGCGAATTTTAAGCGCCGGGCGATCAAACtttt                                                                               >  2:1157798/1‑151 (MQ=255)                                                                        aaaCGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGACCAAAAACTTCTTCCGGCGTGGAGAGGCGGGTCATCAGATATTCAAACGCGCGGGCATTCTGAAAGGCGAATTTTAAGCGCCGGGCGATCAAACTTTTGGCa                                                                           <  2:1003762/150‑1 (MQ=255)                                                                            gTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGACCAAAAACTTCTTCCGGCGTGGAGAGGCGGGTCATCAGATATTCAAACGCGCGGGCATTCTGAAAGGCGAATTTTAAGCGCCGGGCGATCAAACTTTTGGCAATAcc                                                                      <  2:918369/151‑1 (MQ=255)                                                                               cATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGACCAAAAACTTCTTCCGGCGTGGAGAGGCGGGTCCTCAGATATTTAAACGCGCGGGCATTCTGCAAAGCGACTTTTAAGCGACGGGAGATCAAAATTTTGGCCATACCTg                                                                    >  1:176832/1‑150 (MQ=255)                                                                                aTAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGACCAAAAACTTCTTCCGGCGTGGAGAGGCGGGTCATCAGATATTCAAACGCGCGGGCATTCTTAAAAGCGAATTTTAAGCGCCGGGCGATCAAACTTTTTGCAATaaccgg                                                                   >  2:95400/1‑145 (MQ=255)                                                                                  aGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGACCAAAAACTTCTTCCGGCGTGGAGAGGCGGGTCATCAGATATTCAAAAGCGCGGGGCTTCTGAAAGGCGAATTTTTAGCGCCGGGCGATCAAACTTTTGGGAATACCTGgg‑gg                                                               >  1:464841/1‑151 (MQ=255)                                                                                   gcgcCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGACCAAAAACTTCTTCCGGCGTGGAGAGGCGGGTCATCAGAAATTCAAAAGCGCGGGGATTCTGGAAGGGGAATTTTAAGCGCCGGGCGATCAAACTTTTTGCAATACCTGGC‑ggg                                                              >  2:1183957/1‑151 (MQ=255)                                                                                    cgcCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGACCAAAAACTTCTTCCGGCGTGGAGAGGCGGGTCCTCAGATATTCAAACGCGCGGGCATTCTGAAAGGCGAATTTTAAGCGCCGGGCGATTAAACTTTTGGCAATACCCGGC‑ggg                                                              >  2:661542/1‑150 (MQ=255)                                                                                       ccTTCATCTTTTAGCGCCTGAATAGAAAGGGGACCAAAAACTTCTTCCGGCGTGGAGAGGCGGGGCATCAGATATTCAAACGCGCGGGCATTCTGAAAAGCGAAATTTAAACGCCGGGCGGACAAAATTTTTGCAATAACCGGC‑ggggc                                                            >  1:707686/1‑147 (MQ=255)                                                                                           cATCTTTTAGCGCCTGAATAGAAAGGGGACCAAAAACTTCTTCCGGCGTGGAGAGGCGGGTCATCAGATATTCAAACGCGCGGGCATTCTGAAAGGCGAATTTTAAGCGCCGGGCGATCAAACTTTTGGCAATACCTGGC‑GGGCCAaggag                                                      <  1:1111548/151‑1 (MQ=255)                                                                                                  tAGCGCCTGAATAGAAAGGGGACCAAAAACTTCTTCCGGCGTGGAGAGGCGGGTCATCAGATATTCAAACGCGCGGGCATTCTGAAAGGCGAATTTTAAGCGCCGGGCGATCAAACTTTTGGCAATACCTGGC‑GGGCCAAGGAGGAaca                                                 <  1:569856/149‑1 (MQ=255)                                                                                                   aGCGCCTGAATAGAAAGGGGACCAAAAACTTCTTCCGGCGTGGAGAGGCGGGTCATCAGATATTCAAACGCGCGGGCATTCTGAAAGGCGAATTTTAAGCGCCGGGCGATCAAACTTTTGGCAATACCTGGC‑GGGCCAAGGAGGAacaca‑c                                             <  2:1282306/151‑1 (MQ=255)                                                                                                      gcCTGAATAGAAAGGGGACCAAAAACTTCTTCCGGCGTGGAGAGGCGGGTCATCAGATATTCAAACGCGCGGGCATTCTGAAAGGCGAATTTTAAGCGCCGGGCGATCAAACTTTTGGCAATACCTGGC‑GGGCCAAGGAGGAACACA‑Cttt                                          <  1:973613/151‑1 (MQ=255)                                                                                                          gAATAGAAAGGGGACCAAAAACTTCTTCCGGCGTGGAGAGGCGGGTCATCAGATATTCAAACGCGCGGGCATTCTGAAAGGCGAATTTTAAGCGCCGGGCGATCAAACTTTTGGCAATACCTGGC‑GGGCCAAGGAGGAACACA‑CTTTca                                        <  1:552670/149‑1 (MQ=255)                                                                                                                  aGGGGACCAAAAACTTCTTCCGGCGTGGAGAGGCGGGTCCTCAGGTATTTAAACCCGCGGGGATTCTGGAAAGGGTAATTTAAACGCCGGGCGGTCAAAAATTTGTGTATACCTGGC‑GGGCCAAGGAGGAACACA‑CTTTcccccc                                    >  2:1067527/1‑143 (MQ=255)                                                                                                                        ccAAAAACTTCTTCCGGCGTGGAGAGGCGGGTCATCAGATATTCAAACGCGCGGGCATTCTGAAAGGCGAATTTTAAGCGCCGGGCGATCAAACTTTTGGCAATACCTGGC‑GGGCCAAGGAGGAACACA‑CTTTCACCACTTAATGCCGCta                        <  1:1157765/151‑1 (MQ=255)                                                                                                                          aaaaaCTTCTTCCGGCGTGGAGAGGCGGGTCATCAGATATTCAAACGCGCGGGGCTTTTGAAAGGCGAATTTTAAGCGCCGGGCGATCAAAATTTTTGCAATACCTGGC‑GGGGCAAGGAGGAACACA‑CTTTCACCACCTAATGCCGCtaa                       >  2:607955/1‑150 (MQ=255)                                                                                                                          aaaaaCTTCTTCCGGCGTGGAGAGGCGGGTCATCAGATATTCAAACGCGCGGGCATTCTGAAAGGCGAATTTTAAGCGCCGGGCGATCAAACTTTTGGCAATACCTGGC‑GGGCCAAGGAGGAACACA‑CTTTCACCACTTAATGCCGCtaat                      <  1:1175975/151‑1 (MQ=255)                                                                                                                           aaaaCTTCTTCCGGCGTGGAGAGGCGGGTCATCAGAAATTCCAACGCGCGGGCCTTTTTAAAGGGGAATTTTAAACGCCGGGCGATCAAACATTTGGCCAAAACTGGC‑GGGGCAAGGGGGAAAACA‑CTTTCACCACTTAATTccccca                        >  1:1243963/1‑144 (MQ=255)                                                                                                                                   ttcCGGCGTGGAGAGGCGGGGGATCAGATATTCAAACGCGCGGGCATTCTTAAAGGGGAATTTTAAACGCCCGGGGATTAAAACTTTGGCAATACCTGGCGGGGCCAAGGAGGAA‑ACACCTTTCACC‑CTTAAAGGCCCCAATAAACACaa              >  1:73725/1‑150 (MQ=255)                                                                                                                                    tcCGGCGTGGAGAGGCGGGTCATCAGATATTCAAACGCGCGGGCATTCTGAAAGGCGAATTTTTAGCGCCGGGCGATCAAACTTTTGGGCATACCTGGC‑GGGGCAAGGAGGAACACA‑CTTTCACCACTTAAAGCCGCTAATAAAAACaa              >  2:1285225/1‑149 (MQ=255)                                                                                                                                    tcCGGCGTGGAGAGGCGGGTCATCAGATATTCAAACGCGCGGGCATTCTGAAAGGCGAATTTTAAGCGCCGGGCGATCAAACTTTTGGCAATACCTGGC‑GGGCCAAGGAGGAACACA‑CTTTCACCACTTAATGCCGCTAATAAACACaa              <  2:596786/149‑1 (MQ=255)                                                                                                                                    tcCGGCGTGGAGAGGCGGGTCATCAGATATTCAAACGCGCGGGCATTCTGAAAGGCGAATTTTAAGCGCCGGGCGATCAAACTTTTGGCAATACCTGGC‑GGGCCAAGGAGGAACACA‑CTTTCACCACTTAATGCCGCTAATAAACACaa              <  2:1145668/149‑1 (MQ=255)                                                                                                                                       ggCGTGGAGAGGCGGGTCATCAGATATTCAAACGCGCGGGCATTCTGAAAGGCGAATTTTAAGCGCCGGGCGATCAAACTTTTGGCAATACCTGGC‑GGGCCAAGGAGGAACACA‑CTTTCACCACTTAATGCCGCTAATAAACACAAGCGGa         <  1:81377/151‑1 (MQ=255)                                                                                                                                         cGTGGAGAGGCGGGTCATCAGATATTCAAACGCGCGGGCATTCTGAAAGGCGAATTTTAAGCGCCGGGCGATCAAACTTTTGGCAATACCTGGC‑GGGCCAAGGAGGAACACA‑CTTTCACCACTTAATGCCGCTAATAAACACAAGCGGATg       <  1:909841/151‑1 (MQ=255)                                                                                                                                            ggAGAGGCGGGTCATCAGATATTCAAACGCGCGGGCATTCTGAAAGGCGAATTTTAAGCGCCGGGCGATCAAACTTTTGGCAATACCTGGC‑GGGCCAAGGAGGAACACA‑CTTTCACCACTTAATGTCGCTAATAAAAACAAGCGGATGGCg    >  1:1009453/1‑151 (MQ=255)                                                                                                                                              agagGCGGGTCATCAGATATTCAAACGCGCGGGCATTCTGAAAGGCGAATTTTAAGCGCCGGGCGATCAAACTTTTGGCAATACCTGGC‑GGGCCAAGGAGGAACACA‑CTTTCACCACTTAATGCCGCTAATAAACACAAGCGGATGGCGTg  <  2:156599/151‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                      GAATTGCCGGGCCCGCTTTCCAGATCTCATCCAGAAAGACAATTTCTGCTTCCGGCAGGTAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGACCAAAAACTTCTTCCGGCGTGGAGAAGCGGGTCATCAGATATTCAAACGCGCGGGCATTCTGAAAGGCGAATTTTAAGCGCCGGGCGATCAAACTTTTGGCAATACCTGGC‑GGGCCAAGGAGGAACACA‑CTTTCACCACTTAATGCCGCTAATAAACACAAGCGGATGGCGTG  >  NC_000913/3930724‑3931015

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