Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,652,331	T→C	intergenic (+2346/+397)	ydfD → / ← ynfP	lysis protein/protein YnfP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,652,331	0	T	C	100.0%	151.3 / NA	43	intergenic (+2346/+397)	ydfD/ynfP	lysis protein/protein YnfP
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (22/21);  total (22/21)

TTATGGTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACA‑CAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGAC‑TGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAACC  >  NC_000913/1652085‑1652562                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                        ttATGGTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACA‑CAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGAC‑TGAGCAAAGAGCATATCATTGATGTTCACCTCa                                                                                                                                                                                                                                      <  2:723352/250‑1 (MQ=255)          gttgctggttATACTGGAGGTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTGTCACGGTACCAAAAACA‑CAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGAC‑TGAGCAAAGAGCATATCAGTGATGTTCACCTCAGAGAGTATgg                                                                                                                                                                                                                            <  2:613113/243‑1 (MQ=255)                gTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACA‑CAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGAC‑TGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGcagaac                                                                                                                                                                                                                      <  1:693567/251‑1 (MQ=255)                 tttATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACA‑CAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCCCCCCGGAGGCAGATGACATTAACGAGAC‑TGAGCCAAAAACATATCATTTATGTTCACCTCAGAGGGTATGGcagaac                                                                                                                                                                                                                      >  1:392825/1‑250 (MQ=255)                               cGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGAACAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACA‑CAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGAC‑TGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAAt                                                                                                                                                                                                        <  1:618697/250‑1 (MQ=255)                                  ccTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACA‑CAGGCGGGAAACCATCGGGGCATTCATCTTATTTAGCCAGACATCG‑CCGCTCTCCGGGGTCAGATGACATTAAAGAGAC‑TGAGCAACGAGCAAATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGCAAAACAAAAATGCAc                                                                                                                                                                                                    >  1:633112/1‑251 (MQ=255)                                   cgggTGAGCTGTGTTCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACA‑CAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGTATTCAGATGACATTAACGAGAC‑TGAGCAAAGAGCATATCATTGATGTTCACCTCACAGAGTATGGCAGAACAGAAAAACAAAAATGCAcc                                                                                                                                                                                                   <  1:222825/249‑1 (MQ=255)                                        gAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACA‑CAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGAC‑TGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTg                                                                                                                                                                                               <  1:562590/250‑1 (MQ=255)                                         aGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACA‑CAGGCGGGAACCAATAGGGGCATTTATCTTATTAAGCCAGAAATCG‑ACGCTCTACGGAGGCAGATGACATTAACGAGAC‑TGCGCAAAGAGCATCTCATTTTTGTTCACCTCAGAGAGGATGGCAGAACAGaaaaaaagaaataaagagggatc                                                                                                                                                                                             >  2:164622/1‑233 (MQ=255)                                          ggctggtgcaccgggATTGGATGACATAAATCTGGAAAATTGCACAATGATGGACCGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACA‑CAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGAC‑TGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCTCCTTTGttt                                                                                                                                                                                            <  2:562524/238‑1 (MQ=255)                                             gggtgaaCTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACA‑CAGGCGGGAACCAATCGGGGCATTCATCTTATTAAGCCAGCAAACG‑GCGCCCTCCGGAGGCAGATGACAATAACGA‑ACATGCGCAAAGAGCATATCATTTATTTTCACCTCAGAGAGTATGGCAGAACAACAAAACAAAAATTCCCCTTTTTTTTTc                                                                                                                                                                                         >  2:225965/7‑251 (MQ=255)                                                            ggaggaCATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACA‑CAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGAC‑TGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAg                                                                                                                                                                          <  1:239313/251‑1 (MQ=255)                                                                    tAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACA‑CAGGCGGGAACCAAAAGGGTCATTCATCTTTTTAAACCCGCAAACG‑ACGCCCTCCGGAGTCCAATTACAATAACGAGAA‑TGGACAAAGAGCCTCTCATTGTTGTTCTCCCCAGAGAGTATGGCCGGACAAGAAAACAAAAAGTCACCTTTTTTTTTCAACCTgtaaaggccacgggagt                                                                                                                                                                   >  2:314507/1‑234 (MQ=255)                                                                        aatgaagattGGAACGATGAGTATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGCAACAAAA‑CATCAGGCGGTAATCAATCGGGTCATTCAGCTTATTAAGCCAGTAATCG‑ACGCTCTCCGGAGTCAGATGACAGTAACGATAT‑GGAGGAAAGAGCATAGCAATGATGTGCATCTCAGAGAGTAGGGCAGAACAGAAAAACAAAAATGCACCTTGGTTTTTCAACCTGAAGGGTCAGCGAGAGTtaaaa                                                                                                                                                              <  2:53712/242‑1 (MQ=255)                                                                          tGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACA‑CAGGCGGGAACCCATCGGGTCATTCATCTTATTAAGCCAGCCATCG‑ACGCTCCCCCGAGTCAGATGACATTAACGAGAC‑TGAGCAAAGAGCATAACATTTATGTTCACCACAGAGAGTATGGCAGAACAGAAAAAAAAAAATGCACCTTTTTTTTTTAAACTGAAGGGTCAGCGGGAGAaaaaaa                                                                                                                                                             >  1:610194/1‑250 (MQ=255)                                                                            aaaaaTGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACA‑CAGGCGGTAACCAAACGGGTCAATCATCTTATTAAACCAGCAATCG‑AAGCTCTCCCGAGTCCAATGACATTAAAGAGAC‑TGAGCCAAGAAGATTACATTGTTGTTTACCTCAGAGAGGAAGGCAGAAAACAAAAACAAAAAAGCACCTTTTTGTTTCAACCTTAAGTGGCAGCGAGAGGaaaaaagt                                                                                                                                                           >  2:130988/1‑248 (MQ=255)                                                                                  ggAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACA‑CAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGAC‑TGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTg                                                                                                                                                     <  1:330425/250‑1 (MQ=255)                                                                                  ggAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACA‑CAGGCGGGAAACAAATGGGGCATTCATCTTATTAAGCCAGCAATAG‑ACGCTCTCCGGGGACAGATGACATTAAAGAGAA‑TGGGCAAAGAGCATATCATTGGTGTTCACCTCAGGGAGTATGGCAGAAAAGAAAAAACAAAATGCACCTTTGTTTTTTAACCTGAAGTGGCAGCGAGAGTAAAAAAATATGGTg                                                                                                                                                     >  1:112731/1‑250 (MQ=255)                                                                                  ggAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACA‑CAGGCGGGAAAAAATCGGGGCATTCATATTTTTAAGACAGCCAACG‑ACGCCCTCCGGAGTCAGATGACATTAACGAGAA‑TGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCGGAACAGAAAAACAAAAAAACACCTTTTTTTTTCCACCCGGAGTGTCCGCGAGAGGAAAAAATTATGGTgg                                                                                                                                                    >  1:276381/1‑250 (MQ=255)                                                                                        atgatgATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACA‑CAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGAC‑TGAGCAAAGAGCATATCATTGATGTTCACCTCCGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCCttt                                                                                                                                              >  1:438605/1‑251 (MQ=255)                                                                                            ggagcAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACA‑CAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGAC‑TGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCGTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTAcc                                                                                                                                          <  2:633112/247‑1 (MQ=255)                                                                                                 aGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACA‑CAGGCGGGCACCAATCTGGTCCTTCCTCTTATTAAGCCAGCAATCGGAC‑CGCTCCGGGGTCGGAGCACATTTACGAAAA‑TGAGCCAAACGCATATCATTGATGTTTACCTCAGAGAGAATGGAAAAAACGAAAAACAAAAATGCACCCTTGTTTTTCAACCTGTAGTGTTAGCGAGCGTAAAAAATTATGGGGACCAATTTAACGtttt                                                                                                                                     >  1:634603/1‑249 (MQ=255)                                                                                                  ggAGACATTTAGCAAAAGACCGTTTCACGGTACCAAAAACA‑CAGGCGGTAACCAAACGGGGCATTAATCTTTTTAAAACAGCAAACG‑ACACCCTCCCGAGGCAAAATAAAATAAAAAGAC‑TGAGCAAAAAGCAAATAATTGATGTTCCCCTCAAAGAGGATGGCAGAAAAGAAAAACAAAAATACACCCTTGTTTTTTAACCTGAAGGTTCAGCGAAAGAAAAAAAATATTGTGACCCTTTTAACGTTgtc                                                                                                                                    >  2:247132/1‑249 (MQ=255)                                                                                                            atcAAAAGACCGTTTCACGGTACCAAAAACA‑CAGGCGGGAACCAAGCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGAC‑TGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAgg                                                                                                                           <  2:500043/248‑1 (MQ=255)                                                                                                              cAAAAGACCGTTTCACGGTACCAAAAACA‑CAGGCGGGGACCAAACGGGTCACTCCTCTTATTAAACCAGCAATCG‑ACGCTCTCCGGGGGCAAATGACATCAAAGGAAC‑TTAGCAAAGAGCATAACATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAAAGCAACTTTGTTTTTTAACCTAAAGTGTCCGAGAGAGAGAAAAAATATGGTGACCATTTTACCGTTGACCCAATAAGGccg                                                                                                                        >  1:307296/1‑249 (MQ=255)                                                                                                               aaaaGACCGTTTCACGGTACCAAAAACA‑CAGGCGGGAACCAATCGGGGCATTCATCTTATTAAGCCAGCAATCG‑ACGCCCCCCGGAGTCAGAAGACATTAACGAGAC‑TGAGCAAAAAGCATATCATTGATGTTCACCTCAGAGAGTATGGCCGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTTTTGTGACCATTTTTCCGTTGACTCAATAAGGc                                                                                                                          >  1:276504/1‑248 (MQ=255)                                                                                                                     ccGTTTCACGGTACCAAAAACA‑CAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGAC‑TGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATTGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTggg                                                                                                                 <  2:634603/251‑1 (MQ=255)                                                                                                                      cGTTTCACGGTACCAAAAACA‑CAGGCGGGAAACAATCGGGGCATTAATCTTTTTAAGCCAACAATTG‑ACGCCCCCCGGAGGCAGATGACATTAACGAGAC‑TGAGCAAAAAGCATATCATTGATGTTCACCTCAGAGAATATGGCAGGACAGAAAAACAAAAATTCAACCTTGTTTTTTAACCTGAAGTGTCAGCGAGAGTAAAAAATTTTGGTGACAATTTTTCCGTTGACTCAAAAAGGGAGATGggggg                                                                                                                >  1:115760/1‑250 (MQ=255)                                                                                                                                 aCCAAAAACA‑CAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGAC‑TGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATa                                                                                                      <  1:315539/250‑1 (MQ=255)                                                                                                                                                              gTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGAC‑TGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCAATAAAACGTGCCGGACTCCGCCATCGAAAATc                                                                         >  1:307436/1‑251 (MQ=255)                                                                                                                                                                          attaAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAAAC‑TGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGCGAATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAATGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGAAGATGTGGGACGCAGCGATAAAACGCGCCGGACTCCGCCATCGAAAATCATATCAgctg                                                               >  2:591163/1‑246 (MQ=255)                                                                                                                                                                              aGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGAC‑TGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGGGCCGGACTCCGCCATCGAAAATCATATAAGTCGAGACAt                                                         >  2:624729/1‑251 (MQ=255)                                                                                                                                                                                    cAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGAC‑TGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTAt                                                   >  2:302638/1‑251 (MQ=255)                                                                                                                                                                                        cG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGAC‑TGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAAAGCACCTTTGTTTGTCATACTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACACCGCCATCGAAAATCATATCAGTCGAGACTTACTCATGCct                                               <  2:622394/251‑1 (MQ=255)                                                                                                                                                                                                  ccGGAGTCAGATGACAGTAACGAGAC‑TTAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGt                                         <  1:476573/248‑1 (MQ=255)                                                                                                                                                                                                    ggAGTCAGATGACATTAACGAGAC‑TGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTg                                    <  1:311075/251‑1 (MQ=255)                                                                                                                                                                                                     gggTCAGATGCCATTAATGAGAC‑TGAGCCAAGAGCATATCATTGATGTTCACCTCAGAGAGAATGGCAGAACAGAAAAACAAAAATGCAACTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCCTTTTACCGTTGACTCAATAAGGACGATGTGGGACGCAGCGAAAAAACGTGCCGGACTCAGCCATCGACAATCATATCAGTCGAGCCATACTTATGCCTGCTGGTCTCTg                                    >  2:724017/3‑250 (MQ=255)                                                                                                                                                                                                                cATTAACGAGAC‑TGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCt                        <  1:314507/251‑1 (MQ=255)                                                                                                                                                                                                                    aaCGAGAC‑TGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAAc                     <  1:473205/250‑1 (MQ=255)                                                                                                                                                                                                                       gagaC‑TGAGAAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTAGGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCgg                 <  2:547110/251‑1 (MQ=255)                                                                                                                                                                                                                         gaC‑TGAGCAAAGAGCATATCAGTGATGTTCCCCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCa               <  2:438605/251‑1 (MQ=255)                                                                                                                                                                                                                                  aaaGAGCATATCATTGATGTTCACCTCAGCGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGc       >  2:623250/1‑251 (MQ=255)                                                                                                                                                                                                                                       gCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAAcc  >  1:401651/1‑251 (MQ=255)                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                        TTATGGTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACA‑CAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGAC‑TGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAACC  >  NC_000913/1652085‑1652562

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