BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	677,120	G→A	A161A (GCG→GCA)	ybeS →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	677,120	0	G	A	99.6%	1051.8 / NA	272	A161A (GCG→GCA)	ybeS	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (0/0); major base A (147/124); minor base C (0/1); total (147/125)
Fisher's exact test for biased strand distribution p-value = 4.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 280 total aligned reads displayed.

GTTATTAGCAGAAAGATTGCGGTGGTTACAGGAGGAAAATACC‑GGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCGAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGATTTTTACCAAATGTTCGCTC  >  NC_000913/677030‑677208

gtTATTAGCAGAAAGATTGCGGTGGTTACAGGAGGAAAATACC‑GGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCaaa                                                                                        <  1:5809544/93‑1 (MQ=255)

 ttattaGCAGAAAGATTGCGGTGGTTACAGGAGGAAAATACC‑GGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAg                                                                                       <  2:5541652/93‑1 (MQ=255)

 ttattaGCAGAAAGATTGCGGTGGTTACAGGAGGAAAATACC‑GGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAg                                                                                       >  2:1331660/1‑93 (MQ=255)

  tattaGCAGAAAGATTGCGGTGGTTACAGGAGGAAAATACC‑GGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCaaacg                                                                                      >  2:5884835/1‑91 (MQ=255)

  tattaGCAGAAAGATTGCGGTGGTTACAGGAGGAAAATACC‑GGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAgg                                                                                      <  2:1131549/93‑1 (MQ=255)

  tattaGCAGAAAGATTGCGGTGGTTACAGGAGGAAAATACC‑GGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAgg                                                                                      >  2:4876801/1‑93 (MQ=255)

   attaGCAGAAAGATTGCGGTGGTTACAGGAGGAAAATACC‑GGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAggg                                                                                     <  1:6049594/93‑1 (MQ=255)

    ttaGCAGAAAGATTGCGGTGGTTACAGGAGGAAAATACC‑GGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAgggg                                                                                    <  1:1055936/93‑1 (MQ=255)

    ttaGCAGAAAGATTGCGGTGGTTACAGGAGGAAAATACC‑GGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAgggg                                                                                    <  1:3154198/93‑1 (MQ=255)

     taGCAGAAAGATTGCGGTGGTTACAGGAGGAAAATACC‑GGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGa                                                                                   <  2:6209497/93‑1 (MQ=255)

     taGCAGAAAGATTGCGGTGGTTACAGGAGGAAAATACC‑GGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGa                                                                                   >  2:3538494/1‑93 (MQ=255)

      aGCAGAAAGATTGCGGTGGTTACAGGAGGAAAATACC‑GGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGaa                                                                                  >  1:970033/1‑93 (MQ=255)

      aGCAGAAAGATTGCGGTGGTTACAGGAGGAAAATACC‑GGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGaa                                                                                  <  2:1088346/93‑1 (MQ=255)

      aGCAGAAAGATTGCGGTGGTTACAGGAGGAAAATACC‑GGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGaa                                                                                  <  2:1617152/93‑1 (MQ=255)

       gCAGAAAGATTGCGGTGGTTACAGGAGGAAAATACC‑GGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAAt                                                                                 <  1:2964635/93‑1 (MQ=255)

       gCAGAAAGATTGCGGTGGTTACAGGAGGAAAATACC‑GGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAAt                                                                                 <  2:1141426/93‑1 (MQ=255)

       gCAGAAAGATTGCGGTGGTTACAGGAGGAAAATACC‑GGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAAt                                                                                 <  2:3826778/93‑1 (MQ=255)

       gCAGAAAGATTGCGGTGGTTACAGGAGGAAAATACC‑GGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAAt                                                                                 <  2:3200352/93‑1 (MQ=255)

        cAGAAAGATTGCGGTGGTTACAGGAGGAAAATACC‑GGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGaaat                                                                                <  1:4116554/93‑3 (MQ=255)

            aaGATTGCGGTGGTTACAGGAGGAAAATACC‑GGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGtttt                                                                            >  1:5644058/1‑93 (MQ=255)

            aaGATTGCGGTGGTTACAGGAGGAAAATACC‑GGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGtttt                                                                            <  1:3860420/93‑1 (MQ=255)

             aGATTGCGGTGGTTACAGGAGGGAAATACC‑GGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCTAAGGGTAATGttttt                                                                           >  2:2320194/1‑93 (MQ=255)

             aGATTGCGGTGGTTACAGGAGGAAAATACC‑GGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGCATGttttt                                                                           >  1:2373439/1‑93 (MQ=255)

             aGATTGCGGTGGTTACAGGAGGAAAATACC‑GGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGttttt                                                                           <  1:1218354/93‑1 (MQ=255)

             aGATTGCGGTGGTTACAGGAGGAAAATACC‑GGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGttttt                                                                           <  1:2639030/93‑1 (MQ=255)

             aGATTGCGGTGGTTACAGGAGGAAAATACC‑GGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGttttt                                                                           >  2:4524442/1‑93 (MQ=255)

             aGATTGCGGTGGTTACAGGAGGAAAATACC‑GGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGttttt                                                                           >  2:2955203/1‑93 (MQ=255)

             aGATTGCGGTGGTTACAGGAGGAAAATACC‑GGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGttttt                                                                           <  2:5130719/93‑1 (MQ=255)

              gATTGCGGTGGTTACAGGAGGAAAATACC‑GGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTa                                                                          >  1:3296840/1‑93 (MQ=255)

               aTTGCGGTGGTTACAGGAGGAAAATACC‑GGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTaa                                                                         <  2:24155/93‑1 (MQ=255)

                ttGCGGTGGTTACAGGAGGAAAATACC‑GGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAAc                                                                        <  2:1041677/93‑1 (MQ=255)

                 tGCGGTGGTTACAGGAGGAAAATACC‑GGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACt                                                                       <  2:872292/93‑1 (MQ=255)

                   cGGTGGTTACAGGAGGAAAATACC‑GGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTc                                                                     >  2:5451286/1‑93 (MQ=255)

                   cGGTGGTTACAGGAGGAAAATACC‑GGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTc                                                                     <  2:6491535/93‑1 (MQ=255)

                     gtggtTACAGGAGGAAAATACC‑GGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCa                                                                   <  2:6505304/93‑1 (MQ=255)

                     gtggtTACAGGAGGAAAATACC‑GGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCa                                                                   <  2:5239623/93‑1 (MQ=255)

                     gtggtTACAGGAGGAAAATACC‑GGAGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCa                                                                   <  1:3966525/93‑1 (MQ=255)

                       ggtTACAGGAGGAAAATACC‑GGGGAGATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGa                                                                 <  2:2149288/93‑1 (MQ=255)

                       ggtTACAGGAGGAAAATACC‑GGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGa                                                                 >  2:5321939/1‑93 (MQ=255)

                       ggtTACAGGAGGAAAATACC‑GGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGa                                                                 >  2:4252939/1‑93 (MQ=255)

                         ttACAGGAGGAAAATACC‑GGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAGAGGGGAATGTTTTTAACTTCCAGAcc                                                               >  2:1541678/1‑93 (MQ=255)

                         ttACAGGAGGAAAATACC‑GGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGAcc                                                               >  2:3969378/1‑93 (MQ=255)

                          tACAGGAGGAAAATACC‑GGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCa                                                              >  1:5940890/1‑93 (MQ=255)

                           aCAGGAGGAAAATACC‑GGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCAt                                                             <  2:2061018/93‑1 (MQ=255)

                           aCAGGAGGAAAATACC‑GGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCAt                                                             >  1:1103022/1‑93 (MQ=255)

                           aCAGGAGGAAAATACC‑GGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCAt                                                             <  1:1944970/93‑1 (MQ=255)

                           aCAGGAGGAAAATACC‑GGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCAt                                                             <  1:5021147/93‑1 (MQ=255)

                             aggaggAAAATACC‑GGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTc                                                           >  1:4632111/1‑93 (MQ=255)

                               gaggaAAATACC‑GGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTg                                                         >  1:2885628/1‑93 (MQ=255)

                                aggaAAATACC‑GGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGc                                                        >  1:789410/1‑93 (MQ=255)

                                 ggaAAATACC‑GGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCa                                                       <  1:3570569/93‑1 (MQ=255)

                                   aaaaTACCGGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCAt                                                      <  2:5626722/93‑1 (MQ=255)

                                    aaaTACC‑GGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCt                                                    <  1:5590463/93‑1 (MQ=255)

                                    aaaTACC‑GGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCt                                                    <  1:3707179/93‑1 (MQ=255)

                                      aTACC‑GGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGc                                                  >  1:2747227/1‑93 (MQ=255)

                                      aTACC‑GGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGc                                                  >  2:3014858/1‑93 (MQ=255)

                                      aTACC‑GGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGGTTTTAACTTCCAGACCATTCTGCATCTGc                                                  >  1:4927121/1‑93 (MQ=255)

                                       tACC‑GGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGcc                                                 <  2:3397635/93‑1 (MQ=255)

                                         cc‑GGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCg                                               <  1:1709629/93‑1 (MQ=255)

                                         cc‑GGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCg                                               <  2:6173464/93‑1 (MQ=255)

                                          c‑GGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGt                                              <  1:2142457/93‑1 (MQ=255)

                                          c‑GGGGAAAGAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGt                                              >  1:6491662/1‑93 (MQ=255)

                                              ggAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGc                                           >  1:6452134/1‑93 (MQ=255)

                                               gAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGcc                                          >  1:6295716/1‑93 (MQ=255)

                                                aaaTAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCg                                         <  1:755690/93‑1 (MQ=255)

                                                aaaTAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCg                                         <  2:5729054/93‑1 (MQ=255)

                                                 aaTAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGt                                        >  2:6448765/1‑93 (MQ=255)

                                                  aTAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTg                                       >  2:2207920/1‑93 (MQ=255)

                                                  aTAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTg                                       >  1:2518558/1‑93 (MQ=255)

                                                     gACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCaa                                    <  1:5554617/93‑1 (MQ=255)

                                                       cGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCaaaa                                  <  1:5227471/93‑1 (MQ=255)

                                                       cGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCaaaa                                  >  1:3454210/1‑93 (MQ=255)

                                                         aagaagaaCTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATg                                >  2:1831279/1‑93 (MQ=255)

                                                          agaagaaCTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGa                               >  1:4284952/1‑93 (MQ=255)

                                                          agaagaaCTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGa                               <  1:1938240/93‑1 (MQ=255)

                                                          agaagaaCTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGa                               >  1:2250567/1‑93 (MQ=255)

                                                          agaagaaCTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGa                               >  1:3479779/1‑93 (MQ=255)

                                                            aagaaCTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGaca                             >  1:4479314/1‑93 (MQ=255)

                                                            aagaaCTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGaca                             >  1:2558509/1‑93 (MQ=255)

                                                               aaCTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCa                          >  2:93497/1‑93 (MQ=255)

                                                                aCTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCAt                         >  1:5372162/1‑93 (MQ=255)

                                                                 cTGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATc                        <  1:4356482/93‑1 (MQ=255)

                                                                  tGGAATCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCg                       >  2:714866/1‑93 (MQ=255)

                                                                      aTCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGAttt                   <  2:3882894/93‑1 (MQ=255)

                                                                       tCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGAtttt                  >  2:4696893/1‑93 (MQ=255)

                                                                       tCCTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGAtttt                  >  2:6068533/1‑93 (MQ=255)

                                                                       accTTTTTATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGAtttt                  <  1:6607266/92‑1 (MQ=255)

                                                                          tttttATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGATTTTTAc               >  2:4119631/1‑93 (MQ=255)

                                                                          tttttATATGCCATTGCAAAGGGGAATGTTTTTAAATTGCAGACCATTCTGCATCTGCCCGTTGCCATGGAAAATGACACCATCGATTTTTAc               >  2:580835/1‑93 (MQ=255)

                                                                           ttttATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGATTTTTaac              >  2:2616722/1‑91 (MQ=255)

                                                                           ttttATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGATTTTTAcc              >  2:6595720/1‑93 (MQ=255)

                                                                            tttATATGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGATTTTTACCa             <  1:4774420/93‑1 (MQ=255)

                                                                              tataTGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGATTTTTACCaaa           <  2:2742508/93‑1 (MQ=255)

                                                                               atatGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACGCCATCGATTTTTACCaact          >  2:37936/1‑91 (MQ=255)

                                                                                tatGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGATTTTTACCAAATg         <  1:1399249/93‑1 (MQ=255)

                                                                                 atGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGATTTTTACCAAATGt        <  2:1742616/93‑1 (MQ=255)

                                                                                 atGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGATTTTTACCAAATGt        <  2:4327495/93‑1 (MQ=255)

                                                                                 atGCCATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGATTTTTACCAAATGt        >  1:3007319/1‑93 (MQ=255)

                                                                                     cATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGATTTTTACCAAATGTTCGc    >  2:333525/1‑93 (MQ=255)

                                                                                     cATTGCAAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGATTTTTACCAAATGTTCGc    >  2:6151505/1‑93 (MQ=255)

GTTATTAGCAGAAAGATTGCGGTGGTTACAGGAGGAAAATACC‑GGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCGAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGATTTTTACCAAATGTTCGCTC  >  NC_000913/677030‑677208

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 24 ≤` `ATCG/ATCG` `< 30 ≤` `ATCG/ATCG` `< 35 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`