Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,995,269	C→A	A170E (GCG→GAG)	dapF →	diaminopimelate epimerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,995,269	0	C	A	100.0%	196.9 / NA	65	A170E (GCG→GAG)	dapF	diaminopimelate epimerase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (31/34);  total (31/34)

CAACTTCGAACCTTCCGCCGTGCCGTTTCGCGCTAACAAAGCGGAAAAGACCTATATT‑ATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGG  >  NC_000913/3995123‑3995413                                                                                                                                                    |                                                                                                                                                 cAACTTCGAACCTTCCGCCGTGCCGTTTCGCGCTAACAAAGCGGAAAAGACCTATATT‑ATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCgagg                                                                                                                                                <  2:922489/149‑1 (MQ=255)  aaCTTCGAACCTTCCGCCGTGCCGTTTCGCGCTAACAAAGCGGAAAAGACCTATATT‑ATGCGCGCCGCCAAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCgaggc                                                                                                                                               >  1:404710/1‑149 (MQ=255)    cgtctgaCCTTCCGGCGTGCCGTTGCGCGCTAACAAAGCGGAAAAGACCTAGATG‑AGGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGGGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGaggcgg                                                                                                                                             <  2:597385/143‑1 (MQ=255)      tCGAACCTCCCGCCGTGCCGTTTCGCGCTAACAAAGCGGAAACGACCTATATT‑ATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGAGGCGGTa                                                                                                                                           <  1:502971/149‑1 (MQ=255)       cGAACCTTCCGCCGTGCCGTTTCGCGCTAACAAAGCGGAAAAGACCTATATT‑ATGCGCGCCGCCGAGCAGACAATCGTATGCGTCGTGGTGGCGATGGGAAATCCGCATGGCGTGATTCAGGTCGATGATGTCGATACCGAGGCGGTAg                                                                                                                                          <  2:51337/149‑1 (MQ=255)            cTTCCGCCGTGCCGTTTCGCGCTAACAAAGCGGAAAAGACCTATATT‑ATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGAGGCGGTAGAAACg                                                                                                                                     >  1:448092/1‑149 (MQ=255)            cTTCCGCCGTGCCGTTTCGCGCTAACAAAGCGGAAAAGACCTATATT‑ATGCGCACCGCCGAGCAGACGATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCACGTCGATGCTGTCGATACCCAGGCGGTAGAAACg                                                                                                                                     >  1:829617/1‑149 (MQ=255)              tCCGCCGTGCCGTTTCGCGCTAACAAAGCGGAAAAGACCTAGATT‑ATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGAGGCGGTAGAAACGCt                                                                                                                                   <  2:181877/149‑1 (MQ=255)                cgccgTGCCGTTTCGCGCTAACAAAGCGGAAAAGACCTATCTT‑ATGCCCGCCGGCGAGCAGACTATCTTATGCGGCGTGGTGTCTATAGGACATCCGCATTGCGTGATTCAGGTCGATGATGTTGATACCGAGGCGGTAGAAACGCTTg                                                                                                                                 <  2:273134/149‑1 (MQ=255)                cgccgTGCCGTTTCGCGCTAACAAAGCGGAAAAGACCTAGATG‑ATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGAGGCGGTAGAAACGCTTg                                                                                                                                 <  1:157845/149‑1 (MQ=255)                cctcgtgccgtTTCGCGCTAACAAGGCGGAAAAGACCT‑TATTGATGCGCGCCGCCGAGCAGACAACCTTATGCCGCGTGGTGTCGATGGGAAATCCACATTGCGTCATTCAGGTCGATGATGTCGATACCGAGGCGGTAGAAACGCTTg                                                                                                                                 <  2:671962/146‑1 (MQ=255)                   cgggccgtTTCGCGCTAACAAAGCGGAAAAGACCTATATT‑ATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGAGGCGGTAGAAACGCTTGGTc                                                                                                                              <  2:433812/146‑1 (MQ=255)                      gccgtTTCGCGCTAACAAAGCGGAAAAGACCTATATT‑ATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGACATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGAGGCGGTAGAAACGCTTGGTCCTg                                                                                                                           <  1:793689/149‑1 (MQ=255)                        cgtTTCGCGCTAACAAAGCGGAAAAGACCTATATT‑ATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATGCAGGTCGATGATGTCGATACCGAGGCGGTAGAAACGCTTGGTCCTGtt                                                                                                                         >  2:1114028/1‑149 (MQ=255)                         gtTTCGCGCTAACAAAGCGGAAAAGACCTATATT‑ATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGAGGCGGTAGAAACGCTTGGTCCTGTTc                                                                                                                        <  2:364920/149‑1 (MQ=255)                           ttCCCGCTAACAAAGCGGAAAAGACCTATATT‑ATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGAGGCGGTAGAAACGCTTGGTCCTGTTCTg                                                                                                                      <  1:414366/149‑1 (MQ=255)                             cgcgCTAACAAAGCGGAAAAGACCTATATT‑ATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGAGGCGGTAGAAACGCTTGGTCCTGTTCTGGa                                                                                                                    >  1:1169730/1‑149 (MQ=255)                              gcgcTAACAAAGCGGAAAAGACCTATATT‑ATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGAGGCGGTAGAAACGCTTGGTCCTGTTCTGGaa                                                                                                                   <  2:265804/149‑1 (MQ=255)                                    aCAAAGCGGAAAAGACCTATATT‑ATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGAGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCAc                                                                                                             >  1:164300/1‑149 (MQ=255)                                    aCAAAGCGGAAAAGACCTATATT‑ATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGAGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCAc                                                                                                             <  1:337194/149‑1 (MQ=255)                                     cAAAGCGGAAAAGACCTATATT‑ATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGGTACCGAGGCGGTAGAAACGCTGGGTCCTGTTTTGGAAAGCCACg                                                                                                            >  1:113815/1‑149 (MQ=255)                                      aaaGCGGAAAAGACCTATATT‑ATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGAGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGa                                                                                                           >  1:1111836/1‑149 (MQ=255)                                      aaaGCGGAAAAGACCTATATT‑ATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGAGGCGGTAGAAACGCGTGGTCCTGTTCGGGAAAGCCACGa                                                                                                           >  2:856265/1‑149 (MQ=255)                                       aaGCGGAAAAGACCTATATT‑ATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGAGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAg                                                                                                          >  1:297300/1‑149 (MQ=255)                                        aGCGGAAAAGACCTATATT‑ATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGAGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAGc                                                                                                         <  1:362987/149‑1 (MQ=255)                                        aGCGGAAAAGACCTATAGT‑ATGCGCGTCGCCGAGCAGAAAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGAGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAGc                                                                                                         <  2:1105431/149‑1 (MQ=255)                                          cGGAAAAGACCTATATT‑ATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGAGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAGCGt                                                                                                       >  1:1052357/1‑149 (MQ=255)                                                 gACCTCGATT‑GTGCCCGACCCCGTGCAGACAATCTTATAAGTCGGAGTGTCGCTGGGAAATCCGCACAGCGCACGTCAGGACTAGGAGTTCGAGACCGAGGCGGTAGAAGCGCTTGGTCCGGTTCTGGAAAGCCACGAGCGTTTTCCgg                                                                                                <  2:558470/149‑1 (MQ=255)                                                    cTATATT‑ATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGAGGCGGTAGAAACGCTTGGTCCTGTTCTGGCAAGCCACGAGCGTTTGCCGGAgc                                                                                             >  2:62987/1‑149 (MQ=255)                                                     tataTT‑ATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGAGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAgcg                                                                                            >  1:805570/1‑149 (MQ=255)                                                      atatT‑ATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGGGTGGATGGGAAATAAGGATTGCGTGATTCAGGTCGATGATGTCGATACGGAGGAGGTAGAAAAGATTGGTACTGTTCTGGAAAGCAACGAGCGTTTTCCGGAccgc                                                                                           >  2:20792/1‑149 (MQ=255)                                                         tt‑atGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGAGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGcc                                                                                        >  1:1074157/1‑149 (MQ=255)                                                                 cgccgcCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGAGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCg                                                                                 >  1:1151720/1‑149 (MQ=255)                                                                    cgccgAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGAGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGtt                                                                              >  1:405104/1‑149 (MQ=255)                                                                      ccgAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGAGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGtttt                                                                            >  1:1351755/1‑149 (MQ=255)                                                                      ccgAGCAGACAATCTTATGCGGCGTGGTGTCGAAGGGAAATCCGCATCGCGTGAGTCAGGACGATGAGGTCGATACCGAGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAACATCGGtgtc                                                                            >  2:170165/1‑146 (MQ=255)                                                                         aGCAGACAATCTTATGCGCCGTGCTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGAGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATg                                                                         <  2:1240630/149‑1 (MQ=255)                                                                         aGCAGACAATCTTAGGCCGCGTGGGGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGAGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATg                                                                         <  1:575244/149‑1 (MQ=255)                                                                          gCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGAGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTATTCCGGTGCGCGCCAATATCGGTTTTATGc                                                                        <  1:714087/149‑1 (MQ=255)                                                                              aCAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGCTGCTGTCGATCCCGAAGCGGTAGAAACTCTAGGTCCTGTTCGGGAAAGCCTCGAGCGTGTTCCGGAGCTCGCCAATATCGGTTTTATGCAAGt                                                                    >  1:658677/1‑149 (MQ=255)                                                                               ctaTCTTGTGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGAGGCGGTAGAAACGCGTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTg                                                                   <  2:776969/147‑1 (MQ=255)                                                                                  tgtgaTGCGGCGTGGTTGCGAGGGGAACTCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGAGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGtt                                                                <  2:1086743/145‑1 (MQ=255)                                                                                  tCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGAGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGtt                                                                <  2:1194610/149‑1 (MQ=255)                                                                                    ttATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGAGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTACGGTGCGGTCCAATATCGGTTTTATGCAAGTGGTTaa                                                              <  2:297300/149‑1 (MQ=255)                                                                                       tGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGAGGCGGTAGAAACGCTTGGTCCTGTTCTGGAACGCCACGACCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAgcg                                                           <  2:686585/149‑1 (MQ=255)                                                                                          ggcgTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGAGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGa                                                        <  2:1169730/149‑1 (MQ=255)                                                                                          ggcgTGGTGGCGATGGGAAATCCGCATTGCGGGATTCAGGTCGATGATGTCGATACCGAGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGa                                                        <  2:9654/149‑1 (MQ=255)                                                                                          ggcgGGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGAGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGa                                                        <  1:782078/149‑1 (MQ=255)                                                                                             gtggtgTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGAGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCa                                                     <  1:1082654/149‑1 (MQ=255)                                                                                              tggtgTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGAGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCGATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCat                                                    >  1:298594/1‑149 (MQ=255)                                                                                                  gtCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGAGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTc                                                >  1:1114354/1‑149 (MQ=255)                                                                                                    cGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGAGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGt                                              <  2:448092/149‑1 (MQ=255)                                                                                                        gggAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGAGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGGGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTAc                                          >  2:485059/1‑149 (MQ=255)                                                                                                             aTCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGAGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGtt                                     >  1:154252/1‑149 (MQ=255)                                                                                                                 gCATTGCGTGATTCAGGTCGATGATGTCGATACCGAGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCACGTGGTTAAGCGCGAGCATAGTCGTTTACGCGTTTATg                                 >  1:187116/1‑149 (MQ=255)                                                                                                                  cATTGCGTGATTCAGGTCGATGATGTCGATACCGAGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGa                                >  1:449695/1‑149 (MQ=255)                                                                                                                  cATTGCGTGATTCAGGTCGATGATGTCGATACCGAGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGAGTTTATGa                                >  2:275897/1‑149 (MQ=255)                                                                                                                    tggggTGATTCAGGTCGATGAGGTCGATCCCGAGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCACGCGGAAAGTCGCGAGCATATTCGTTTACGCGTTTATGAGc                              <  2:296699/145‑1 (MQ=255)                                                                                                                        gTGATTCAGGTCGATGATGTCGATACCGAGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTgg                          >  2:613279/1‑149 (MQ=255)                                                                                                                         tGATTCAGGTCGATGATGTCGATACCGAGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTCTGAGCGTggg                         >  2:164278/1‑149 (MQ=255)                                                                                                                             tCAGGTCGATGATGTCGATACCGAGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAACATATTCGTTTACGCGTTTATGAGCGTGGGGCAg                     >  1:906937/1‑149 (MQ=255)                                                                                                                               aGGTCGGTGATGTCGGTCCCGAGGCCGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTGAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGga                   <  2:820314/149‑1 (MQ=255)                                                                                                                                   cGATGATGTCGATACCGAGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGaaa               <  2:1111836/149‑1 (MQ=255)                                                                                                                                      tgatgTCGATACCGAGGCGGTAGAAACGCTTGGTCCCGTTCGGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAAccc            <  1:34905/149‑1 (MQ=255)                                                                                                                                       gatgTCGATACCGAGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCa           <  2:164300/149‑1 (MQ=255)                                                                                                                                               tACCGAGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCg   >  2:820830/1‑149 (MQ=255)                                                                                                                                                aCCGAGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCgg  <  2:805570/149‑1 (MQ=255)                                                                                                                                                aCCGAGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCgg  >  2:811837/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                 CAACTTCGAACCTTCCGCCGTGCCGTTTCGCGCTAACAAAGCGGAAAAGACCTATATT‑ATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGG  >  NC_000913/3995123‑3995413

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